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ABSTRACT 


The  effects  of  the  zonal  and  meridional  superposition  of  finite  amplitude  upper- 
and  lower-level  perturbations  are  investigated  using  a  two-level  analytical  model  and 
the  Navy  Operational  Regional  Atmospheric  Prediction  System  (NORAPS). 
Superposition  of  a  jet  streak,  aloft  and  a  finite  amplitude  low-level  perturbation  is  also 
investigated  using  NORAPS.  The  initial  along-stream  phase  difference  is  crucial  to  the 
subsequent  intensification  of  the  low-level  perturbation.  A  westward  tilt  with  height  of 
slightly  less  than  a  quarter  wavelength  is  most  favorable.  The  initial  meridional 
superposition  only  modifies  the  growth  rate  that  is  established  by  the  optimum  zonal 
phase  difference.  The  strength  and  location  of  low-level  baroclinity  and  upper-level 
shear  vorticity,  which  depend  on  the  structure  of  the  mean  fiow,  are  important  factors 
in  determining  the  most  favorable  cross-track  tilt.  In  the  jet  streak  case  with  strong 
upper-level  forcing,  the  lower  perturbation  quickly  moves  to  the  most  favored 
meridional  position  under  tne  left  front  quadrant. 

An  east  coast  cyclogenesis  event  during  19-21  January  1986  is  used  to  iHustrate 
the  effects  of  superposition.  Despite  strong  upper-level  forcing,  rapid  cyclogenesis  at 
the  surface  did  not  occur,  because  the  optimum  superposition  between  the  upper-level 
vorticity  center  and  the  surface  cyclone  was  not  achieved.  Modifications  to  the  initial 
conditions  along  the  Atlantic  coast  produce  a  simulation  with  a  more  favorable 
superposition,  which  results  in  more  rapid  intensification. 
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I.  INTRODUCTION 

Maritime  extratropical  cyclogenesis  events  are  among  the  most  pqwerful  and 
destructive  weather  systems  outside  the  tropics.  The  western  oceans  are  prime  areas 
for  cyclogenesis,  especially  during  the  cold  season  (Miller,  1946;  U.S.  Nav\%  1955;  many 
others).  Because  most  of  these  systems  intensify  steadily,  they  are  relatively  well- 
predicted.  However,  about  20  percent  of  the  cyclogenesis  events  include  a  short 
(12-36  h)  period  of  extremely  rapid  intensification.  This  intensification,  which  is  often 
poorly  forecast,  results  in  a  rapid  rise  in  wind  velocity  and  sea  heights.  Improved 
understanding  of  these  "rapid  deepeners"  will  be  beneficial  to  coastal  forecasters,  the 
commercial  shipping  and  fishing  industries,  and  to  the  U.S.  Navy. 

For  tactical  and  strategic  reasons,  it  is  often  necessary*  or  desirable  for  naval  task 
groups  to  operate  in  regions  where  winter  storms  are  common.  Accurate  and  timely 
forecasts  of  the  onset  and  intensity  of  these  events  are  required  so  that  operational 
commanders  can'  make  decisions  regarding  ship  and  aircraft  deployments  and  select 
appropriate  missions.  Although  hazards  to  mariners  from  high  winds  and  seas  are  of 
primar>'  concern,  many  decisions  are  influenced  by  direct  and  indirect  effects  of 
cyclones.  Important  tactical  variables  that  are  greatly  influenced  by  the  passage  of  a 
strong  storm  include  radar  or  electro-optical  sensor  performance  and  the  acoustic 
characteristics  of  the  upper  ocean. 

Unfortunately,  the  storms  that  are  most  dangerous  and  have  the  greatest  impact 
on  naval  operations  are  the  least  well  understood  and  the  least  reliably  forecast. 
Numerical  weather  prediction  models  have  a  poor  track  record  in  predicting  the 
magnitude  and  timing  of  intensification  in  these  rapidly  deepening  systems  (Sanders 
and  Gyakum,  1980;  Bosart,  1981).  To  better  forecast  these  lows,  it  is  necessary  to 
achieve  a  better  understanding  of  the  interactions  between  the  synoptic-scale, 
mesoscale  and  boundary  layer  processes  that  are  required  for  rapid  intensification. 
Sanders  and  Gyakum  (1980)  arbitrarily  chose  a  deepening  rate  of  1  mb  h'^  for  24  h 
(adjusted  for  latitude  by  the  factor  sin  (p/sin  60°)  as  the  division  between  ordinar\- 
cyclogenesis  and  rapid  deepening.  .A.t  40°  N,  this  is  approximately  equivalent  to  IS  mb 
in  24  h.  This  deepening  rate  has  been  generally  accepted  as  the  lower  end  of  the  rapid 
deepener   scale.     The   most   rapidly  intensifying   storm  during   the   three-year   period 
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studied  by  Sanders  and  Gyakum  had  a  deepening  rate  more  than  three  times  the 
minimum  value. 

A  general  hypothesis  of  rapid  cyclogenesis  has  been  defined  for  the  Experiment 
on  Rapidly  Intensifying  Cyclones  over  the  Atlantic  (ERICA),  which  is  part  of  a  five- 
year  program  of  basic  research  sponsored  by  the  Office  of  Naval  Research  (Hadlock, 
1987).  In  condensed  form,  the  hypothesis  is  that  rapid  cyclogenesis  in  the  western 
Atlantic  depends  on  the  proper  superposition  of  strong  positive  vorticity  advection 
(PVA)  aloft  over  a  pre-existing  low- level  cyclone  in  a  favorable  environment.  There 
are  several  contributing  factors  that  are  related  to  this  hypothesis.  Some  of  these 
factors  may  also  be  present  in  cases  which  are  not  rapid  deepeners.    . 

Jet-level  processes  may  enhance  cyclogenesis  in  several  ways.  Maximum 
divergence  aloft,  which  encourages  upward  vertical  motion  and  convection,  is  expected 
when  the  low-level  center  is  under  the  left-exit  or  right-entrance  quadrants  of  a  jet 
streak.  In  the  exit  region,  the  lower  branch  of  the  ageostrophic  circulation  assists  low- 
level  convergence  into  the  surface  center  and  provides  a  source  of  moisture  for  deep 
convection,  which  encourages  coupling  of  the  upper  and  lower  systems  and  increases 
the  baroclinic  transfer  of  eddy  available  potential  energy  (EAPE)  to  eddy  kinetic  energy 
(EKE). 

Jet  streaks  also  advect  positive  vorticity  into  the  upper  levels  above  the  surface 
system.  Even  with  straight  flow  at  the  jet  level,  an  intense  jet  with  strong  cyclonic 
shear  on  the  poleward  side  can  provide  significant  PVA.  The  magnitude  of  the  PVA 
will  be  increased  even  more  if  the  jet  streak  is  propagating  through  the  base  of  a  short 
wave  trough  or  from  trough  to  ridge. 

A  pre-existing  maritime  cyclone  provides  a  starting  point  for  rapid  deepening. 
Presence  of  a  pre-existing  perturbation  reduces  the  time  required  for  the  growth  of  a 
small-amplitude  perturbation  to  a  significant  amplitude  cyclone.  The  combination  of 
PVA  aloft,  a  surface  low,  and  a  favorable  environment  must  be  maintained  for  some 
time  if  the  new  system  is  to  become  sufficiently  organized  to  achieve  rapid  deepening. 
The  surface  low  circulation  also  organizes  low-level  moisture  over  a  relatively  deep 
layer. 

For  rapid  cyclogenesis,  the  low-level  environment  must  be  highly  baroclinic  and 
have  low  static  stability.  Such  an  environment  is  created  when  a  deep  layer  of  cold  air 
fiows  over  warm  surfaces,  especially  in  regions  with  large  thermal  gradients.  These 
conditions  are  most  frequently  satisfied  following  winter  season  cold  outbreaks  over 
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warm  oceans.  Strong  surface  temperature  gradients  can  be  found  at  the  land-sea 
interface  or  over  the  "north  wall"  of  the  Gulf  Stream.  Cold-air  damming  along  the 
Atlantic  seaboard  can  greatly  intensify  the  low-level  baroclinic  zone  along  the  coast. 
Synoptic  patterns  that  cause  cold  outbreaks  over  the  western  Atlantic  are  also 
favorable  for  subsequent  cold-air  damming  events.  The  optimal  environment  for 
explosive  deepening  in  the  western  Atlantic  is  hypothesized  to  be  one  in  which  cold-air 
damming  intensifies  the  coastal  baroclinic  zone,  while  cold  air  flows  steadily 
southeastward  across  a  strong  sea  surface  temperature  gradient  at  the  edge  of  the  Gulf 
Stream.  Continued  flow  of  cold  air  over  the  warm  water  will  maintain  low-level 
baroclinity  and  lead  to  large  air-sea  fluxes  of  sensible  and  latent  heat.  Thus,  these 
processes  may  work  together  to  broaden  and  strengthen  the  rapid  deepening  domain. 

In  summary,  many  factors  can  play  important  roles  in  rapidly  intensifying 
cyclones.  This  study  will  concentrate  on  the  importance  of  the  "proper  superposition" 
between  the  source  of  PVA  aloft  and  the  existing  surface  low.  For  this  research, 
rapidly  intensifying  cyclones  ("rapid  deepeners")  will  be  deflned  to  have  intensification 
rates  equal  to  or  greater  than  the  Sanders  and  Gyakum  definition. 

In  this  study,  finite  amplitude  ■perturbations  at  lower  and  upper  levels  will  be 
placed  in  a  variety  of  superpositions  in  two  atmospheric  models-a  two-level  analytical 
model  and  a  channel  version  of  the  'Sa.vy  Operational  Regional  Atmospheric  Prediction 
System  (NORAPS).  The  emphasis-  is  on  the  development  of  finiie-ampliiude 
disturbances  that  are  placed  in  a  baroclinically  unstable  mean  flow.  The  perturbations 
are  inserted  at  400  mb  and  SOO  mb  in  both  models  and  the  subsequent  development  at 
both  levels  is  described.  As  only  one  perturbation  wavelength  is  present,  this  approach 
differs  from  classical  experiments  in  which  infinitesimal  perturbations  with  a  spectrum 
of  wavelengths  are  allowed  to  grow  in  a  baroclinically  unstable  mean  flow.  In  this 
study,  the  perturbation  wavelengtii  represents  the  most  unstable  mode  that  eventually 
would  dominate  in  the  classical  experiments  for  the  mean  conditions  in  each  case. 
Different  combinations  of  perturbation  amplitudes  and  positions  are  tested  to 
demonstrate  the  influence  of  superposition  on  cyclogenesis  rates.  One  immediate  goal 
of  the  study  is  to  provide  further  evidence  in  support  of  the  superposition  aspects  ot^ 
the  ERICA  hypothesis  stated  above.  A  better  understanding  of  the  superposition 
aspects  will  ultimately  improve  forecasting  of  these  dangerous  maritime  cyclones. 
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A  review  of  literature  pertaining  to  rapid  cyclogenesis  is  presented  in  Chapter  II. 
The  idealized  model  results  are  in  Chapters  III  (two-level  model)  and  IV  (NORAPS 
channel  model).  In  Chapter  V,  the  conclusions  from  the  model  studies  are  tested  in  a 
case  study  from  the  Genesis  of  Atlantic  Lows  Experiment  (GALE).  Chapter  VI  is  a 
summarv  of  the  conclusions  from  the  research. 
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II.  LITERATURE  REVIEW 

A  number  of  papers  pertaining  to  rapidly  intensifying  cyclones  have  been 
published  since  the  landmark,  study  of  Sanders  and  Gyakum  (1980),  which  first 
separated  rapid  deepening  from  ordinary-  cyclogenesis.  Some  of  the  relevant  works  are 
reviewed  here,  along  with  some  older  references  that  are  also  applicable  to  the  study  of 
rapid  cyclogenesis.  The  references  will  be  addressed  as  they  apply  to  the  key 
requirements  of  the  hypothesis: 

•  upper-level  PVA; 

•  pre-existing  low-level  cyclone;  and 

•  favorable  environment. 

A.       UPPER-LEVEL  PROCESSES 

The  siiperposition  hypothesis  requires  that  a  source  of  positive  vorticity  aloft  be 
properly  positioned  with  respect  to  an  existing  surface  cyclone  that  has  a  horizontal 
scale  that  closely  matches  the  scale  of  maximum  baroclinic  iftstability  of  the  jet-level 
feature.  Both  jet  streaks  and  mobile  short  waves  may  serve  as  sources  of  upper-level 
PVA. 

The  Sanders  and  Gyakum  (1980)  study  included  267  cases  of  rapid  deepening 
that  occurred  in  the  Northern  Hemisphere  during  the  three  cold  seasons  from  autumn 
1976  to  spring  1979.  The  average  position  of  the  surface  low  was  approximately  400  n 
mi  east-northeast  of  a  mobile  500  mb  trough.  This  mobile  upper  trough  was  found  in 
the  southwest  quadrant  of  the  surface  cyclone  in  more  than  75  percent  of  the  cases. 
More  generally,  the  nearest  500  mb  trough  was  in  the  western  semi-circle  in  over  90 
percent  of  the  events.  In  a  study  of  48  cases  of  rapid  deepening  in  the  western  North 
Atlantic  Ocean,  Sanders  (1986)  concluded  that  rapid  cyclogenesis  is  fundamentally  a 
baroclinic  process,  and  that  a  prominent  absolute  vorticity  (^^)  maximum  at  500  mb 
contributes  significantly  to  rapid  deepening. 

This  scenario  of  a  propagating  upper  trough  overtaking  a  deepening  surface 
cyclone  is  consistent  with  the  classical  Petterssen  development  equation.  Petterssen 
(1956)  defined  cyclone  intensification  as  the  change  in  geostrophic  surface  vorticity, 
and  related  it  to  the  change  in  upper-level  vorticity  and  the  change  in  thickness 
between  the   upper-level  pressure   surface  and    1000  mb.     Thus,   it   is  expected   that 
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surface  cyclogenesis  is  associated  with  the  region  of  positive  vorticity  advection  (PVA) 
in  advance  of  the  mobile  500  mb  trough. 

Cyclonic  shear  vorticity  on  the  poleward  side  of  the  jet  enhances  PVA  above 
perturbations  on  the  cold  side  of  the  surface  baroclinic  zone.  This  may  explain  Sanders 
and  Gyakum's  (1980)  observation  that  the  maximum  in  explosive  cyclogenesis  tends  to 
occur  just  north  of  the  strongest  tropospheric  baroclinic  zone.  Based  on  data  from 
nine  winters,  Blackmon  et  al.  (1977)  also  found  that  the  band  of  maximum  baroclinic 
eddy  activity  was  below  the  cyclonic  shear  zone  on  the  poleward  side  of  the  strongest 
jet  streams. 

In  addition  to  upper-level  short  waves,  alongstream  variability  in  the  intensity  of 
the  jet  is  important  in  cyclogenesis.  The  polar  front  jet  is  characterized  by  propagating 
cores  with  higher  speeds.  These  isotach  maxima  were  defined  by  Palmen  and  Newton 
(1969)  as  jet  streaks.  The  jet  streaks  propagate  eastward  much  more  slowly  than  the 
winds  at  jet  level.  Therefore,  parcels  transiting  the  streak  accelerate  in  the  entrance 
region  and  then  decelerate  in  the  exit  region.  The  thermally  direct  circulation  in  the 
entrance  region,  with  divergence  aloft  in  the  right  quadrant  and  convergence  in  the  left, 
'  converts  available  potential  energy  to  kinetic  energy.  As  the  parcel  exits  the  streak. 
kinetic' energy  is  converted  to  potential  energy  by  the  thermally  indirect  circulation.  In 
the  exit  region,  there  is  divergence  aloft  in  the  left  quadrant  and  convergence  in  the 
right. 

Uccellini  and  Johnson  (1979)  demonstrated  the  importance  of  jet  streaks  in  the 
development  of  severe  convection.  The  upper  tropospheric  jet  streak  advects  cool,  dr\' 
air  into  the  upper  part  of  the  cyclone  environment,  enhances  upper-level  divergence, 
and  transports  sensible  heat  downstream  from  the  convective  region.  Additionally,  the 
lower  branch  of  the  ageostrophic  transverse  circulation  below  the  exit  region  of  the 
upper  jet  streak  can  induce  a  low-level  jet  (LLJ)  that  imports  warm,  moist  air  into  the 
convective  region.  This  low-level  jet,  termed  the  "conveyor  belt"  by  Harrold  (1973), 
can  provide  a  connection  between  a  source  of  warm,  moist  air  (such  as  the  warm 
Atlantic  Ocean)  and  the  convective  region.  The  combined  result  of  the  upper  jet  streak 
and  the  LLJ  is  to  create  a  region  of  convective  instability  with  low-level  convergence 
and  upper-level  divergence. 

Shapiro  (1983)  suggested  that  proper  alignment  of  the  upper  and  lower  jets  can 
facilitate  deep  convection  and  couple  the  low-level  system  to  the  jet-level  processes.  As 
an  upper-level  jet  overtakes  a  surface  front,  the  pattern  shown  in  Fig.  2.1a  is  formed. 
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The  LLJ  advects  warm,  moist  air  into  the  area,  while  the  upper  jet  provides  cooler, 
drier  air  aloft,  and  makes  the  lapse  rate  convectively  unstable.  The  ascent  branches  of 
the  upper  and  lower  secondar\'  circulations  are  vertically  aligned  (Fig.  2.1b).  which 
produces  a  deep  convective  plume  just  ahead  of  the  surface  cold  front. 
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Figure  2.1.   Vertically  coupled  upper  and  lower  jet/front  systems  and  associated 
secondan.'  circulations,   (a)  Jet  streak  exit  situated  above  the  surface  front  and 
low-levefjet.   (b)  Cross-section  along  the  line  BB'  of  (a).   (Shapiro,  1983) 


Calland  fl9S3)  studied  a  rapid  cyclogenesis  event  that  occurred  in  the  western 
North  Pacific  during  12-16  Januarv'  1979.  Because  cyclogenesis  occurred  under  almost 
straight  500  mb  flow,  this  was  a  case  in  which  upper-level  support  for  development  (in 
terms  of  PVA)  was  expected  to  be  mJnimal.  However,  vorticity  budget  analyses 
showed  that  horizontal  vorticity  advection  was  important,  even  without  a  significant 
short  wave  trough  aloft.    .A.  large  eddy  vorticity  fiux  at  upper  levels  was  associated  with 
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advection  of  cyclonic  shear  vorticity  by  an  intense  jet  streak.  Rapid  intensification 
occurred  when  the  surface  low  moved  under  the  upper  divergence  associated  with  the 
left-exit  quadrant. 

An  Atlantic  polar  low  that  also  developed  under  relatively  straight  flow  at  500 
mb  was  studied  by  Cook  (1983).  The  incipient  low  formed  south  of  Iceland  on  26 
January  1979  and  deepened  rapidly  as  it  moved  southeastward.  As  in  Calland's  case, 
upper-level  eddy  vorticity  transport  into  the  budget  volume  was  present  despite  the 
lack  of  an  appreciable  short  wave  aloft.  The  rapid  deepening  coincided  with  the 
movement  of  the  divergent  left-exit  quadrant  of  a  propagating  jet  streak  over  the 
surface  low. 

The  Presidents'  Day  storm  of  18-20  February  1979  has  been  extensively  studied 
by  Bosart  (1981),  Bosart  and  Lin.  (1984)  and  Uccellini  et  al.  (1984,  1985).  The  surface 
low  formed  in  an  inverted  trough  along  the  South  Carolina  coast  and  moved 
northward  along  a  coastal  front.  The  cyclone  was  influenced  by  two  upper-level  jet 
streaks.  Genesis  of  the  low  in  the  coastal  trough  was  accompanied  by  an  intensifying 
subtropical  jet  (STJ)  streak  moving  from  Kentucky  to  the  Virginia/North  Carolina 
coast  (Uccellini  et  ai,  1984).  A  strongly  ageostrophic  low-level  jet  over  the  Carolinas 
advected  moisture  northward  and  upward,  and  contributed  to  a  region  of  moderate  to 
heavy  snow  in  the  western  Carolinas.  Rapid  deepening  in  the  Presidents'  Day  storm 
occurred  later,  when  a  polar  front  jet  (PFJ)  streak  and  associated  short  wave  trough 
moved  eastward  over  the  coastal  low. 

Uccellini  (1986)  argues  that  upper-level  baroclinic  processes  also  contributed  to 
the  rapid  intensification  of  the  September  1978  QE  II  storm  (so  named  due  to  the 
damage  suffered  by  the  oceanliner  Queen  Elizabeth  II).  Gyakum  ( 1983a. b)  and  Anthes 
e[  al.  (1983)  attribute  the  QE  II  deepening  to  diabatic  processes,  but  Uccellini  showed 
evidence  that  an  intensifying  short  wave  trough/jet  streak  system  upstream  o{  the 
surface  low  could  have  contributed  to  the  rapid  development. 

Numerical  studies  have  also  emphasized  the  role  of  upper-level  jets  in  rapid 
cyclogenesis.  Mullen  (1979)  composited  NMC  analyses  for  the  genesis  periods  of  22 
eastern  North  Pacific  polar  lows.  Hodur  (1984)  used  this  composite  jet,  with  maximum 
winds  approximately  twice  the  values  found  in  a  mean  winter  jet,  to  study 
intensification  of  Pacific  polar  lows.  He  found  that  the  linear  growth  rates  associated 
with  the  mean  winter  jet  were  approximately  0.54  d'^  while  growth  rates  associated 
with  the  genesis  jet  were  1.2  d''.    Hayes  (1985)  suggests  that  intense  jets  crossing  the 
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coast  from  continent  to  sea  are  conducive  to  development  of  surface  lows  with 
horizontal  scales  typical  of  maritime  cyclones.  Sinclair  and  Elsberry  (1986)  describe  a 
case  of  cyclogenesis  in  a  numerical  simulation  that  occurred  on  the  anticyclonic  side  of 
the  jet.  Nevertheless,  jet-level  processes  appear  to  have  been  important,  since  the 
development  took  place  under  the  divergence  associated  with  the  right-entrance 
quadrant  of  a  propagating  jet  streak. 

The  importance  of  jet  processes  was  also  highlighted  by  Ebersole  (1984)  and  Peak 
ei  al.  (1985).  Forecasts  made  from  all  available  data  were  compared  to  forecasts  made 
without  satellite-retrieved  information.  For  the  Calland  (1983)  case,  the  "no  satellite" 
forecasts  were  surprisingly  superior.  The  difference  was  traced  to  discrepancies  in  the 
position  and  strength  of  the  upper-level  jet  stream.  This  sensitivity  to  the  strength  of 
the  jet  stream  reinforces  the  idea  that  jet  streaks  play  a  crucial  role  in  explosive 
cyclogenesis. 

B.       THE  PRE-EXISTING  CYCLONE 

The  second  key  factor  in  the  hypothesis  is  the  presence  of  a  pre-existing  maritime 
cyclone  at  low  levels.  In  a  composite  study  of  rapidly  intensifying  Atlantic  cyclones, 
Rogers  and  Bosart  (1986)  showed  that  rapid  deepeners  typically  evolve  from' shallow 
lower-tropospheric  systems.  The  presence  of  this  low  between  the  Appalachians  and 
the  Gulf  Stream  allows  rapid  cyclogenesis  to  begin  immediately  when  a  baroclinically 
unstable  jet  streak  of  proper  scale  is  superposed. 

In  the  absence  of  a  pre-existing  surface  low,  cyclogenesis  will  still  occur  when  a 
baroclinically  unstable  jet  streak  crosses  the  coast  to  a  favorable  low-level  environment. 
This  unstable  basic  state  encourages  exponential  growth  of  small  perturbations,  with 
the  wavelength  of  maximum  growth  being  shorter  over  the  ocean  than  over  contments. 
The  decrease  in  the  most  unstable  wavelength  is  due  tO"  reduced  friction,  absence  of 
orography,  reduced  static  stability,  and  fluxes  of  sensible  and  latent  heat  (Tokioka. 
1973;  Gambo,  1976;  Staley  and  Gall,  1977;  Sandgathe,  1981). 

This  exponential  growth  from  perturbation  to  incipient  low  requires  a  finite 
amount  of  time.  Since  the  "rapid  deepening  zone"  has  a  limited  east-west  extent,  the 
favorable  superposition  of  upper  support  and  conducive  environment  has  a  limited 
lifetime  for  each  case.  Without  a  pre-existing  system,  the  available  time  within  the 
zone  is  only  sufficient  to  develop  a  new  finite  amplitude  low  from  the  small  initial 
disturbance.    The  new  low  has  then  already  propagated  out  of  the  rapid  deepening 
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domain  before  rapid  cyclogenesis  can  occur.  Therefore,  a  pre-existing  maritime  cyclone 
that  is  in  the  proper  position  is  a  prerequisite  for  explosive  deepening. 

Two  winter  synoptic  patterns  were  identified  by  Miller  (1946).  Type  A  cyclones 
form  as  a  wave  on  a  cold  front  that  divides  warm,  moist  air  over  the  ocean  from  cold, 
dry  air  associated  with  a  large  anticyclone  centered  over  the  Midwest.  These  cyclones 
typically  form  over  the  ocean  and  move  northeastward  along  the  existing  low-level 
baroclinic  zone. 

Miller's  Type  B  cyclogenesis  takes  place  southeast  of  an  older  low  situated  in  the 
Midwest  or  Great  Lakes  region.  The  new  center  forms  near  the  coast,  often  along  the 
warm  front  of  the  old  system.  A  wedge  of  cold  continental  air  propagates  southward 
along  the  eastern  slope  of  the  Appalachians  from  New  England  and  retards  the  warm 
front  over  the  coastal  plain.  This  cold-air  damming  against  the  Appalachians 
(Richwien,  1980)  does  not  slow  the  northward  movement  of  the  front  over  the  ocean, 
and  the  front  becomes  oriented  NE-SW  along  the  land-sea  boundar>'.  The  secondary- 
cyclone  can  then  form  as  a  wave  on  the  convoluted  warm  front  to  the  south,  or  along 
the  coastal  front.  Kocin  and  Uccellini  (l985a,b)  studied  a  24-year  sample  of  18  east 
coast  snowstorms.  Fifteen  of  these  cases  eventually  qualified  as  .rapid  deepeners.  In 
nearly  every  case,  cold-air  damming  and  coastal  frontogenesis  occurred. 

Orographic  processes  can  also  favor  cyclogenesis  east  of  the  Appalachians.  Cold 
air  damming  associated  with  the  Appalachians  is  often  a  precursor  of  Type  B 
cyclogenesis.  The  mountains  can  also  influence  cold  frontal  passages  from  west  to 
east.  The  surface  front  is  impeded  by  the  relatively  high  orography  in  the  Carolinas, 
but  proceeds  more  rapidly  across  Georgia.  This  causes  a  wave-like  inflection  in  the 
low  level  front  {Type  A  cyclogenesis).  In  the  middle  and  upper  troposphere,  the 
frontal  zone  continues  its  eastward  traverse  above  the  mountains,  which  results  in 
static  destabilization  that  promotes  convective  activity  and  enhances  cyclogenesis. 

Westerly  flow  over  a  mountain  results  in  an  induced  ridge  over  upwind  slopes 
and  a  corresponding  trough  on  the  downwind  side  (Holton,  1979;  Hayes,  1985).  The 
downstream  trough  can  enhance  cyclogenesis.  A  similar  effect  may  be  caused  by  cold- 
air  damming.  A  thermal  ridge  between  the  Appalachians  and  the  Atlantic  with  a 
thermal  trough  over  the  warm  water  would  effectively  move  the  area  of  enhanced 
cyclogenesis  eastward. 

Recent  studies  of  western  Atlantic  explosive  cyclogenesis  events  all  discuss  rapid 
deepening  from  a  pre-existing  cyclone.    The  Presidents'  Day  storm  of  1979  started  as  a 
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wave  on  a  coastal  front  over  the  Carolinas  (Bosart,  1981).  Of  the  15  cases  of  rapid 
deepening  documented  by  Kocin  and  Uccellini  ( 1985a, b),  12  are  classed  as  Miller 
Type  B  (secondary  low  forming  on  the  warm  front/coastal  front)  and  the  remaining 
three  {Miller  Type  A)  moved  to  a  position  just  off  the  east  coast  before  rapid 
intensification  commenced.  Finally,  the  QE  II  storm  started  as  a  shallow  baroclinic 
system  over  the  New  Jersey  coast  (Gyakum,  1983a, b). 

C.       THE  FAVORABLE  ENVIRONMENT 

The  third  requisite  for  rapid  cyclogenesis  is  that  the  superposition  between  the 
pre-existing  surface  low  and  the  upper-level  jet  streak,  takes  place  in  an  environment 
conducive  to  cyclogenesis.  The  overwhelming  tendency  for  rapid  deepening  to  occur 
over  oceans  (Sanders  and  Gyakum,  1980;  Roebber,  1984)  provides  ample 
circumstantial  evidence  that  the  most  favorable  conditions  are  maritime.  Several 
physical  processes  have  been  suggested  to  explain  this  correlation. 

During  the  cold  season,  the  ocean  acts  as  a  huge  reservoir  of  heat  and  moisture 
for  the  atmosphere.  •  Cold  contmental  air  flowing  over  the  relatively  warm  water  gams 
heat  and  moisture  through  upward  fluxes  of  sensible  and  latent  heat.  Gall  and 
Johnson  (1971)  calculated  the  generation  of  available  potential  energy  by  sensible 
heating  during  North  Atlantic  cyclogenesis  to  be  at  least  34  percent  of  the  boundary' 
layer  frictional  dissipation.  Bunker  (1976)  and  Bunker  and  Worthington  (1976)  showed 
that  the  upward  air-sea  fluxes  during  winter  are  a  maximum  over  the  warm  water  south 
of  the  Gulf  Stream.  These  sensible  heat  fluxes  are  due  to  the  large  air-sea  temperature 
difTerence.  Additionally,  cold  season  mean  winds  are  higher  than  warm  season  values 
(U.S.  Navy,  1955),  which  enhances  the  heat  and  moisture  fluxes  from  the  ocean. 

The  direct  effect  of  the  upward  energy  flux  is  to  destabilize  the  lowest  layers  of 
the  air.  Sensible  heat  warms  the  bottom  of  the  air  mass  and  decreases  static  stability. 
Concurrently,  the  flux  of  latent  heat  increases  convective  instability.  Nuss  and  Anthes 
(1987)  found  that  fluxes  of  sensible  and  latent  heat  in  the  warm  sector  updraft 
contribute  significantly  to  cyclone  intensification. 

A  number  of  analytical  and  numerical  studies  of  general  cyclogenesis  have 
emphasized  the  important  effect  of  decreased  low-level  stability.  When  the  Richardson 
number  (R:)  of  the  mean  flow  is  small  (low  static  stability  and  high  vertical  wind 
shear),  the  growth  rate  increases  and  the  wavelength  of  maximum  instabiUty  decreases 
(Tokioka,  1973;  Gambo,  1976;  Gall,  1976a;  Staley  and  Gall,  1977;  Sandgathe.   19SI). 
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Low  R-  and  increased  convection  induced  by  air-sea  fluxes  are  also  important  in  the 
genesis  of  polar  lows  (Reed,  1979;  Mullen,  1979;  Sardie  and  Warner,  1983;  Hodur, 
1984;  Orlanski,  1984). 

In  case  studies  of  rapid  cyclogenesis,  warming  and  moistening  of  the  boundary 
layer  appear  to  be  important.  Both  Calland  (1983),  for  a  North  Pacific  case,  and  Cook 
(1983),  for  a  North  Atlantic  polar  low,  noted  a  dramatic  decrease  in  low-level  static 
stability  during  the  early  stages  of  cyclogenesis.  In  the  Presidents'  Day  case,  Bosart 
(1981)  attributed  the  location  of  the  strongest  surface  winds  (in  the  northern  semi- 
circle) to  cold  air  flowing  over  warm  water,  which  sustained  superadiabatic  lapse  rates 
near  the  surface.  Gyakum  (1983a)  cites  increasingly  deep  layers  of  potential  instability 
as  being  important  in  the  evolution  of  the  QE  II  storm. 

One  important  consequence  of  destabilization  by  air-sea  fluxes  is  a  marked 
increase  in  cumulus  convection.  The  same  processes  that  make  the  environment 
favorable  for  rapid  deepening  (sensible  heating  of  the  lowest  layer,  an  unlimited 
moisture  source,  latent  heat  flux)  also  result  in  convective  instability.  Thus,  it  seems 
that  "rapid  cyclogenesis  and  cumulus  convection  are  closely  related.  Cumulus 
convection  has  been  shown  to  enhance  baroclinic  energy  conversions  from  zonal 
available  potential  energ}'  (ZAPE)  to  eddy  available  potential  energy  (HAPE)  and  from 
EAPE  to  eddy  kinetic  energy  (Gall,  1976b;  Hayashi  and  Golder,  1981;  Robertson  and 
Smith,  1983). 

A  number  of  studies  have  connected  rapid  cyclogenesis  with  the  effects  of 
convection.  Pyke  (1965)  and  Simpson  (1969)  proposed  that  increased  cumulus 
convection  can  contribute  to  low-level  convergence,  which  increases  the  vorticity. 
Bosart  (1981)  stated  that  a  CISK-type  process  could  not  be  ruled  out  as  a  contributing 
factor  in  the  1979  Presidents'  Day  storm,  while  Bosart  and  Lin  (1984)  discovered  that 
the  vertical  advection  of  vorticity  was  a  key  factor  in  the  same  storm.  Gyakum 
( 1983a, b)  suggested  that  convection  was  important  in  the  evolution  of  the  QE  II  storm, 
and  that  a  CISK  mechanism  was  responsible  for  a  notable  increase  in  potential 
vorticity  during  the  rapid  deepening.  Anthes  et  al.  (1983)  numerically  simulated  the 
QE  II  storm,  and  found  that  cumulus  convection  and  air-sea  fluxes  of  heat  and 
moisture  were  important  in  the  development.  Rogers  and  Bosart  (1986)  also  concluded 
that  cumulus  convection  enhances  intensification. 

Convection  may  play  an  important  role  in  forming  the  "bridge"  between  the  pre- 
existing low  and  the  upper-level  jet.     Mullen  (1979)  and  Hodur  (1984)   found  that 
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convection  was  required  to  explain  the  vertical  extent  of  polar  lows.  However,  not  all 
cases  of  rapid  cyclogenesis  include  significant  cumulus  convection.  Mullen  (1982) 
documented  rapid  deepening  in  a  relatively  dry  air  mass  over  south-central  Canada. 
Although  the  low-stability,  high-moisture  environment  of  most  rapid  deepeners 
encourages  cumulus  convection,  rapid  deepening  can  occur  without  these  conditions. 

Another  mechanism  that  apparently  makes  the  maritime  environment  conducive 
to  rapid  cyclogenesis  is  enhanced  baroclinity  in  the  boundary'  layer.  The  Richardson 
number,  which  is  used  as  a  stability  parameter  in  several  studies  mentioned  above,  is 
also  a  measure  of  baroclinity.  The  vertical  wind  shear  term  in  the  denominator  is 
proportional  to  the  horizontal  temperature  gradient.  Consequently,  baroclinity  in  the 
lower  troposphere  contributes  to  lower  values  of  the  Richardson  number. 

Sanders  and  Gyakum  (1980)  presented  evidence  that  rapid  cyclogenesis  occurs 
over  the  largest  gradients  in  sea-surface  temperature  (SST).  They  also  suggested  that 
larger  deepening  rates  in  the  Atlantic  are  related  to  SST  gradients  that  are  double  those 
of  the  Pacific.  Sanders  (1986)  found  that  the  most  rapid  deepening  in  the  strongest 
Atlantic  cyclones  occurs  as  the  surface  lows  cross  the  extremely  strong  SST  gradient 
along  the  northern  edge  of  the  Gulf  Stream.  Pandolfo  (1985)  suggested  that  the 
generation  of  cyclonic  potential  vorticity  at  low  levels  might  be  significantly  enhanced 
by  SST  gradients  across  Gulf  Stream  meanders  and  eddies.  Sandgathe  (1981)  found 
that  locating  an  upper  jet  over  the  zone  of  largest  SST  gradient  was  conducive  to 
cyclogenesis.  .  ■ 

Roebber  (1984)  identified  a  local  maximum  in  occurrence  of  rapid  cyclogenesis  in 
the  Pacific  around  150°\V.  No  warm  current  is  present  in  this  region.  However, 
Roebber  maintains  that  the  maximum  is  still  due  to  increased  baroclinity  associated 
•with  outbreaks  of  Arctic  air  through  the  Bering  Straits  and  across  warmer  water. 

Since  low-level  baroclinity  is  important  to  rapid  cyclogenesis  over  the  open 
ocean,  similar  conditions  may  be  even  more  important  in  coastal  regions.  The 
proximity  of  warm  water  to  cold  land  surfaces  produces  extremely  large  horizontal 
temperature  gradients  in  a  conditionally  unstable  environment.  Baroclinically  unstable 
perturbations  developing  in  this  favorable  environment  can  be  vertically  coupled  with 
the  short  waves  or  jet  streaks  aloft. 


25 


D.      SUPERPOSITION 

Sanders  (1986)  showed  that  rapid  intensification  takes  place  as  a  500  mb  vorticity 
center  is  superposed  to  the  west  of  the  incipient  surface  cyclone.  The  relative  position 
between  the  developing  surface  low  and  the  vorticity  center  changes  during  the  life- 
cycle  of  the  cyclone.  In  composites  of  rapid  deepeners,  Sanders  showed  that  the 
average  distance  between  the  500  mb  vorticity  maximum  and  the  surface  low  pressure 
center  decreases  and  the  vorticity  center  moves  from  northwest  to  southwest  of  the 
surface  cyclone  and  then  to  the  south  as  the  low  develops. 

Farrell  (1985)  used  an  initial  value  approach  to  show  that  jet-level  and  surface 
features  exhibit  significant  growth  when  properly  superposed.  An  important 
requirement  for  rapid  growth  is  that  there  exist  a  "Fourier  overlap"  between  the  upper 
and  lower-level  disturbances.  That  is,  the  upper-level  forcing  must  have  significant 
energy  in  the  horizontal  scale  that  matches  the  surface  low.  If  the  superposition 
between  relatively  large-amplitude  features  of  similar  horizontal  scale  is  favorable, 
robust  growth  can  occur  even  when  the  waves  are  neutral  or  stable  to  exponential 
growth. 

In  summary,  rapid  cyclogenesis  is  a  phenomenon  that  depends  on  interactions 
between  processes  at  all  levels  of  the  troposphere,  and  on  several  spatial  scales.  A 
single  dominant  process  during  rapid  intensification  has  not  been  identified,  and  it  may 
vary  in  different  conditions.  Additionally,  a  unique  factor  that  differentiates  rapid 
cyclogenesis  from  "normal"  cyclogenesis  has  not  been,  established.  The  ERICA 
hypothesis  (rapid  cyclogenesis  requires  the  proper  superposition  of  PVA  aloft  and  an 
existing  surface  cyclone  in  a  favorable  environment)  represents  current  scientific 
understanding  of  the  phenomenon.  This  study  investigates  the  nature  of  the  "proper 
superposition"  factor.  The  importance  of  superposition  is  first  discussed  in  terms  of  a 
frictionless,  adiabatic  model  that  can  be  solved  analytically. 
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III.  TWO-LEVEL  MODEL  ANALYTICAL  STUDY 

A.       THE  TWO-LEVEL  MODEL 

To  examine  the  most  basic  interactions  between  low-level  and  upper-level 
perturbations,  an  analytical  study  is  conducted  using  a  simple  two-level  model.  The 
quasi-geostrophic  equations  in  pressure  coordinates  are  combined  to  form  a  finite 
difference  potential  vorticity  equation  following  Haltiner  and  Williams  (1980). 

This  simplified  linear  model  provides  a  method  to  investigate  different 
combinations  of  upper-level  and  lower-level  perturbations  in  different  mean  conditions. 
The  static  stability  and  the  vertical  shear  in  the  mean  flow  (which  is  directly 
proportional  to  the  meridional  temperature  gradient)  can  be  changed  to  study  the 
impact  on  the  development  of  waveform  disturbances.  The  wavelength,  amplitude  and 
phase  relationship  of  the  initial  upper  and  lower  perturbations  can  also  be  varied. 

By  using  the  two-level  model,  the  most  basic  baroclinic  instability  growth 
mechanism  can  be  investigated.  Increased  growth  rates  are  expected  for  decreases  in 
static  stability  and  increases  in  mean  shear,  which  can  be  related  to  increases  in  jet 
strength.  By  inputting  finite  amplitude  waves  of  different  magnitudes  and  phases  at  the 
two  levels,  the  relative  growth  rates  can  be  studied.  The  effects  of  an  intense  jet  streak 
or  shortwave  aloft  moving  over  a  weak  surface  system  can  be  estimated  and  compared 
to  the  growth  which  results  when  the  same  upper-level  feature  is  superposed  over  a 
relatively  vigorous  surface  cyclone.  Finally,  the  effects  of  different  "alongstream' 
superpositions  can  be  estimated  by  altering  the  phase  relationship  between  the  upper 
and  lower  perturbations. 

The  phase  relationship  between  upper  and  lower-level  perturbations  is  described 
intuitively  by  Hoskins  er  al.  (1985),  in  which  the  authors  argue  logically  that  the 
induced  perturbation  velocities  at  the  two  levels  will  reinforce  each  other,  which  keeps 
the  perturbations  in  step  and  enhances  gro\\th.  They  also  show  that  the  perturbations 
are  locked  on  to  each  other  with  a  phase  difference  that  depends  on  the  spatial  scale  of 
the  disturbances.  If  the  phase  is  changed  slightly,  the  amount  of  reinforcement  also 
changes,  which  alters  the  relative  phase  speeds  and  forces  the  perturbations  back  to 
their  equilibrium  phase  relationship.  The  two-level  model  is  used  to  investigate  this 
behavior    for    a    relatively     small    synoptic-scale    disturbance    in    a    low    stability 
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Figure  3.1     The  vertical  structure  of  the  two-level  model. 

The  vertical » structure  of  the  model  is.  shown  in  Fig.  Z.\..  The  atmosphere  is 
divided  into  four  layers,  with  the  wind  specified  at  the  intermediate  levels  (1  and  3)  and 
vertical  motion  specified  at  the  middle  level.  For  this  study,  p^  =  1000  mb  and 
Ap  =  400  mb. 

Assuming  waveform  solutions  for  the  streamfunction  at  the  upper  and  lower 
levels,  and  for  vertical  velocity  in  the  middle  troposphere,  the  frequency  equation  for 
the  system  can  be  determined.   The  streamiunctions  at  levels  1  and  3  are: 

T//,^3  =  Re  la^4';;3e'^^^-=*^^  +  a-^^-^e'^^'^-^"^^^   .  (3.1> 


and  the  phase  speed  is: 

where  ^  =  wave  number, 

"k  =  2fQ"  (Ap)"(y-)  is  a  stability  parameter, 

fQ  =  Coriolis  parameter, 

Ap  =  pressure  interval  (pj  -  Pq  =  2Ap), 
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^2  =  static  stability  at  level  2,  and 
U|  3  =  basic  flow  at  levels  1  and  3. 

The  shortwave  cutofi"  is  defined  by  the  relation  X  =  n^^.  For  n  greater  than  \i^ 
(shorter  wavelength)  the  phase  speed  will  be  real  and  the  wave  will  be  neutral.  The 
system  will  be  unstable  when  \i*  <  X^,  provided  (U|.U3)^  is  sufTiciently  large.  In  this 
case,  the  phase  speed  can  be  written  as  c  =  c    ±  ic:  where 

U,iU,_J|^  (3.3) 

and 

Then  (3.1)  can  be  rewritten: 

To  simplify  the  mathematics  involved  in  solving  the  two-level  model,  two 
assumptions  are  made.  First,  the  equations  are  solved  on  an  f-plane  (P  =  0).  Second, 
the  mean  wind  is  set  equal  to  zero  by  making  Ui  +  U-j  =  0.  This  condition  is 
satisfied  by  making  the  wind  at  level  1  westerly  and  the  wind  at  level  3  easterly  with  the 
same  speed.   Then  Cj.  =  0  and 

C,  =  U,(^,)      .  ,  (3.6) 

The  streamfunction  field  at  level  I  can  be  related  to  that  at  level  3: 


,±_  ,T,±  At^±  i  Vx  ^-/j,* 


^±^^±H-^'  — — 1±_^  (3.7) 

A 

Assuming  sinusoidal  initial  conditions  for  the  streamfunctions,  the  coefiicients 


a      and  a'  are: 
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where  \|/|^  and  Vj^  are  the  amplitudes  of  the  streanifunctions  at  time  zero. 

Inserting  equations  (3.6)  through  (3.9)  into  (3.5),  the  streamfunctions  at  levels  1 
and  3  can  be  determined  for  any  combination  of  wave  number,  basic  flow  and  static 
stability.  Under  the  P  =  0  assumption,  any  wave  longer  than  the  shortwave  cutoff 
will  be  unstable,  provided  Uj  is  not  equal  to  zero. 

By  combining  the  vorticity  equations  at  levels  1  and  3  and  averaging  over  x,  an 
eddy  kinetic  energy  (EKE)  equation  can  be  formed: 

d(EKE)         f 


dt     = -^  ^2'(^,-V'3)  .  O-IO) 

The  first  law  of  thermodynamics  can  be  used  to  form  an  equation  for  the  change 
in  eddy  available  potential  energy  (EAPE): 


d(EAPE)    A,.    ..  ,,,     ,  ^,av.  ^V3^  .    ^ 


dt 


=  ^(u,-U3)(i^,-i^3)(a7'-^af )  +  ^  ^2'(i^,-V3)  •         (3.11) 


The  right  side  term  of  (3.10)  is  found  with  the  opposite  sign  on  the  right  side  of 
(3.11).  This  term  is  the  conversion  of  EAPE  to  EKE.  Since  the  (Vj/j  -  vi/-^)  term  is 
proportional  to  the  temperature  perturbation,  this  term  is  the  two-level  model  variation 
of  the  (d'T' baroclinic  instability  conversion  of  EAPE  to  EKE. 

d(EKE) 


dt 


oc -cjT.  (3.12) 


Thus,  for  a  direct  circulation  (warm  air  rising  and  cold  air  sinking),  EAPE  is  converted 
to  EKE. 

A  second  term  contributes  to  the  change  in  EAPE.  (U1-U3)  is  the  shear  in  the 
mean  flow,  and  is  proportional  to  the  meridional  temperature  gradient.  The  x- 
derivative  of  the  eddy  streamlunction  can  be  expressed  as  the  meridional  wind,  v'. 
Thus,  this  term  is  proportional  to  -^T;^y(v'T')  and  represents  the  conversion  from 
zonal  available  potential  energy  (ZAPE)  to  EAPE. 

d(EAPE)       dT^ZT 


«:  --?-  V  T    +  CJ  T 


__o,__^,    ^^,    .  ^3J3^ 
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For  a  linear  model  such  as  this,  the  mean  flow  is  assumed  to  be  constant. 
Therefore,  the  (U|-U-j)  shear  remains  constant  and  provides  an  unlimited  source  of 
ZAPE.  If  the  perturbation  is  unstable,  this  ZAPE  is  continually  converted  to  EAPE, 
and  then  is  available  for  conversion  to  EKE.  Thus  unstable  perturbations  grow 
indefinitely  in  this  model. 

B.       BAROCLINICALLY  UNSTABLE  CASES 

The  validity  of  the  model  is  first  tested  against  the  theories  of  baroclinic 
instability  in  a  standard  atmosphere.  The  static  stability  was  chosen  to  correspond  to  a 
constant  lapse  rate  of  approximately  6.25  K  km'^  For  this  study,  the  growth  rate, 
exp  (Jic-t),  can  be  written 


Growth  rate  =  exp  (;zU^t^  ^^    ^   )  •  .  (3.14) 

For  a  shear  of  20  ms*^  between  800  mb  and  400  mb,  the  maximum  growth  rate  is  2.38 
day'^  at  a  wavelength  of  4000  km  (Fig.  3.2).  The  shortwave  cutolf  for  this  static 
stability  is  at  2591  km. 

The  growth  of  perturbations  in  a  standard  atmosphere  is  investigated  in 
Appendix  A.  The  model  faithfully  reproduces  expected  results.  When  a  4000  km 
perturbation  is  inserted,  it  grows  more  than  fivefold  in  the  first  48  hours  and  shifts  so 
that  a  westward  tilt  with  height  is  established.  For  a  2600  km  perturbation  (slightly 
longer  than  the  shortwave  cutoff)  a  similar  pattern  develops,  but  the  growth  rate  is 
noticeably  smaller.  The  importance  of  the  shear  in  the  mean  tlow  is  also  verified  in 
Appendix  A.  As  expected  from  (3.14),  the  growth  rate  is  directly  related  to  the  shear. 
When  the  shear  is  doubled  with  a  4000  km  perturbation,  fivefold  growth  occurs  in  less 
than  24  hours. 

As  maritime  cyclogenesis  is  generally  a  small  synoptic-scale  phenomenon  which 
occurs  in  a  low  stability  environment,  the  unstable  cases  of  interest  are  run  using  a 
static  stability  lower  than  that  of  the  standard  atmosphere  (equivalent  to  a  lapse  rate  of 
8  K  km'h.  As  the  static  stability  decreases  (k  increases),  the  growth  rate  increases  and 
the  wavelength  of  maximum  growth  shifts  to  shorter  waves,  as  can  be  seen  by 
comparing  the  low  stability  growth  rate  to  the  standard  atmosphere  growth  rate  in 
Fig.  3.2.  The  maximum  growth  rate  occurs  near  2800  km,  and  this  wavelength  is 
chosen  for  the  perturbations. 
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Figure  3.2.   Growth  rate  for  a  standard  atmosphere  (lapse  rate  6.25  K  km'^)  and  for  a 
low  stability  environment  (lapse  rate  8  K  km"^)  when  the  mean  flow  shear  is  20  m  s'  . 

The  amplitudes  of  the  initial  streamfunction  perturbations  at  levels  1  and  3  are 
varied  in  Cases  1  through  3  to  simulate  the  superpositions  of  strong  upper  features 
over  weak  or  strong  surface  cyclones.  Similarly,  weak  upper  features  in  Cases  4  and  5 
are  superposed  over  strong  and  weak  lower-levei  waves.  In  all  of  these  cases,  the  shear 
in  the  mean  flow  is  20  ms"'  between  levels  1  and  3.  with  westerly  wmds  aloft  and 
easterlies  near  the  surface. 

Cases  1  through  3  investigate  the  differences  in  development  of  low-level 
perturbations  which  have  diHering  initial  amplitudes.  The  initial  amplitude  of  the  SOO 
mb  streamfanction  in  Case  1  is  1  m~s"'.  The  corresponding  upper-level  perturbation  is 
an  order  of  magnitude  larger  (10  nrs  ).  The  initial  perturbations  are  in  phase  (there  is 
no  tilt  with  height).  The  ensuing  growth  of  both  perturbations  is  shown  in  Fig.  3.3. 
The  striking  development  in  this  experiment  is  the  extremely  rapid  growth  at  the  lower 
level.    While  the  400  mb  perturbation  triples  in  size  in  36  hours,  the  amplitude  of  the 


low-level  streamfunction  increases  from  1  m~s      to  32.5  m's    .    As  will  be  seen,  this 
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Figure  3.3.   Amplitudes  of  Case  1  400  mb  (\j/j,  solid)  and  SOO  mb  (^3,  dashed) 
perturbations  through  36  h. 


rapid  growth  to  an  amplitude  nearly  as  large  as  that  at  the  other  level  is  a 
characteristic  of  ever\'  case  in  which  one  of  the  initial  amplitudes  is  much  smaller  than 
the  other. 

For  Cases  2  and  3  the  level  1  perturbation  is  again  10  m."s"^  Case  2  is  initiated 
without  a  perturbation  at  the  lower  level.  Case  3  investigates  an  initial  lower 
perturbation  which  is  as  large  as  the  upper  feature.  As  in  Case  1,  the  initial 
perturbations  are  in  phase.  The  groAMh  of  the  800  mb  perturbations  for  Cases  2  and  3 
is  compared  to  that  for  Case  1  in  Fig.  3.4a.  The  corresponding  phase  relationships  are 
shown  in  Fig.    3.4b.    A  positive  phase  difference  indicates  a  westward  tilt  with  heighi:. 

Although  Case  2  is  initiated  without  an  SOO  mb  perturbation,  a  low-level  feature 
still  develops  and  the  pattern  is  not  appreciably  different  from  case  1.  After  only  9 
hours  the  upper  feature  lags  the  lower  by  approximately  70°  in  both  cases.  The  low- 
level  ampUtude  in  Case  2  actually  exceeds  that  from  Case  1  after  6  hours.  These 
results  indicate  that  the  presence  of  a  weak,  closed  circulation  at  the  surface  does  not 
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Figure  3.4.   (a)  Amplitudes  of  800  mb  perturbations,  and  (b)  phase  difTerences  between 
400  mb  and  800  mb  perturbations  through  36  h  for  Cases  1  (solid),  2  (dotted),  3a 
(dashed),  and  3b  with  an  initial  westward  tilt  of  90°  (hea\y  dashed). 
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appreciably  enhance  development  when  a  strong  in-phase  jet-level  feature  is 
superposed.  Indeed,  a  weak  surface  low  in  the  wrong  position  (Case  1)  will  actually  be 
detrimental  to  development  compared  to  a  case  without  an  initial  low-level 
perturbation  (Case  2). 

In  Case  3  the  initial  low-level  perturbation  is  as  large  as  the  upper  feature.  The 
amplitudes  at  both  levels  after  36  hours  are  somewhat  larger  than  for  the  cases  with 
little  or  no  initial  low-level  disturbance,  although  the  overall  intensification  at  the  lower 
level  is  actually  smaller  than  in  cases  1  and  2.  As  neither  level  dominates  in  Case  3,  the 
wave  grows  steadily  at  both  levels.  There  is  no  period  of  extremely  rapid  growth  as  the 
smaller  wave  grows  to  the  same  order  of  magnitude  as  that  of  the  other  level,  as  there 
was  in  Case  1.  After  only  9  hours  in  that  case,  the  lower  perturbation  amplitude  was 
more  than  five  times  the  initial  value.  Additionally,  the  phase  difference  between  the 
two  levels  after  9  hours  is  nearly  a  quarter-cycle.  On  the  other  hand,  when  the  lower 
perturbation  is  as  large  as  at  the  upper  level  (Case  3)  the  increase  in  the  low-level 
amplitude  early  in  the  experiment  is  not  nearly  as  dramatic.  The  growth  in  9  hours  is 
less  than  half  that  experienced  in  Case  I.  Also,  the  phase  separation  between  the 
upper  and  lower  waves  early  in  the  period  .is  much  smaller. 

It  might  be  expected  that  the  Case  3  growth,  which  starts  with  large 
perturbations  at  both  levels,  should  far  exceed  the  growth  in  cases  where  the  low-level 
perturbations  are  small.  However,  Case  3  growth  is  significantly  smaller  than  that 
expected  (3.4  day'M  from  the  growth  rate  equation  (3.14).  The  Case  3  perturbation 
actually  takes  23  hours  to  double  in  size.  The  energetics  equations  (3.12)  and  (3.13) 
help  explain  this  apparent  discrepancy. 

The  conversions  from  ZAPE  to  EAPE  and  EAPE  to  EKE  are  proportional  to  the 
correlations  between  the  eddy  quantities  v',  co'  and  T'.  The  sizes  of  these  individual 
terms  are  important,  but  the  energy  conversions  also  depend  on  the  spatial 
relationships.  Improper  phasing  can  result  m  an  indirect  circulation  (cold  air  rising  as 
warm  air  sinks)  which  means  EKE  is  converted  to  EAPE  and  the  perturbation  is 
actually  decaying. 

Synoptic  experience  indicates  that  a  perturbation  which  is  vertically  stacked  is 
not  favorable  to  further  development.  This  "barotropic"  state  is  most  often  associated 
with  the  mature  and  decaying  stages  of  synoptic-scale  cyclones.  A  more  favorable 
pattern  for  baroclinic  growth  is  a  perturbation  which  tilts  westward  with  height. 
"against"  the  mean  How  shear.    It  can  be  seen  in  Fig.  3.4b  that  this  favorable  tilt  is 


quickly  established  in  Cases  1  and  2.  However,  it  takes  15  hours  for  a  45°  phase 
difierence  to  develop  in  Case  3.  This  lack  of  favorable  superposition  explains  the 
failure  of  the  Case  3  perturbation  to  grow  as  rapidly  as  expected. 

Cases  1  and  3  are  repeated  with  an  initial  westward  tilt  of  a  quarter-wavelength. 
In  both  cases  growth  is  faster,  but  the  difference  is  especially  noticeable  in  Case  3 
(heavy  dashed  line  in  Fig.  3.4).  The  growth  rate  in  this  case  becomes  that  predicted  by 
(3.14).  The  perturbation  doubles  in  amplitude  in  only  13  hours  because  the 
superposition  is  favorable  for  baroclinic  growth  from  the  beginning  of  the  experiment. 

In  each  of  these  cases  the  phase  difference  clearly  approaches  the  same  value 
(approximately  67°).  This  westward  tilt  with  height  is  the  equilibrium  phase  difference 
described  by  Hoskins  ec  al.  (1985).  For  larger  (smaller)  phase  differences,  the  total 
perturbation  velocity  causes  the  disturbances  to  move  toward  (away  from)  each  other 
until  the  equilibrium  tilt  is  established. 

For  Case  4,  the  size  of  the  initial  upper  and  lower  perturbations  are  reversed, 
with  the  10  m  s"^  amplitude  at  800  mb  and  the  1  m^s  amplitude  at  400  mb.  The 
growth -at  the  lower  level  (Fig.  3.5)  is  ver\'  similar  to  Case  I,  even  though  the  initial 
perturbation  at  the  upper  level  is  small.  If  the  large  perturbation  is  initially  at  the 
lower  level,  the  spectacular  growth  is  at  the  upper  level,  and  the  amplitude  qilickly 
approaches  that  of  the  800  mb  wave.  As  in  Case  1,  the  tilt  of  the  wave  after  36  hours 
is  toward  the  west  \^ith  height,  with  a  phase  of  about  67°.  This  is  to  be  expected  in 
the  two.-level  model,  since  friction  is  neglected  and  the  static  stability  is  assumed- 
constant  with  height.  The  equations  for  the  streamfunctions  at  levels  1  and  3  are 
symmetric,  so  if  the  initial  conditions  are  simply  inverted,  with  no  other  changes  to  the 
mean  flow  or  stability,  the  result  must  be  growth  rates  which  are  also  reversed. 

The  development  of  the  perturbations  in  Cases  1  and  2  can  be  thought  of  as 
Petterssen  Type  B  events,  (Petterssen  and  Smebye,  1971)  in  which  an  upper-level 
positive  vorticity  maximum  is  superposed  over  a  region  which  is  favorable  for 
baroclinic  instability,  but  which  has  little  or  no  initial  vorticity  at  the  surface.  The 
Case  4  development  does  not  fit  this  pattern.  A  large  perturbation  at  the  surface  under 
an  upper  How  which  is  nearly  zonal  causes  extremely  rapid  growth  at  the  upper  level. 
This  case  is  more  like  a  Petterssen  Type  A  disturbance,  in  which  self-development 
causes  upper-level  amplification  as  the  surface  disturbance  grows  and  creates 
dilTerential  temperature  advection.  Resulting  thickness  changes  cause  an  upper-level 
ridge; trough  system  which  is  situated  such  that  positive  vorticity  advection  aloft  will 
cause  further  intensification  of  the  surface  low. 
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Figure  3.5.   Amplitudes  of  400  mb  and  SOO  mb  perturbations  through  36  h  for  Case  4. 
Vj  (solid)  and  v^  (dotted);   and  for  Case  5,   y-^  (dashed). 

If  both  the  upper  and  lower-level  initial  waves  are  small,  growth  is  ven.'  small 
compared  to  cases  in  which  either  or  both  of  the  features  is  large.  In  case  5,  the 
amplitudes  of  both  initial  waves  are  set  at  1  m'^s'^  Although  the  phase  shift  between 
the  growing  waves  is  similar  to  the  rapid  growth  cases,  the  growth  is  much  slower  than 
in  cases  1  through  4.  This  is  not  surprising,  since  the  horizontal  and  vertical  velocity 
and  temperature  perturbations  are  all  quite  small.  Thus,  the  conversions  from  ZAPE 
to  EAPE  and  EAPE  to  EKE  remain  small  and  subsequent  growth  is  limited. 

It  can  be  concluded  from  these  five  cases  that  growth  through  baroclinic 
instability  v/ill  result  when  any  perturbation  is  inserted  in  a  mean  flow  which  contains 
sufilcient  vertical  shear  (horizontal  temperature  gradient).  The  magnitude  o^  the 
gro\^ih  depends  on  the  size  of  the  initial  perturbations.  If  the  eddy  is  large  at  either 
400  mb  or  SOO  mb  (or  at  both  levels),  subsequent  growth  is  large.  These  finite-size 
perturbations  can  be  thought  of  as  linear  model  representations  of  observed  synoptic 
features  commonly  associated  with  rapid  cyclogenesis.   Actual  cases  of  cyclogenesis  are 
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rarely,  if  ever,  attributed  to  the  growth  of  small-amplitude  perturbations  in  flow  which 
is  initially  undisturbed.  Rather,  cyclogenesis  is  typically  associated  with  the  coupling  of 
large-amplitude  jet-level  perturbations  (shortwaves  and  jet  streaks)  with  low-level 
cyclonic  circulations  in  a  favorable  environment. 

The  two-level  model  indicates  that  the  vertical  location  of  the  initial  large 
perturbation  is  not  critical  to  the  resulting  growth.  Similar  intensification  will  result  if 
the  initial  large  wave  is  at  either  400  mb  or  800  mb,  which  implies  that  a  strong  surface 
wave  under  steady  flow  aloft  can  induce  rapid  cyclogenesis  as  effectively  as  a 
propagating  jet  streak  or  upper-level  shortwave.  The  addition  of  friction  and  vertically- 
varying  static  stabihty  to  the  model  may  alter  this  conclusion,  but  the  two-level  model 
alone  does  not  favor  upper-level  processes  as  more  effective  in  inducing  cyclogenesis. 

C.       PHASE  DIFFERENCES  IN  THE  INITIAL  PERTURBATIONS 

Two  other  sets  of  experiments  are  made  with  the  Case  1  and  3  initial  conditions 
to  further  study  the  effect  of  phase  differences  between  the  initial  perturbations.  Using 
the  Case  1  initial  conditions,  the  small  low-level  initial  feature  was  displaced  relative  to 
the  large  upperrlevel  feature  and  the  development  was  studied  for  a  full  range  of  phase 
differences  (Fig.  3.6).   The  initial  conditions  can  be  written 

\|/j(x,0)  =  10  cos  nx  ,  (3.15) 

and 

M/3{x,0)  =  cos  n(x  +  6)  .  (3.16) 

where  5  is  the  fraction  of  a  wavelength  phase  difference  between  the  two  perturbations. 
When  Cases  1,  3,  4  and  5  from  the  previous  section  are  integrated  with  westward 
tilts  (90°  phase  difference)  in  the  initial  conditions,  the  growth  rates  increase  over  the 
experiments  with  the  perturbations  in  phase.  A  full  range  of  initial  tilts  are  tested  to 
determine  the  most  favorable  and  least  favorable  superpositions.  In  all  cases,  the 
phase  differences  converge  to  67°  (Fig.  3.6b).  The  only  appreciable  differences  are  in 
the  amplitudes  of  the  perturbations,  where  the  fastest  growing  case  (5  =  115°)  exceeds 
the  least  favorable  case  (6  =  -65'^)  by  20  percent  after  36  hours. 
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Figure  3.6.    Ca)  .Amplitudes  of  SCO  mb  perturbations  and  (b)  phase  dilTerences  in  case  1 
through  36  h  with  initial  phase  difTerences  of  6=  -135°  (dashed),   6=  -65°  (dotted). 
6=  45°  (dot-dashed),  and  6=   115°  (solid,). 


39 


A  closer  look  at  the  growth  during  the  early  part  of  the  period  is  useful.  In  every 
case  the  phase  of  the  low-level  perturbation  quickly  shifts  so  that  a  westward  tilt  is 
established.  Beyond  18  hours  the  phase  difference  is  approximately  67°  for  ever\'  case. 
As  expected  from  baroclinic  instability  theory-,  maximum  growth  occurs  when  the  low- 
level  wave  is  initially  placed  in  a  westward  tilt  position.  When  the  800  mb  perturbation 
leads  the  400  mb  wave  by  90°  to  135°,  growth  in  the  first  several  hours  is  significantly 
greater  than  when  the  initial  perturbations  are  unfavorably  phased.  However,  after  the 
first  few  hours,  the  phase  relationship  does  not  appear  to  be  very  important  for  a  large 
perturbation  superposed  over  a  small  one.  Regardless  of  the  initial  positioning,  the  tilt 
becomes  favorable  for  baroclinic  growth  within  six  hours.  The  energy  concentrated  at 
400  mb  dominates  the  gro'^^lh,  quickly  establishing  a  lower-level  perturbation  at  a 
favorable  position.  All  perturbations  then  intensify  rapidly,  with  the  growlh  in  the 
least  favorable  case  only  four  hours  behind  that  in  the  most  favorable  case. 

It  is  hypothesized  that  the  initial  phase  relationship  between  the  upper  and  lower 
features  will  be  much  more  important  when  the  perturbations  are  approximately  the 
same  size  (Case  3).  Neither  wave  dominates,  and  it  takes  longer  for  a  shift  to  a 
favorable  westward  tilt  with  height  to  be  established.  In  this'  case  the  level  3 
perturbation  can  be  written 

\j/3(x,0)  =  cos  n(x  +  6)  .  (3.17) 

Again  the  initial  low-level  wave  is  displaced  relative  to  the  upper-level  perturbation.  If 
the  tilt  is  "downshear"  the  waves  are  actually  damped  during  the  early  period  (Fig.  3.7). 
If  the  initial  tilt  is  westward  with  height,  the  waves  at  both  levels  grow  from  the  outset. 
After  9  hours,  the  amplitude  of  the  waves  in  the  most  favorable  case  has  grown  nearly 
65  percent  over  the  initial  amplitude,  while  the  amplitude  m  the  least  favorable  case 
has  decreased  37  percent. 

In  contrast  to  the  small  amplitude  case  in  Fig.  3.6b,  the  phase  difference  changes 
slowly  during  the  early  development  when  large  amplitude  perturbations  are  present  at 
both  levels  (Fig.  3.7b).  As  the  upper  perturbation  induces  an  800  mb  disturbance 
downstream,  the  lower  perturbation  induces  a  400  mb  wave  upstream.  The  actual 
pattern  is  a  summation  of  the  two,  and  the  disturbances  gradually  move  to  the 
equilibrium  position.  In  the  least  favorable  case  (6  =  -65°)  the  phase  difference 
remains  negative  through  the  first  44  hours.    For  other  initial  phase  differences,  the  tilt 
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Figure  ..^   (a)  Amplitudes  of  800  mb  perturbations  and  (b)  phase  dilTerences  in  Case  1 
through  o6  h  with  initial  phase  diHerenccs  of  5=  -135°  (dashed).   6=  -6^°  (hea\T 
dashed).   0=  -45"  (dotted).  6=  45°  (dot-dashed),  and  6=   lI5°{soHd)  ^    ' 
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becomes  westward  with  height  somewhat  sooner.  When  the  orientation  becomes 
favorable  for  growth,  the  wave  begins  to  grow  at  both  levels.  The  amplitudes  after  36 
hours  are  markedly  different,  which  demonstrates  the  importance  of  the  initial 
superpositioning  when  both  features  are  comparatively  large.  If  the  initial  tilt  is 
favorable  {6  =  90*^  to  6  =  135°),  the  36  hour  growth  is  nearly  sevenfold,  but  the 
perturbation  continually  decays  to  less  than  15  percent  of  the  initial  value  for  an  initial 
tilt  of -65°. 

The  most  favorable  initial  phase  difference  is  115°  for  both  Case  1  and  Case  3. 
This  superposition  represents  a  balance  between  the  conversion  of  EAPE  to  EKE 
{greatest  at  180°)  and  the  conversion  of  ZAPE  to  EAPE  (maximized  at  90°).  As  the 
phase  difTerence  shifts  towards  67°,  this  case  experiences  the  greatest  early  growth. 
The  67°  phase  difTerence  is  where  the  change  in  EAPE  (conversion  from  ZAPE  minus 
loss  to  EKE)  is  the  most  positive. 

It  is  concluded  that  the  importance  of  the  initial  perturbation  superpositioning 
depends  on  the  relative  size  of  the  features.  If  the  energy  at  the  jet-level  greatly 
exceeds  that  near  the  surface,  the  lower-level  wave  will  quickly  shift/redevelop  so  that  a 
westward  tilt  with  height  is  estabhshed.  If  both  the  jet-level  and  the  surface 
circulations  are  strong,  the  superpositioning  is  critical  to  further  development.  In  this 
case,  the  upper  wave  cannot  quickly  overpower  the  low-level  feature  and  establish  a 
favorable  westward  tilt.  If  the  initial  phasing  is  unfavorable,  the  upper  and  lower 
waves  interfere  and  the  amplitudes  at  both  levels,  decrease.  Eventually  a  westward  tilt 
is  established  and  growth  begins  through  baroclinic  instability.  By  this  time,  the 
magnitude  of  the  unfavorably  phased  perturbation  is  significantly  smaller  than  that  of 
the  wave  which  had  a  favorable  tilt  at  'the  initial  time.  Therefore,  the  efliciencv  of  the 
energy  transfers  from  ZAPE  to  EAPE  and  EAPE  to  EKE  will  be  higher  for  the 
favorable  tilt  case.  After  36  hours,  the  amplitude  of  the  6  =  115°  case  perturbation  is 
nearly  50  times  the  amplitude  of  the  6'=  -65°  wave.  Comparing  Fig.  3.7a  to  Fig.  3.6a, 
when  the  Case  3  large  amplitude  perturbations  are  favorably  placed,  the  lower-level 
growth  is  significantly  greater  than  that  for  the  most  favorably  phased  Case  1  example 
with  a  large  perturbation  aloft  over  a  small  perturbation  at  the  surface. 

D.       NEUTRAL  WAVE  DEVELOPMENT 

The  static  stability  in  the  two-level  model  is  constant  with  height  and  in  time.  It 
is  recognized  that  the  value  used  in  Cases  1  through  5,  which  corresponds  to  a  lapse 
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rate  of  8  K  km'  ,  may  be  unrealistically  low  for  an  average  through  the  entire 
troposphere.  Although  low  stability  values  are  often  associated  with  oceanic 
cyclogenesis  (Sanders  and  Gyakum,  1980;  Bosart,  1981;  Sardie  and  Warner,  1983:  and 
many  others),  these  values  are  usually  confined  to  the  lower  layers  of  the  atmosphere. 
Nevertheless,  the  wavelengths  of  many  explosively  deepening  cyclones 'are  much 
shorter  than  the  4000  km  predicted  for  an  average  mid-latitude  atmosphere. 

Recent  papers  by  Farrell  (1982,  1984,  and  1985)  and  Shepherd  (1985)  discuss  the 
short-term  transient  growth  of  waves  which  are  baroclinically  stable.  The  development 
of  neutral  waves  is  also  investigated  using  the  two-level  model,  and  the  results  are 
presented  in  Appendix  A.  For  the  standard  atmosphere  static  stability,  the  shortwave 
cutoff  is  2591  km.  In  this  f-plane  study,  any  perturbation  longer  than  the  cutotT  is 
unstable.  Several  perturbation  wavelengths  shorter  than  the  critical  value  are  tested. 
The  neutral  waves  do  not  maintain  a  favorable  vertical  connection  between  the 
perturbations.  The  phase  continually  changes  throughout  the  experiment,  and  the 
amplitudes  of  the  neutral  waves  grow  (decay)  when  the  phase  dilTerence  is  positive 
(negative).  When  the  wave  tilts  "against"  the  shear,  the  conversion  from  ZAPE  to 
EAPE  to  EKE  continues  in  "a  manner  similar  to  baroclinic  instability.  However,  the 
phase  will  eventually  shift  to  an  ea.stward  tilt.  As  soon  as  the  phase  difference  changes 
sign,  the  amplitude  of  the  wave  begins  to  decay.  The  closer  the  wavelength  is  to  the 
critical  value,  the  longer  the  period  before  decay  begins  and  the  larger  the  maximum 
amplitude  of  the  perturbation. 

Different  initial  phases  are  also  tested  on  a  neutral  case  (wavelength  2000  km). 
As  in  the  unstable  cases,  introduction  of  initial  phase  differences  when  one 
perturbation  is  much  larger  than  the  other  does  not  greatly  affect  subsequent 
development.  Low- level  perturbations  in  the  most  favorable  case  (6  =  180°)  reach  a 
maximum  amplitude  only  2  m"s'^  greater  than  in  the  least  favorable  case. 

In  the  unstable  cases,  initial  phase  differences  were  critical  when  the  upper  and 
lower  features  were  the  same  size.  In  an  analogous  neutral  case  the  initial  tilt  of  the 
perturbation  is  again  crucial  to  the  development  of  the  wave  (Fig.  3.8).  The  least 
favorable  tilt  for  a  neutral  wave  is  when  the  initial  perturbation  is  vertically  stacked. 
As  can  be  seen  in  Fig.  3.8b,  the  neutral  wave  continually  shifts  toward  a  more  eastward 
tilt.  If  the  perturbations  are  initially  in  phase,  the  phase  difference  immediately 
becomes  negative  and  the  wave  converts  EKE  to  E.-XPE  and  continues  to  decrease  in 
amplitude  until  the  phase  difference  passes  180  degrees.    At  that  time,  the  next  wave  in 
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Figure  3.8.   (a)  Amplitudes  of  800  mb  perturbations  and  (b)  phase  difTerences  through 

'^''- h  for  a  neutral  case  (L  =  2000  km)  ^ "  '"        "         ^ 

90°  (dotted),  and  6=   180°  (solid). 


96  h  for  a  neutral  case  (L  =  2000  km)  with  initial  phase  differences  of  6=  0°  (dashed). 
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the  upper-level  wave  train  is  in  position  to  support  growth  and  the  perturbations  grow 
back  to  their  initial  size. 

The  most  favorable  superposition  for  the  neutral  case  occurs  when  the  initial 
perturbations  are  180  degrees  out  of  phase.  In  this  case,  the  early  development  results 
in  wave  growth  of  nearly  100  percent  before  the  phase  dilTerence  becomes  negative  and 
the  amplitude  again  decreases.  The  maximum  amplitude  of  the  oscillating  perturbation 
in  the  most  favorable  case  is  double  that  produced  when  initial  superposition  is  least 
favorable. 

E.       CONCLUSIONS 

A  two-level  linear  analytical  model  is  used  to  investigate  basic  relationships 
between  upper-level  and  lower-level  perturbations.  While  this  type  of  linear  model  is 
too  simple  for  operational  use,  it  has  great  utility  for  the  study  of  the  most  basic 
relationships  between  waves  at  an  upper  level  (approximating  the  jet  level)  and 
perturbations  present  above  the  planetary  boundary  layer. 

The  energy  transformations  through  which  the  potential  energy  of  the  basic  state 
(evident  in  the  shear  between  levels  1  and  3)  is  fed  to  the  kinetic  energy  of  the 
perturbations  can  be  equated  to  the  baroclinic  instability  transfers  first  identified  by 
Chamey  (1947)  and  Eady  (1949).  Since  the  model  is  linear,  the -eddies  do  not  interact 
with  the  mean  flow,  nor  is  the  available  potential  energy  of  the  mean  flow  decreased 
during  integration.  An  infinite  source  of  energy  is  therefore  available  for  exponential 
growth  of  perturbations  at  both  levels. 

As  long  as  perturbations  remain  relatively  small,  the  linearity  assumption  is 
reasonable  and  variations  in  growth  rates  among  the  tested  cases  are  indicative  of  what 
might  be  found  in  actual  events.  Since  the  model  is  linear,  the  size  of  the  initial 
features  can  be  scaled  as  desired  to  ensure  that  they  remain  comparatively  small 
throughout  each  integration.  The  major  interest  here  is  not  the  actual  amplitude  of  the 
perturbations,  but  the  relative  size  of  the  growth  rates  at  the  two  levels  and  for 
different  cases. 

Many  of  the  conclusions  reached  using  the  two-level  model  are  as  expected  from 
theories  derived  from  more  complicated  models  and  observational  studies.  For 
example,  energy  transformations  and  perturbation  growth  rates  are  significantly  larger 
when  static  stability  is  decreased,  and  the  growth  rate  is  directly  proportional  to  the 
shear  between  levels  1  and  3.    Baroclinic  instability  causes  growth  whenever  the  wave 
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tilts  westward  with  height  (for  mean  flows  which  become  more  westerly  with  height). 
That  is,  the  growing  waves  "oppose"  the  shear.  Study  of  shorter  neutral  waves  reveals 
that  significant  growth  can  take  place  before  the  tilt  first  becomes  eastward  with 
height. 

Several  two-level  model  conclusions  are  important  to  the  study  of  explosive 
cyclogenesis,  and  will  be  further  examined  in  following  sections: 

•  When  a  relatively  large  perturbation  is  superposed  over  a  small  one  (or  even  a 
lower  level  which  is.  without  perturbation),  the  lower  level  perturbation 
experiences  a  period  of  extremely  rapid  growth, 

•  The  phase  relationship  between  the  initial  perturbations  is  critical  if  the 
perturbations  are  about  the  same  size  at  both  levels.  If  the  superposition  is 
such  that  the  wave  tilts  eastward  with  height,  the  early  development  will  be 
damped  at  both  levels  until  a  favorable  tilt  is  established.  When  the 
perturbations  are  favorably  superposed,  growth  through  baroclinic  instability 
can  begin  at  the  outset.  Superposition  of  the  upper  and  lower  features  is  also 
important  for  the  neutral  case  if  both  perturbations  are  relatively  large.  The 
phasing  is  not  very  important  if  the  energy  is  primarily  at  one  level.  This  level 
will  dominate  the  other  and  a  favorable  westward  tilt  will  quickly  be  established. 

•  For  a  given  static  stability  and  mean  flow,  maximum  growth  occurs  when  a 
relatively  large  upper-level  perturbation  is  favorably  phased  with  respect  to  a 
well-developed  surface  feature. 
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IV.  NORAPS  CHANNEL  MODEL  EXPERIMENTS 

A.       THE  MODEL 

The  conclusions  reached  in  the  two-level  model  study  of  Chapter  III  are  further 
tested  using  a  numerical  model  that  provides  a  more  realistic  representation  of  the 
atmosphere.  The  model  chosen  is  the  \av\'  Operational  Regional  Atmospheric 
Prediction  System  (NOR.'\PS),  which  is  described  in  detail  by  Hodur  (1982,  1984). 
NOR.'XPS  is  a  primitive  equation  limited  area  model  in  (T  coordinates.  Momentum  and 
mass  variables  are  placed  on  a  type  C  staggered  grid  (Arakawa  and  Lamb.  1977).  The 
model  uses  split-explicit  time  integration,  fourth-order  advection,  and  fourth-order 
diffusion.  The  grid  is  globally  relocatable  with  variable  horizontal  resolution.  The 
number  and  location  of  the  <t  surfaces  can  be  chosen  by  the  user.  One-way  interactive 
boundar\'  conditions  are  provided  from  the  corresponding  NOGAPS  (Navy- 
Operational  Global  Atmospheric  Prediction  System)  run.  The  full  model  includes  a 
Deardorfr(1972)  planetary  boundary  layer  scheme,  drv-  convective  adjustment,  large- 
scale  precipitation,  and  a  modified  Kuo-type  cumulus  parameterization  (Kuo.  1965). 

The  experiments  m  this  chapter  use  a  channel  version-  of  XOR.APS  to  make  the 
numerical  experiments  analogous  to  the  analytical  study  of  Chapter  IIL  For  these 
experiments,  the  grid  spacing  is  set  at  100  km  in  the  x  and  y  directions.  The  channel  is 
6400  km  wide  with  free-slip,  impermeable  boundaries  on  the  northern  and  southern 
sides  and  cyclic  boundary'  conditions  in  the  east-west  direction.  The  channel  length 
depends  on  the  wavelength  of  the  longest  wave  in  each  experiment.  There  are  12  sigma 
levels  in  the  vertical  between  the  surface  and  50  mb.  Frictionless  and  adiabatic  options 
are  chosen,  and  p  is  set  to  zero. 

Because  the  experiments  in  Chapter  III  were  made  with  a  linear  model. 
NOR.A.PS  is  run  in  a  quasi-linear  mode  in  this  chapter.  The  effects  on  the  mean  flow 
due  to  nonlinear  interactions  with  the  perturbations  are  removed  by  restoring  the  mean 
flow  after  each  step  of  the  integration.  That  is.  the  zonal  average  winds,  temperatures 
and  heights  at  each  level  are  replaced  with  the  initial  zonal  mean  values  of  u.v.T  and  cp 
during  each  time  step.  This  method  allows  development  and  movement  of  all 
perturbations  without  changing  (reducing)  the  initial  sources  of  zonal  available 
potential  energy  (ZAPE)  and  zonal  kinetic  energy  (ZKE).    Wave-wave  interactions  do 
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occur,  but  are  assumed  to  be  negligible,  especially  for  perturbations  longer  than  global 
wave  number  20  (Hodur,  1984). 

A  preliminary  growth-rate  experiment  is  conducted  to  determine  the  most 
unstable  wavelength  for  conditions  that  parallel  those  of  the  two-level  model  cases. 
The  mean  wind  is  a  zonal  jet  stream  that  is  centered  in  the  middle  of  the  channel  (Fig. 
4.1).  The  width  of  the  stream  is  exaggerated  in  the  central  2000  km- wide  swath  to 
reduce  any  possible  effects  of  shear  vorticity  and  to  require  a  constant  north-south 
temperature  gradient.  This  imitates  the  conditions  in  the  two-level  model,  which  did 
not  include  variations  in  the  y-direction.  The  vertical  shear  in  the  central  region  is 
constant  with  height  from  U  =  0  at  the  surface  to  U  =  45  m  s'^  at  200  mb.  The  winds 
above  200  mb  are  constant  with  height.  This  wmd  profile  makes  the  shear  between 
400  mb  (U^gg  =  27  m  s'^)  and  800  mb  (Uggg  =  7  m  s'^)  equal  to  the  20  m  s"^  used  in 
the  two-level  cases.  The  mean  temperature  lapse  rate  (8  K  km*^)  is  equivalent  to  the 
low-stability  lapse  rate  of  Chapter  III. 

For  the  growth-rate  experiment,  initial  sinusoidal  perturbations  are  inserted  in  a 
103  X  65  grid  for  model  wave  numbers  1  through  10  (wavelengths  10,000  km  down  to 
1000  km).  The  perturbations  consist  of  two  half-sine  waves  in  (T-coordinates  so  there  is 
a  double  maximum  in  amplitude  in  the  vertical  (Fig.  4.2).  In'  this  experiment,  the 
initial  amplitudes  of  the  upper  and  lower  perturbations  are  the  same.  The  upper 
ma.ximum  is  centered  at  400  mb  (as  in  the  two-level  model)  and  extends  from  187  mb 
to  613  mb.  The  lower  maximum  is  centered  at  800  mb  (again  paralleling  the  two-level 
study)  and  extends  from  588  mb  to  the  surface.  These  perturbations  have  a  sine- 
squared  profile  in  the  cross-channel  direction,  with  zero  amplitude  at  the  northern  and 
southern  walls  and  maximum  amplitude  at  mid-channel. 

After  the  perturbations  are  inserted,  the  model  is  initialized  and  integrated  for 
84  h.  Growth  rates  (Fig.  4.3)  are  determined  for  each  wave  number  by  analyzing  the 
change  in  the  v-component  with  time.  This  curve  is  reassuringly  similar  to  the  low 
stability  growth  rate  curve  for  the  two-level  model  (Fig.  3.2).  The  wavelength  of 
maximum  growth  is  near  2800  km  (as  in  the  two-level  model),  and  the  maximum 
growth  rate  is  3.28  day'^  which  compares  to  3.38  day'^  in  the  analytical  study. 
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Figure  4.1.   Cross-section  of  the  mean  zonal  wind  (U)  field  for  NOR^^PS  channel 
model  Cases  1  through  5. 
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Figure  4.2.    Vertical  and  cross-channel  shape  of  initial  upper  and  lower  perturbations. 
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Figure  4.3.   Growth  rates  observed  in  the  NOR.APS  channel  model  for  lapse  rate 
S  K  km'^  and  the  wind  profile  of  Figure  4.1. 


B.        BASIC  NORAPS  EXPERIMENTS 

I.  E.xperiments  Without  Initial  Phase  Differences 

Based  on  the  growth  rate  curve,  a  2S00  km  perturbation  is  chosen  for  the 
remainmg  NOR.A.PS  channel  model  experiments.  The  x-dimension  of  the  channel  is 
reduced  to  one  wavelength  to  save  integration  time.  Thus,  each  experiment  includes 
only  wavelengths  equal  to  or  shorter  than  the  most  unstable  perturbation  (model  wave 
number  one).  Experiments  analogous  to  all  Chapter  3  cases  are  made  with  che 
NOR.\PS  channel  model  (Table  1).  .\s  in  the  analytical  study,  an  order  of  magnitude 
diflerence  in  initial  ampUtude  is  used  to  diflerentiate  between  strong  and  weak 
perturbations.  In  the  two-level  model,  the  perturbations  were  inserted  as  small  changes 
to  the  streamfunction  field.  In  this  chapter,  the  initial  waves  are  introduced  as  v- 
component   winds.     The   case   numbers   are   the   same   as    in   Chapter   3    for   ease   in 
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comparison,  except  that  Case  4B  in  this  chapter  is  an  additional  case.  The  magnitudes 
of  the  upper  and  lower-level  v-prime  values  in  Table  1  (0,  1  and  10  cm  s'^)  are  the 
maximum  amplitudes  of  the  upper  and  lower  perturbations.  As  they  did  in  the 
analytical  study,  the  maxima  occur  in  the  center  of  the  channel  at  400  mb  and  800  mb. 
Since  the  mean  How  is  purely  zonal,  the  v-field  at  any  time  can  be  used  to  establish  the 
position  and  amplitude  of  the  developing  perturbations. 


TABLE  1 

NORAPS  Channel  Model  Cases 

Case 

Upper 

v'  (cm  s'^) 

Lower  v'  (cm  s'^) 

1 

10 

1 

2 

10 

0 

3 

10 

10 

4A 

1 

10 

4B 

0 

10 

5 

1 

1 

Cases  1  through  3  explore  the  response  of  weak,  zero  and  strong  lower 
perturbations  under  strong  upper-level  forcing.  Cases  4A  and'4B  test  weak  and  zero 
upper  waves  with  strong  low  perturbations,  while  weak  waves  at  each  level  are  included 
in  Case  5.  As  in  the  first  set  of  two-level  cases,  these  perturbations  are  inserted  at  the 
same  zonal  position.  That  is,  there  is  initially  no  vertical  phase  difference  between  the 
perturbations. 

In  all  cases,  the  perturbations  develop  similar  profiles  of  perturbation 
amplitude  versus  height  (Fig.  4.4).  The  peak  values  (10  cm  s"M  at  the  initial  time  are 
at  400  and  800  mb.  Around  600  mb.  the  initial  perturbations  overlap  and  the 
amplitude  is  a  sum  of  contributions  from  the  upper  and  lower  features.  The  only 
growth  during  the  first  12  h  of  the  integration  is  at  the  surface  as  the  perturbation  is 
re-distributed  among  the  levels.  By  36  h,  the  basic  profile  of  the  growing  wave  is 
evident.   The  maximum  amplitude  is  at  the  surface,  as  expected  for  a  frictionless  model. 
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Figure  4.4.   Case  3  perturbation  amplitude  versus  height.    Initial  conditions  (bold). 
T  =   12  h  (dotted),  t  =  24  h  ^dashed),  t  =  36  h  (dot-dashed),  i  =  4S  h  (thin  solid). 
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and  a  secondary  maximum  is  found  near  380  mb.  A  similar  profile  develops  in  ever>' 
case,  even  when  the  initial  perturbation  is  absent  at  the  lower  (Case  2)  or  upper  level 
(Case  4B).  If  the  lower  perturbation  maximum  amplitude  had  been  originally  placed  at 
the  surface  rather  than  800  mb,  the  24  h  adjustment  period  would  have  been  much 
shorter,  and  the  amplitudes  of  the  resulting  low-level  perturbations  would  have  been 
slightly  larger.  However,  the  relative  growth  of  the  waves  in  the  different  cases  would 
have  been  unchanged. 

The  phase  of  the  perturbation  at  each  level  is  determined  by  locating  the 
maximum  v-component  and  expressing  the  distance  from  a  fixed  reference  point  as  a 
phase  angle,  where  360°  equals  one  wavelength  (2800  km).  As  expected,  a  westward 
tilt  with  height  is  quickly  established  and  is  maintained  for  the  duration  of  the 
experiment  (Fig.  4.5).  The  plots  of  phase  versus  height  can  be  quite  irregular  during 
the  first  18  to  24  h  of  the  experiments.  Again,  had  the  initial  low-level  amplitude 
maximum  been  placed  at  the  surface,  these  oscillations  would  have  been  reduced.  For 
comparison  with  the  two-level  study,  the,  phase  difierences  between  400  mb  and  800  mb 
are  calculated  by  fitting  a  least-squares  line  to  the  phases  at  all  levels  below  350  mb. 

The  first  conclusion  from  the  two-level  study  (that  lower  perturbations  will 
experience  a  period  of  rapid  growth  when  a  vigoious  upper  feature  is  superposed)  is 
verified  by  these  NORAPS  integrations.  In  all  frictionless  NORAPS  cases,  the 
amplitude  of  the  lower-level  (specifically,  at  the  surface)  perturbation  exceeds  that  at 
the  upper  level  after  24  h  of  integration,  even  if  the  initial  conditions  contain  no  lower- 
level  perturbation.  The  development  through  48  h  of  the  Case  3  perturbations  at  both 
levels  is  summarized  in  Fig.  4.6.  As  in  the  two-level  study,  a  westward  tilt  with  height 
is  quickly  established  and  the  phase  difference  eventually  reaches  an  equilibrium  value 
that  is  maintained  for  the  remainder  of  the  experiment.  This  preference  for  a  certain 
phase  difference  is  described  intuitively  by  Hoskins  et  al.  (1985)  and  was  summarized 
in  Chapter  III.  Unlike  the  two-level  study,  the  400-800  mb  phase  difference  exhibits  a 
damped  oscillation  as  the  perturbations  approach  equilibrium  separation.  As  the  upper 
and  lower  perturbations  adjust  toward  positions  where  they  provide  maximum  positive 
feedback,  they  overshoot  these  positions  and  are  then  forced  back  until  eventually  a 
balance  is  reached.  In  the  two-level  model,  the  equilibrium  phase  separation  between 
400  and  800  mb  was  approximately  67°.  In  all  NORAPS  cases,  the  eventual 
equilibrium  phase  difference  is  55°.  If  the  phase  difference  between  400  mb  and  the 
surface  (where  the  low-level  perturbation  is  largest)  is  calculated,  the  most  favorable 
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separation  is  nearly  73°.    Therefore,  the  equilibrium  separation  that  is  eventually 
established  is  nearly  the  same  in  both  models. 

The  amplitudes  of  the  initial  perturbations  in  the  NORAPS  study  are 
relatively  large  compared  to  that  of  the  streamfunction  perturbations  in  the  two-level 
study.  Since  linear  theory  assumes  infinitesimal  perturbations,  the  validity  of  finite- 
amplitude  waves  in  this  "quasi-linear"  model  may  be  suspect.  If  the  quasi-linear 
approach  is  valid,  the  amplitudes  at  any  time  should  be  proportional  to  the  amplitude 
of  the  initial  perturbations.  Using  Case  3  as  a  basis  for  comparison,  additional 
integrations  are  made  in  which  the  initial  amplitudes  are  multiplied  by  10  (v'  = 
1  m  s"  )  and  50  (v'  =  5  m  s'^).  The  upper  and  lower  perturbation  amplitudes  for  the 
largest  case  are  shown  in  Fig.  4.7.  Except  for  the  scale  of  the  y-axis,  this  case  almost 
exactly  parallels  the  small  amplitude  Case  3  (Fig.  4.6a).  Therefore,  the  quasi-linear 
approach  is  satisfactory-  for  initial  perturbations  of  this  size. 

Cases  I,  2  and  3  combine  a  vigorous  upper  perturbation  with  small,  zero  and 
large  lower  features,  respectively  (Fig.  4.8).  Only  minor  differences  exist  between  these 
cases  and  the  analogous  two-level  model  cases  (Fig.  3.4).  In  the  NORAPS  simulations, 
the  Case  1  amplitude  remains  larger  than  Case  2  throughout  the  integration.  Whereas 
the  Case  1  and  2  amplitudes  at  36  hours  are  only  about  50  percent  of  that  in  Case  3, 
the  Case  1  and  2  perturbations  in  the  two-level  study  grew  to  about  90  percent  of  the 
amplitude  of  Case  3.  In  the  two-level  model,  the  Case  3  phase  difference  shifted  to  the 
equilibrium  separation  in  about  30  h.  Using  NORAPS,  the  Case  3  tilt  becomes  nearly 
the  same  as  the  Case  1  tilt  after  about  12  h. 

In  summary,  similar  relationships  exist  in  both  models.  When  both 
perturbations  are  relatively  strong  (Case  3),  initially  unfavorable  phase  diflerences  in 
both  models  are  slower  to  shift  to  the  equilibrium  separation  compared  to  Cases  I  and 
2.  In  the  two-level  study,  the  Case  3  phase  difference  was  slow  to  reach  the  optimal 
westward  tilt,  which  resulted  in  growth  that  was  slower  than  expected.  During  the  first 
16  h  of  the  NORAPS  Case  3  integration,  the  perturbation  tilt  is  eastward  or  slightly 
westward  and  the  amplitude  of  the  lower  perturbation  remains  nearly  constant  durmg 
this  period  of  adjustment.  In  either  model,  rapid  growth  occurs  when  a  favorable 
westward  tilt  is  established.  One  difference  between  the  two  experiments  is  that  the 
zero  initial  phase  difference  superposition  is  not  as  unfavorable  for  growth  in  NOR.APS 
as  was  indicated  bv  the  two-level  model. 
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Figure  4.7.   Amplitudes  of  upper  and  lower  perturbations  when  Case  3  initial 
perturbation  amplitudes  are  multiplied  by  50. 


5S 


CASES  1.  2  and  3 

Lower  Perturbation  Amplitude 


50- 

V 

40- 

/ 

E     30- 

Case  3  / 

3 

/ 

,' 

1  -i 

a. 
1 

> 

y^      Case  1 

Case  2 

10  - 

"^                      -''    •••■■■' 

"        ''         1          '          1          '          1          '■         I 

a 

1         '         1 

Phase  Difference 


-W  120 n 

©  \ 

©  \ 

a>  90-J 
© 

3>  60- 

C 
© 

©  30- 


© 

CO 

(0 

x: 

QL     -30 


't^ 


Case  2 


Case  3 


16  24  32 

Time  (h) 


40 


48 


Figure  4.S.   (a)  Amplitudes  of  low-level  perturbations  and  (b)  piiase  dilTcrences  through 
4S  h  for  Cases  1  (dashed),  2  (dotted),  and  3  (solid). 


59 


The  lower  perturbation  amplitudes  and  phase  differences  for  the  remaining 
cases  are  shown  in  Fig.  4.9.  In  Case  4A,  an  initially  strong  lower  perturbation  is 
superposed  below  a  weak  upper  perturbation.  In  Case  4B,  the  upper  level  is  without  an 
initial  perturbation.  The  developments  in  Cases  4A  and  4B  are  similar  to  Cases  1  and 
2,  with  the  small  perturbation  case  (4A)  growing  only  slightly  faster  than  the  no- 
perturbation  case  (4B).  Case  5  initially  contained  small  perturbations  at  both  levels. 
Without  large  values  in  the  eddy  wind,  mass  and  temperature  fields,  the  energy 
conversions  in  Case  5  are  small  compared  to  those  of  Cases  1  through  4  and  the  net 
growth  of  the  perturbation  is  very  small. 

2.  Effect  of  Initial  Zonal  Phase  Differences 

An  important  conclusion  from  the  two-level  study  is  that  the  initial  zonal 
positions  of  upper  and  lower-level  perturbations  are  critical  when  both  perturbations 
are  relatively  strong.  For  example,  initial  perturbations  that  were  vertically  stacked 
(zero  initial  phase  difference)  proved  to  be  unfavorable  for  growth.  In  the  previous 
section,  it  was  indicated  that  the  in-phase  position  did  not  -seem  to  be  as  detrimental  to 
Case  3  growth  in  NORAPS.  Some  differences  between  the  two  models  that  might 
contribute  to  differences  in  the  NORAPS  study  inclu'^e:- 

•  numerical  versus  analytical  solutions; 

•  "quasi-linear"  versus  linear  approximations; 

•  vertical  resolution  that  is  greatly  increased  and  vertically-var\-ing  static  stability 
in  NORAPS;  and 

•  mass  and  temperature  fields  that  are  fixed  in  the  two-level  model  are  free  to 
adjust  to  each  other  in  NORAPS. 

To  test  the  effects  of  changes  in  the  initial  phase  difference  with  NORAPS,  the 
upper  perturbation  (10  cm  s'^  in  Cases  1  and  3,  1  cm  s"^  in  Case  4A)  is  placed  at  the 
same  location  for  each  experiment,  and  the  lower  perturbation  is  moved  east  and  west 
along  the  centerline  in  the  channel  until  the  entire  range  of  possible  zonal  phase 
relationships  is  tested.  Based  on  the  two-level  model,  the  effect  of  the  initial  phase  of 
the  superposition  is  minimal  when  one  perturbation  is  much  larger  than  the  other 
(Cases  1  and  4A).  In  the  NORAPS  cases,  the  phase  differences  also  quickly  shift  to  a 
westward  tilt  in  all  cases  and  converge  to  nearly  the  same  value  by  9  h.  The  differences 
in  growth  of  the  lower  wave  (Fig.  4.10)  for  Case  1  (upper  v'  =  10  cm  s"^  lower  v'  = 
1  cm  s    )  are  largest  between  initial  phase  differences  of  6  =  55*^  (most  favorable)  and 
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Figure  4.9.   As  in  Figure  4.S,  except  for  Cases  4A  (solid),  4B  (dashed),  and  5  (dotted). 
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6  =  -135°  (least  favorable).  The  amplitude  is  reduced  by  about  25  percent  in  the  least 
favorable  case.  While  the  equilibrium  phase  separations  that  are  eventually  reached 
are  nearly  the  same  in  the  NORAPS  and  two-level  models,  the  initial  superpositions 
that  result  in  the  greatest  (and  smallest)  growth  are  quite  difTerent.  In  the  two-level 
study  (Fig.  3.6),  the  most  (least)  favorable  initial  superposition  was  6  =  115°  (-65°). 
The  shghtly  larger  response  to  the  initial  tilt  in  the  NORAPS  model  may  be  due  to  the 
tendency  for  even  more  rapid  low-level  growth  in  NORAPS  than  in  the  two-level 
model.  As  the  lower  perturbation  grows  (even  in  the  6  =  -135°  case),  it  induces 
upper-level  development  to  the  west,  which  interferes  with  the  existing  large  upper-level 
perturbation  and  slows  the  establishment  of  the  most  favorable  westward  tilt.  In  Case 
4A  (small  upper  perturbation  superposed  on  a  large  low-level  feature)  the  effect  of 
changing  the  initial  phase  difference  is  even  less  important.  The  amplitudes  at  both 
levels  in  the  least  favorable  case  are  only  about  12  percent  less  than  the  most  favorable 
case  (not  shown). 

Where  strong  features  exist  at  both  levels  (Case  3)  in  the  two-level  model,  the 
phase  difference  is  very  slow  to  shift  to  a  favorable  westward  tilt  if  it  is  initially 
unfavorable  (Fig.  3.7).  In  contrast,  the  phase  differences  in  the  NOR^^PS  Case  3 
display  extremely  rapid  shifts  (Fig.  4.11b)  and  the  upper  and  lower  perturbations 
continually  overshoot  the  equilibrium  value.  Thus,  the  phase  differences  are  still  slow 
to  converge,  which  has  the  same  effect  on  the  amplitudes  as  in  the  two-level  study. 
That  is,  cases  that  start  near  the  equilibrium  value  (6  =  55°)  grow  quite  rapidly,  since 
the  large  perturbations  are  positioned  in  such  a  way  as  to  allow  maximum  energy 
conversion  from  zonal  available  potential  energy  (ZAPE)  to  eddy  available  potential 
energy  (EAPE)  and  from  EAPE  to  eddy  kinetic  energy  (EKE).  On  the  other  hand,  the 
phase  difference  in  the  least  favorable  case  (6  =  -135°)  oscillates  wildly  for  nearly  36  h 
before  finally  converging  on  a  55*^  westward  tilt.  Consequently,  the  perturbations  are 
in  unfavorable  superpositions  as  often  as  they  are  favorably  positioned,  and  the 
amplitude  of  the  6  =  -135°  wave  at  48  h  is  about  equal  to  the  initial  amplitude. 

Energy  budget  terms  for  selected  Case  3  runs  are  shown  in  Table  2.  The 
energy  budget  equations  are  those  of  Robertson  and  Smith  (1983)  as  modified  for  the 
closed  domain  by  Hodur  (1984).  The  dominant  budget  terms  are  the  baroclinic 
instability  processes:  CA— conversion  of  zonal  available  potential  energy  (ZAPE)  to 
eddy  available  potential  energy  (EAPE);  and  CE-conversion  of  EAPE  to  eddy  kinetic 
energy  (EKE).    The  next  most  important  terms  are  the  barotropic  term  (CK),  which 
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Figure  4.11.   As  in  Figure  4.10,  except  for  Case  3. 
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TABLE  2 

Energy  Conversions, for  Case  3 

in  units  of  10'- 

'  W  m-'. 

See  te\t 

tor  expianat 
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Tilt 

Time  (h) 

CA 

CE 

CZ 
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measures  the  conversion  of  zonal  kinetic  energy  (ZKE)  to  EKE,  and  the  conversion 
(CZ)  of  ZAPE  to  ZKE  by  vertical  redistribution  of  mass.  However,  the  magnitudes  of 
CK  and  CZ  are  only  about  5  percent  of  the  leading  terms.  The  baroclinic  terms  are 
positive  at  all  times  (instability)  and  generally  increase  exponentially  through  the 
integration.  However,  the  magnitudes  of  the  leading  terms  vary  widely  for  different 
initial  superpositions.  For  example,  CA  and  CE  in  the  least  favorable  case 
(6  =  -135^)  are  only  about  5  percent  of  the  most  favorable  case  {5  =  55°). 

Thus,  both  the  linear  model  and  NORAPS  results  demonstrate  that  the  initial 
phase  difference  is  crucial  when  upper-level  and  lower-level  perturbations  are  both 
strong.  The  most  favorable  initial  tilt  difTers  between  the  two  models  and  the  phase 
difference  changes  slowly  in  the  two-level  model  and  rapidly  in  NORAPS.  While 
shifting  from  a  115°  phase  difference  to  the  equilibrium  position,  the  two-level  model 
perturbations  grow  rapidly  as  a  result  of  consistently  large  baroclinic  energy  transfers. 
In  NORAPS,  the  same  initial  superposition  results  in  oscillations  around  the 
equilibrium  phase  difference  that  are  larger  than  the  oscillations  that  occur  when  the 
initial  ph-ise  difTerence  equals  the  equilibrium  value.  While  the  most  favorable  initial 
tilts  are  different  in  the  two  models,  both  agree  in  general  with  synoptic"  Experience  and 
theoretical  studies  (Eady,  1949;  etc.)  that  a  westward  tilt  with  height  of  about  a 
quarter-wavelength  provides  optimal  conditions  for  intensification. 

3.  Summary'  of  Zonal  Tilt  Results 

The  use  of  a  quasi-linear,  frictionless  numerical  model  with  a  constant  (low) 
static  stability  throughout  the  troposphere  confirms  the  three  most  important 
conclusions  reached  in  the  two-level  analytical  study. 

•  When  a  large  perturbation  is  present  at  the  jet  level  and,  there  is  little  or  no 
energy  at  lower  levels,  a  low-level  perturbation  will  be  induced  and  will  undergo 
a  period  of  extremely  rapid  growth.  When  this  low-level  perturbation  is  small 
relative  to  the  upper-level  perturbation,  the  initial  phase  of  the  superposition  is 
almost  irrelevant.  That  is,  the  growth  of  a  new  perturbation  near  the  most 
favorable  position  relative  to  the  large  upper-level  perturbation  dominates  over 
the  contribution  from  the  initial  low-level  perturbation. 

•  For  a  given  mean  flow  and  static  stability,  a  more  rapid  growth  occurs  when 
the  perturbation  at  the  lower  level  is  also  relatively  large  and  the  initial 
superposition  leads  to  a  westward  tilt  of  55°  with  height. 
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•  The  initial  phase  difference  of  the  superposition  is  critical  when  both  upper  and 
lower  perturbations  are  nearly  the  same  strength,  with  maximum  (minimum) 
growth  when  the  phase  difference  is  55°  (-135°). 

Thus  far,  only  along-stream  (zonal  tilt)  superposition  has  been  investigated. 
The  NOR.A.PS  model  will  be  used  to  examine  the  cross-track  superposition  in  the  next 
section. 

C.       EFFECT  OF  CROSS-CHANNEL  SUPERPOSITION 

1.  Procedure 

The  superposition  questions  investigated  in  the  previous  sections  were 
restricted  to  along-channel  (zonal)  tilts.  In  every  case,  the  upper  and  lower 
perturbations  were  centered  under  a  broad  jet  stream  in  the  middle  of  the  channel. 
Real-world  disturbances  are  not  confined  to  a  position  under  the  jet  stream.  It  is 
common  for  cyclogenesis  to  occur  south  of  the  jet  stream  and  for  the  intensifying 
system  "to  cross  under  the  jet  at  some  point  in  its  development.  The  average  position 
of  rapidly  deepening  surface  cyclones  is  in  the  cold  air  to  the  north  of  the  associated 
500  mb  forcing  (Sanders  and  Gyakum,  1980;  Rogers  and  Bosart.  1986)  but  this 
relationship  is  not  as  well  defined  as  the  westward  tilt  with  height.  Other  studies  have 
found  rapidly  deepening  cyclones  located  south  of  (or  under)  the  jet  stream  (Jannuzzi, 
1985;  Sinclair  and  Elsberr\',  1986).  The  added  dimension  of  a  cross-track  tilt  will  be 
investigated  in  this  section. 

To  study  the  effect  of  y- variations  in  the  perturbation  locations,  it  is  necessarv' 
to  use  a  more  realistic  mean  flow  than  that  from  the  previous  sections  (Fig.  4.1).  The 
jet  stream  shown  in  Fig.  4.12  has  the  same  vertical  structure  in  mid-channel  as 
previously,  in  that  the  speed  increases  linearly  with  height  from  0  m  s'^  at  the  surface 
to  45  m  s'^  at  200  mb  and  remains  constant  above  this  level.  To  satisfy  thermal  wind 
constraints,  the  horizontal  temperature  gradients  at  each  level  in  mid-channel  are  also 
equivalent  to  those  in  previous  sections.  For  example,  the  surface  temperature 
gradient  is  approximately  1.1  K,T00  km.  The  mean  winds  north  and  south  of  the  jet 
decrease  symmetrically  to  the  boundaries  so  that  horizontal  shear  effects  are 
introduced. 

The  perturbations  shown  in  Fig.  4.2  are  modified  slightly  by  making  the  y- 
dependency  a  function  of  sin"*  y  instead  of  sin"  y,  which  concentrates  perturbation 
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Figure  4.12.   As  in  Figure  4.1,  except  for  the  "basic"  case  cross-track  tilt  experiment. 
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energy  nearer  the  central  latitude  of  the  wave.  The  upper-level  perturbation  is  still 
centered  at  the  mid-point  of  the  channel,  but  the  lower-level  feature  is  placed  north  or 
south  of  niid-channel  as  well  as  east  or  west  of  the  upper  wave.  The  Case  3  amplitudes 
(10  cms'  at  both  levels)  are  used  in  this  section,  since  growth  is  most  sensitive  to 
superposition  when  both  perturbations  are  relatively  strong.  The  perturbations  are 
placed  in  a  zonal  superposition  that  is  favorable  for  growth  (55^  westward  tilt  from  800 
mb  to  400  mb)  and  the  lower  wave  is  positioned  between  800  km  north  and  800  km 
south  of  the  jet. 

Many  factors  may  influence  the  most  favorable  cross-track  position  for  the 
surface  low  pressure  center.  Not  all  of  them  can  be  included  in  this  frictionless, 
adiabatic  model  that  does  not  include  along-stream  variations  in  the  mean  flow  (jet 
streaks).  Two  factors  that  can  be  investigated  are  the  relative  contributions  of  the 
surface  baroclinity  and  shear  vorticity  effects  to  growth  of  the  low-level  perturbation. 
After  the  first  set  of  experiments  using  the  basic  jet  (Fig.  4.12),  tests  are  conducted 
using  a  mean  wind  with  enhanced  boundary  layer  shear  (strong  surface  baroclinity)  and 
using  a  mean  wind  profile  with  enhanced  mid-tropospheric  shear.  A  fourth  experiment , 
is  conducted  with  a  mean  flow  jet  that  slopes  from  south  to  north  with  height  rather 
chan  being  vertically  stacked  as  in  Fig.  4.12.  The  four  cases  are  designated  "basic," 
"baroclinic,"  "shear"  and  "tilt"  for  ease  in  reference.  In  all  cases,  the  mean  lapse  rate 
(8  K  km*^)  is  equivalent  to  that  used  in  previous  experiments. 

2.  Basic  Jet  Case 

Perturbations  centered  at  400  mb  and  800  mb  are  inserted  in  the  mean  flow 
shown  in  Fig.  4.12  with  the  upper  disturbance  in  the  center  of  the  channel  (north- 
south)  and  the  800  mb  wave  placed  55^^  to  the  east  at  a  variety  of  positions  (-800.  -600. 
-400,  -200,  0,  200,  400.  600.  800  km)  to  the  north  and  south.  Further  integrations  are 
then  made  to  isolate  the  most  favorable  meridional  tilt.  After  24  h,  the  largest 
amplitude  perturbations  are  found  at  the  surface  in  every-  integration.  Based  on  the 
relative  sizes  of  these  surface  disturbances  at  48  h  (Fig.  4.13),  the  most  favorable 
meridional  superposition  is  slightly  north  (60  km)  of  the  center  of  the  channel  for  this 
mean  flow.  The  maximum  growth  rate  is  1.62  d"^  The  diflerences  in  48  h  amplitude 
among  the  different  north-south  positions  are  relatively  small  since  the  zonal  phase 
difference  is  favorable  for  growth  in  every  case.  The  growth  in  the  least  favorable  case 
(800  km  south)  is  only  20  percent  less  than  that  for  the  60  km  north  case.    Experiments 
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Figure  4.13.   Amplitudes  of  surface  perturbations  after  4S  h  for  cases  with  low-level 
perturbations  initially  displaced  from  mid-channel.    Basic  case  (solid),  baroclinic  case 
(dot-dashed),  shear  case  (dashed),  tilt  case  (dotted). 
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with  the  initial  low-level  perturbation  more  than  1000  km  from  the  most  favorable 
meridional  position  indicate  that  the  growth  rate  decreases  even  more  rapidly  with 
increasing  separation  as  the  perturbation  is  removed  from  the  strongest  tropospheric 
baroclinic  zone. 

With  the  addition  of  y-tilt  in  the  problem,  the  most  favorable  east-west  tilt 
may  be  other  than  55°.  Several  combinations  of  zonal  and  meridional  tilts  in  the 
vicinity  of  55°,  60  km  north  are  tested,  and  it  is  found  that  this  combination  results  in 
maximum  growth  for  this  particular  mean  flow. 

The  perturbations  are  of  course  free  to  move  meridionally  and  zonally  relative 
to  the  initial  locations.  Just  as  the  zonal  phase  dilTerences  gradually  converged  on  the 
equilibrium  value  of  55°,  the  lower  perturbations  in  these  integrations  move  toward  the 
most  favorable  north-south  position  as  they  grow  (Fig.  4.14).  After  48  h.  the  average 
positions  of  the  maximum  and  minimum  surface  pressure  centers  for  all  the  basic  jet 
cases  are  within  150  km  of  the  most  favorable  position  at  60  km  north  of  mid-channel. 
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Figure  4.14.    Latitudinal  positions  of  surface  disturbances  after  48  h  for  the  basic  case. 

The  leading  energy  budget  terms  for  several  of  the  experiments  in  this  case  are 
presented  in  Table  3.  The  magnitudes  of  all  energy  transfers  for  initial  positions 
greater  than  400  km  north  or  600  km  south  of  the  center  of  the  channel  are  smaller 
than  for  the  initial  positions  nearer  mid-channel.    The  baroclinic  instability  terms  C.A 
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TABLE  3 

Energy  Conversions 

for  the  Basic  Jet  Case 

in  units  of  10'"^ 

W  m-2. 

Tilt 

Time  (h) 

CA 

CE 

CZ 

CK 

400  km  N 

12 

0.87 

0.43 

0.02 

-0.02 

24 

3.47 

2.36 

0.03 

0.08 

36 

10.23 

6.34 

-0.36 

0.45 

48 

29.90 

16.19 

-0.24 

1.48. 

200  km  N 

12 

0.89 

0.49 

0.03 

-0.01 

24 

3.48 

2.35 

0.05 

0.08 

36 

10.42 

6.42 

-0.33 

0.46 

48 

30.31 

16.56 

-0.25 

1.46 

60kmN 

12 

0.88 

0.54- 

•  0.03 

0.0 

24 

3.48 

2.31 

• 

0.06 

0.08 

.  36 

10.48      . 

6.42 

-0.32 

0.44 

48     ■ 

30.25 

16.57 

-0.23 

1.39 

200  km  S 

12 

0.84 

0.59 

0.03 

0.0 

24 

3.37 

"2.15- 

0.10 

0.07 

36 

10.19 

6.14 

-0.31 

0.39 

48 

28.96 

15.97 

-0.21 

1.23 

400  km  S 

12 

0.79 

0.62 

0.03 

0.01 

24 

3.22 

2.01 

0.13 

0.05 

36 

9.73 

5.79 

-0.28 

0.33 

48 

27.16 

15.05 

-0.19 

1.05 

600  km  S 

12 

0.74 

0.62 

0.03 

0.01 

24 

3.01 

1.82 

0.16 

0.03 

36 

9.03 

5.30 

-0.27 

0.25 

48 

24.75 

13.71 

-0.16 

0.84 

and  CE  are  by  far  the  largest  energy  conversions,  with  magnitudes  more  than  20  times 
as  large  as  the  barotropic  term  (CK)  and  the  conversion  of  ZAPE  to  ZKE  (CZ).  At 
the  earliest  time  period  {12  h),  the  largest  conversion  from  ZAPE  to  EAPE  (CA)  is 
slightly  north  of  the  most  favorable  position,  while  the  largest  CE  (conversion  from 
EAPE  to  EKE)  actually  occurs  for  perturbations  500  km  souih  of  mid-channel. 
However,  the  perturbation  at  60  km  north  grows  fastest  due  to  a  combination  of  all 
the  energy  transfers.  As  this  perturbation  continues  to  grow  larger  than  those  in  the 
other  experiments,  it  naturally  has  the  largest  magnitude  energy  conversions  for  all 
times  greater  than  12  h. 

3.  Enhanced  Surface  Baroclinitv  Case 

One  factor  that  may  be  expected  to  be  extremely  important  in  rapid 
cyclogenesis  is  surface  baroclinity.  Enhanced  boundan.^  layer  baroclinity  has  been 
suggested  as  an  ingredient  in  a  number  of  cases  of  rapidly  intensifying  cyclones 
(Sanders  and  Gyakum,  1980;  Bosart  and  Lin,  1984).  A  mean  How  that  features  a 
surface  temperature  gradient  q^  1.95  K,  100  km  (approximately  75  percent  greater  than 
that  in  the  basic  case)  is  similarly  tested  for  the  most  favorable  cross-channel  tilt.  The 
maximum  velocity  of  the  jet  is  maintained  at  45  m  s"^  but  the  vertical  wind  shear 
between  850  mb  and  200  mb  is  reduced.  Thus,  there  is  a  more  rapid  increase  m  mean 
wind  velocity  between  the  surface  and  850  mb  than  in  the  basic  jet  case  (Fig.  4.15). 
Thermal  wind  considerations  require  that  this  vertical  shear  be  balanced  by  a  large 
temperature  gradient  in  the  lowest  levels. 

The  meridional  position  of  the  fastest  growing  wave  in  this  "baroclinic"  case 
shifts  to  390  km  south  of  mid-channel  (Fig.  4.13).  As  in  the  basic  case,  all 
perturbations  grow  since  the  zonal  tilt  is  westward  with  height.  However,  the  amount 
of  growth  by  48  h  decreases  more  rapidly  with  distance  from  the  most  favorable 
position  at  390  km  south.  For  example,  the  growth  in  the  800  km  north  case  is  50 
percent  less  than  the  maximum  growth.  The  23  percent  increase  in  the  maximum 
growth  rate  (2.0  d'^)  over  the  basic  jet  case  is  attributable  to  the  increased  ZAPE  in  the 
baroclinic  jet  (Table  4).  In  ever\'  case,  the  perturbations  converge  toward  a  position 
about  500  km  south  of  mid-channel  by  48  h  (Fig.  4.16).  As  in  the  basic  case,  the 
southern  part  of  a  disturbance  placed  to  the  north  of  the  most  favorable  latitude  will 
grow  fastest,  which  shifts  the  perturbation  slightly  to  the  south.  This  process  continues 
until  the  wave  reaches  the  latitude  of  maximum  erowth. 
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Figure  4.15.   As  in  Figure  4.1,  except  for  the  "baroclinic"  case  cross-track  tilt 
experiment. 


TABLE  4 

Total  ZAPE  (10^ 

J 

48 

")  and  maximum  e 
h  for  the  cross-cr 

rowth 
lannel 

rates  (d'^)  during  the  first 
cases. 

Case 

ZAPE 

Ma.ximum  Growth  Rate 

Basic 

45.6 

1.62 

Baroclinic 

57.3 

2.00 

Shear 

.51.2 

1.64 

Tilt 

47.5 

1.51 
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Figure  4.16.   As  in  Figure  4.14,  except  for  the  baroclinic  case. 

The  sensitivity  to  small  zonal  and  meridional  displacements  about  this  55'^  390 
km  south  position  in  the  "baroclinic"  case  was  tested.  The  maximum  growth  rate 
occurs  for  a  superposition  in  which  the  low-level  perturbation  is  initially  45°  (350  km) 
east  of  and  390  km  south  of  the  upper-level  wave.    However,  this  slightly  smaller  zonal 
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phase  difTerence  results  in  only  a  1  percent  increase  in  growth  rate  over  the  55^/390  km 
south  case. 

Changes  in  the  mean  flow  may  slightly  alter  the  most  unstable  wavelength. 
Particularly,  increased  low-level  baroclinity  may  decrease  the  wavelength  of  maximum 
growth.  The  experiments  in  this  study  were  all  conducted  with  2800  km  perturbations. 
The  results  of  experiments  with  a  shorter  wavelength  should  show  only  minor 
differences,  and  the  major  conclusions  should  remain  the  same  as  in  the  2800  km  case. 

The  southward  shift  in  the  most  favorable  position  in  the  baroclinic  case  may 
be  due  to  an  enhanced  frontogenesis  mechanism  when  the  surface  temperature  gradient 
is  initially  strong.  The  surface  frontal  zone  forms  to  the  south  of  a  zonal  jet  stream 
(Reed  and  Danielson,.1959).  Stronger  frontogenesis  south  of  the  jet  results  in  greater 
baroclinic  energy  conversion  in  that  region.  Since  vertical  motions  and  temperature 
perturbations  are  especially  enhanced  in  frontal  regions,  the  conversion  from  EAPE  to 


EKE  (CE),  which  is  proportional  to  -coT',  should  be  especially  important  in  the  region 
of  the  front.  The  leading  energy  terms  for  a  number  of  integrations  in  this  case  are 
shown  in  Table  5.  The  energy  terms  for  cases  with  initial  perturbations  farther  north 
are  all  smaller  than  those  in  the  table.  Early  in  the  integrations  (12  h)  when  the 
perturbations  are  still  near  their  initial  positions,  the  maximum  CA  is  in  the  390  km 
south  case  and  the  maximum  CE  is  in  the  800  km  south  case.    The  importance  of  the 


(oT'  contribution  can  best  be  illustrated  by  the  ratio  of  CE  to  CA  at  12  h.  The  farther 
south  the  perturbation  is  placed,  the  more  important  CE  becomes.  The  ratio  increases 
from  0.32  for  the  800  km  north  case  to  0.71  for  the  800  km  south  case.  By  48  h,  all  the 
perturbations  are  south  of  the  jet  and  the  CE  to  CA  ratio  is  0.47  for  every  case,  as 
expected  for  a  linear  system,  after  the  most  unstable  mode  has  become  dominant. 

Thus,  the  most  favorable  position  (390  km  south)  represents  a  balance 
between  the  case  with  the  most  rapid  early  growth  and  the  case  with  the  largest  energy 
growth  later  in  the  integration.  At  12  h,  the  100  km  south  case  has  the  largest 
amplitude.  After  36  h,  the  600  km  south  case  has  the  largest  conversion  to  EKE,  which 
indicates  that  this  perturbation  will  eventually  dominate. 

4.  Enhanced  Upper-Tropospheric  Shear  Case 

Another  important  factor  in  the  meridional  location  of  developing  cyclones  is 
the  relative  vorticity  present  in  the  mean  flow.  For  this  experiment,  the  relative 
vorticity  of  the  mean  wind  is  due  only  to  horizontal  shear,  since  the  mean  wind  field 
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TABLE  5 

Energy  Conversions 

for  the  BarocUnic  Case  in  units  of  10'"^ 

W  m"^. 

Tilt 

Time  (h) 

CA 

CE 

CZ 

CK 

200  km  N 

12 

1.44 

0.72 

-0.03 

-0.03 

• 

24 

5.34 

2.98 

-0.01 

0.05 

36 

18.60 

8.98 

-0.84 

0.21 

48 

60.76 

28.57 

-1.08 

0.21 

Center 

12 

1.51 

0.84 

-0.05 

-0.02 

24 

5.74 

3.11 

0.01 

0.06 

36 

20.32 

9.66 

-0.85 

0.19 

48 

67.04 

31.41 

-1.07 

-0.01 

200  km  S 

12 

1.56 

0.94 

-0.05 

0.00 

, 

24 

6.10 

3.20 

0.03 

0.06' 

36 

21.83 

10.28 

-0.86 

0.14 

48 

72.82 

34.12 

-1.03 

-0.33 

390  km  S 

12 

1.58 

1.01 

-0.05 

0.01      . 

24 

6.31 

3.21 

"  0.05 

0.04 

.    36 

22.76 

10.62 

-0.84 

0.06 

48 

76.42 

35.80 

-0.94 

.  -0.69 

600  km  S 

12 

1.57 

1.06 

-0.05 

0.01 

24 

6.37 

3.18 

0.06 

0.03 

36 

23.16 

10.75 

-0.81 

-0.04 

48 

78.15 

36.62 

-0.86 

-1.15 

800  km  S 

12 

1.53 

1.08 

-0.06 

0.02 

24 

6.26 

3.06 

0.08 

0.01 

36 

22.84 

10.53 

-0.78 

-0.15 

48 

77.22 

36.20 

-0.74 

-1.55 
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has  no  curvature.  The  effect  of  the  enhanced  shear  vorticity  is  tested  by  increasing  the 
mean  wind  above  500  mb  by  a  factor  up  to  1.33  in  the  central  third  of  the  channel. 
The  new  maximum  wind  is  60  m  s"^  at  mid-channel  from  200  mb  to  the  top  of  the 
model  (Fig.  4.17).  Since  the  vertical  wind  shear  in  the  lower  troposphere  is  unchanged, 
the  surface  temperature  gradient  and  low-level  baroclinity  are  the  same  as  in  the  basic 
jet  case.  This  increased  shear  may  result  in  barotropic  instability  at  upper  levels  south 
of  the  jet.  Although  the  magnitude  of  the  barotropic  energy  conversion  term  (CK)  is 
larger  than  in  previous  cases,  it  is  still  only  10  percent  of  the  baroclinic  terms,  CA  and 
CE  (Table  6).  As  in  the  baroclinic  case,  the  changes  in  the  mean  flow  may  slightly 
alter  the  most  unstable  wavelength. 

As  in  the  previous  cases,  an  upper-level  perturbation  is  inserted  at  mid- 
channel  in  this  "shear"  case  and  lower  perturbations  are  tested  in  a  variety  of 
meridional  positions  relative  to  the  jet  stream.  All  integrations  result  in  reasonable 
growth,  regardless  of  the  initial  meridional  position  of  the  lower  perturbation  (Fig. 
4.13).  However,  the  greatest  growth  occurs  in  the  case  in  which  the  perturbation  is 
initially  200  km  north  of  the  lower  jet.  That  is,  the  increased  middle  and  upper-level 
shear  results  in  a  northward  shift  of  the  most  favorabi'e  postidn.  Further  sensitivity 
tests  indicate  maximum  growth  for  this  mean  flow  is  achieved  when  the  superposition 
includes  a  50*^  westward  tilt  and  a  200  km  north  displacement. 

The  positions  of  the  surface  perturbations  after  48  h  demonstrate  pronounced 
northward  shifts  in  these  "shear"  cases  (Fig.  4.18).  In  ever>'.case,  the  surface 
perturbations  converge  toward  a  position  approximately  350  km  north  of  the  jet 
stream.  As  in  the  previous  experiments,  perturbations  initially  placed  near  this  position 
have  the  largest  energy  conversions  at  48  h  (Table  6).  Early  in  the  integrations  (12  h) 
the  largest  conversion  from  ZAPE  to  EAPE  (CA)  coincides  with  the  most  favorable 
position  at  200  km  north  of  the  jet  for  this  experiment.  This  case  is  also  similar  to  the 
baroclinic  case  in  that  the  location  with  the  largest  amplitude  at  48  h  is  between  the 
location  with  the  greatest  early  growth  (largest  perturbations  at  12  h  are  for  the  mid- 
channel  case)  and  the  location  with  the  maximum  energy  conversion  at  later  times. 

Increasing  the  jet  strength  also  results  in  a  small  (less  than  5  percent)  increase 
in  the  maximum  lower  perturbation  gro\^i:h  rate  compared  to  the  basic  case.  This 
increase  in  the  intensity  of  the  surface  low  is  smaller  than  might  have  been  expected, 
since  the  ZAPE  for  this  case  is  12  percent  greater  than  in  the  basic  case  (Table  4). 
This  increased  ZAPE  is  introduced  in  the  upper  troposphere  and  the  relatively  small 
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Figure  4.17.   As  in  Figure  4.1,  except  for  the  "shear"  case  cross-track  tilt  experiment. 
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TABLE  6 

Ener] 

2,y  Conversions  for  the  Shear  Case 

m  units  of  10"* 

Wm-2. 

Tilt 

Time  (h) 

CA 

CE 

CZ 

CK 

600  km  N 

12 

0.90 

0.50 

0.09 

0.00 

24 

4.72 

2.92 

0.12 

0.20 

36 

14.73 

8.76 

-0.76 

1.11 

48 

45.62 

25.39 

-0.03 

4.54 

400  km  N 

12 

0.93 

0.56 

0.10 

0.01 

24 

4.73 

2.93 

0.12 

0.20 

36 

14.96 

8.89 

-0.77 

1.13 

48 

45.90 

25.68 

-0.50 

4.58 

200  km  N 

12 

0.94 

0.56 

0.11     ■ 

0.02 

* 

24 

4.70 

2.90 

0.13 

0.20 

36 

14.94 

8.84 

-0.76 

1.10 

"48 

45.25 

25.41 

-0.08 

4.49 

Center 

12 

0.91 

6.64 

0.11 

0.03 

24 

4.54 

2.78 

0.15 

0.18 

36 

14.49 

8.53 

-0.75 

1.03 

48 

43.26 

24.31 

-0.10 

4.26 

200  km  S 

12 

0.86 

0.65 

0.12 

0.05 

24 

4.30 
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36 

13.74 
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0.92 

48 

40.22 

22.67 

-0.10 

3.88 

400  km  S 

12 

0.79 

0.63 

0.12 

0.06 

24 

3.95 

2.39 

0.17 

0.13 

36 

12.60 

7.35 

-0.68 

0.80 
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36.11 

20.37 

-0.06 
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Figure  4.18.   As  in  Figure  4.14,  except  for  the  shear  case. 


vertical  shear  below  500  mb  is  not  conducive  to  transport  of  the  increased  energy 
toward  the  surface.  It  is  also  surprising  that  the  magnitudes  of  the  upper  perturbations 
(as  measured  by  the  v-component  of  the  wind)  do  not  increase  relative  to  the  basic 
case.  This  may  be  a  reflection  of  the  reduced  vertical  energy  flux,  or  it  may  be  a 
feature  of  the  quasi-linear  approach,  in  which  the  mean  zonal  wind  is  restored  to  the 
original,  value  after  each  time  step.  Although  the.  upper  level  v'  values  are  essentially 
the  same,  the  northward  flux  of  zonal  momentum  is  much  larger  for  the  enhanced  jet.- 
In  summan-',  increasing  the  jet  velocity  without  increasing  lower  tropospheric 
baroclinity  tends  to  shift  the  most  favorable  position  for  surface  low  development  to 
the  north  of  mid-channel,  but  does  not  result  in  significantly  larger  perturbation 
growth. 

5.  Tilted  Jet  Case 

The  final  experiment  in  this  section  is  conducted  with  a  mean  flow  jet  that 
slopes  northward  with  height  (Fig.  4.19).  Although  the  winds  above  200  mb  are  the 
same  as  the  basic  case,  the  meridional  location  of  the  maximum  wind  speed  tilts  from 
the  center  of  the  channel  at  200  mb  to  1000  km  south  of  mid-channel  at  the  surface. 
This  1,100  slope  is  more  realistic  than  the  symmetric  cross-sections  used  in  the 
previous    cases   (Reed,    1955;    Reed    and    Danielson.    1959).     The    maximum    surface 
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Figure  4.19.    As  in  Figure  4.1,  except  for  the  "tilt"  case  cross-track  tilt  experiment. 
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temperature  gradient  is  the  same  as  in  the  basic  case,  but  it  is  centered  380  km  south  of 
the  center  of  the  channel. 

As  in  the  other  cases,  the  westward  tilt  makes  all  positions  favorable  for 
growth.  The  meridional  position  of  the  initial  perturbation  that  results  in  the  largest 
amplitude  at  48  h  is  230  km  south  of  the  jet  (Fig.  4.13).  This  position  is  under  the 
sloping  jet  core  between  the  maximum  surface  temperature  gradient  and  the  positive 
shear  vorticity.  Based  on  the  "baroclinic"  case  results,  the  southward  shift  of  the 
baroclinic  zone  should  force  a  similar  shift  in  the  meridional  position  of  the  most 
favorable  growth  case.  However,  the  upper-tropospheric  shear  is  nearly  the  same  as 
the  basic  case,  which  favors  surface  low  intensification  just  north  of  mid-channel.  The 
relatively  large  separation  between  the  tilts  favored  by  the  shear  and  baroclinity  factors 
results  in  growlh  rates  in  this  "tilt"  case  that  are  the  smallest  (maximum  1.51  d"^)  of 
any  in  this  section,  even  though  the  ZAPE  in  the  mean  tlow  is  slightly  larger  than  in 
the  basic  case  (Table  4).  The  most  favorable  case  experiences  net  growth  that  is  41 
percent  greater  than  for  the  least  favorable  case  (800  km  north).  The  largest 
conversions  to  EKE  early  in  the  integration  is  for  the  case  with  the  perturbation  at  600 
km  south  (Table  7).  After  12  h,  the  largest  baroclinic  energy  conversions  are  for  the 
perturbation  that  starts  just  north  of  mid-channel.  In  the  tilted  case,  the  perturbations 
have  a  decreased  tendency  to  move  toward  the  most  favorable  latitude  compared  to 
previous  experiments  (Fig.  4.20).  The  conflicting  influences  of  the  low-level  baroclinity 
.to  the  south  and  positive  relative  vorticity  aloft  to  the  north  reduce  the  lateral 
movement  of  the  perturbations.  After  48  h,  the  surface  pressure  centers  are  still 
distributed  over  a  range  of  785  km,  whereas  the  separation  is  only  270  km  in  the  basic 
case  (Fig.  4.14).  The  superposition  that  results  in  the  largest  amplitude  perturbation  at 
48  h  is  a  compromise  between  the  latitude  of  the  greatest  early  growth  (600  km  south) 
and  the  latitude  of  the  wave  that  will  eventually  dominate  (50  km  north). 

6.  Summar>'  of  Cross-Channel  Tilt  Results 

Several  conclusions  can  be  drawn  from  these  experiments  in  which  different 
meridional  tilts  are  tested  within  a  simple  model  (frictionless,  adiabatic.  and  quasi- 
linear). 

•  When  upper-level  and  lower-level  perturbations  of  approximately  the  same 
amplitude  are  superposed,  the  most  important  factor  determining  subsequent 
growth   is   the   along-stream  (east-west)   tilt   between   them.     Changes   in   the 
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TABLE  7 

Energy  Conversions  for  the  Tilt  Case  in  units  of  lO''* 

Wm-2. 

Tilt 

Time  (h) 

CA 

CE 

CZ 

CK 

400  km  N 

12 

0.70 

0.33 

0.07 

-0.03 

24 

2.36 

1.77 

0.15 

0.04 

36 

6.22 

3.93 

-0.12 

.    0.31 

48 

17.00 

9.31 

0.03 

0.97 

200  km  N 

12 

0.71 

0.38 

0.08 

-0.03   . 

c* 

24 

2.39 

1.75 

0.16 

0.04 

36 

6.40 

3.97 

-0.11 

0.31 

48 

17.34 

9.57 

0.04 

0.94 

Center  . 

12 

0.71 

0.43 

0.08 

-0.02    • 

• 

24 

2.41 

1.70 

0.18 

0.04 

36 

6.44 

3.91 

'-0.11 

0.29 

48 

17.28 

9.53 

0.06 

0.87 

230  km  S 

12 

0.68 

0.47 

0.07 

-0.01 

24 

2.39 

1.59 

0.21 

0.03 

36 

6.35 

3.78 

-0.12 

0.25 

48 

16.66 

9.22 

0.09 

0.73 

400  km  S 

12 

0.65 

0.49 

0.07 

-0.01 

24 

2.33 

1.50 

0.23 

0.02 

36 

6.18 

3.60 

-0.12 

0.21 

48 

15.97 

8.82 

0.11 

0.60 

600  km  S 

12 

0.62 

0.51 

0.07 

0.00 

24 

2.27 

1.40 

0.25 

0.01 

36 

5.95 

3.36 

-0.14 

0.16 

48 

15.01 

8.29 

0.13 

0.46 
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Figure  4.20.   As  in  Figure  4.14,  except  for  the  tilt  case. 
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cross-flow  (north-south)  position  by  less  than  1000  km  only  slightly  modify  the 

growth  rates  established  by  the  zonal  superposition. 

Lower   perturbations   that    are   laterally   displaced   from   the   most-   favorable 

latitude  shift  toward  this  position  as  they  grow.    The  adjustment  in  zonal  tilt  to 

the  most  favorable  value  (all  cases  converged  to  the  equilibrium  tilt  of  55^  by 

48  h)  appears  to  be  more  rapid  than  the  convergence  in  the  meridional  tilt. 

The  most  favorable  cross-track  superposition  depends  on  the  characteristics  of 

the  mean  flow.     Increased  low-level  barocUnity  enhances   the  conversion  of 

EAPE  to  EKE  south  of  the  jet.  and  moves  the  most  favorable  position  of  the 

lower  wave  toward  the  south.    Enhanced  horizontal  shear  in  the  middle  and 

upper-troposphere  (increased  cyclonic  vorticity  north  of  the  jet)  moves  the  most 

favorable  position  of  the  lower  wave  toward  the  north. 

The  position  of  maximum  48   h  growth  in  these  experiments  is   a  balance 

between  the  position  where  early  (0-12  h)  growth  is  largest  and  the  position  that 

will  eventually  dominate  (greatest  energy  conversions  for  times  >  36  h). 

A  jet  core  that  slopes  from  south  to  north  with  height  results  in  a  region  of 

maximum  surface  temperature  gradient  that  is  laterally  separated  from  the  zone 

of  cyclonic  vorticity.    Growth  rates  in  this  case  are  smaller  than  for  a  vertically 

stacked  jet  with  the  same  velocity  and  surface  temperature  gradient. 
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D.       JET  STREAK  SUPERPOSITION 

Along-stream  maxima  in  jet  velocity  (jet  streaks)  are  important  factors  in  upper- 
level  forcing  of  cyclogenesis.  Jet  streaks  provide  an  additional  source  of  positive 
vorticity  advection,  and  superpose  a  "checkerboard"  convergence/divergence  pattern  on 
the  basic  flow.  For  example,  divergence  aloft  in  the  left  exit  and  right-entrance 
quadrants  can  enhance  existing  upward  vertical  motion  if  properly  positioned  with 
respect  to  the  surface  disturbance.  Curvature  efiects  may  also  be  added  if  the  jet  streak 
is  associated  with  an  upper-level  shortwave  trough.  In  this  section,  a  jet  streak  (rather 
than  the  400  mb  perturbation)  is  superposed  over  a  fmite-amplitude  low-level 
perturbation.  Growth  rates  at  both  levels  are  investigated  for  a  variety  of  zonal  and 
meridional  superpositions. 

.  1 .  Procedure 

A  sinusoidal  variation  in  the  x-direction  with  wavelength  2800  km  is  added  to 
the  mean  flow  in  the  basic  jet  case  (Fig.  4.12).  The  corresponding  v-component  is 
derived  by  specifying  that  the  total  flow  be  non-divergent.  A  balanced  mass  field  is 
then  derived.  Although  the  flow  is  initially  non-divergent  and  balanced,  these 
constraints  do  not  apply  during  the  integration.  The  resulting  jet-level  (200  mb)  flow 
pattern  has  a  maximum  velocity  of  52  m  s'^  and  a  mid-channel  minimum  velocity  of  36 
m  s  (Fig.  4.21).  As  the  jet  streak  provides  significant  forcing  to  lower  levels,  the  400 
mb  perturbation  that  was  previously  inserted  at  the  beginning  of  each  experiment  is 
not  necessary.  Whereas  the  width  of  the  perturbation  in  the  y-direction  was  the  same 
at  both  levels  in  previous  experiments,  the  y-scaie  of  the  jet  streak  is  different  from  that 
of  the  800  mb  perturbation.  Specifically,  v'  in  the  jet  streak  is  maximum  500  km  north 
and  south  of  the  center  of  the  channel,  and  is  zero  along  the  centerline.  This  y-scale 
mismatch  may  slightly  reduce  the  efficiency  of  the  energy  transfers  to  the  lower  levels. 
As  in  previous  sections,  the  model  is  run  in  a  frictionless.  adiabatic,  quasi-linear  mode. 
The  quasi-linear  approach  retains  the  zonal  mean  wind  throughout  the  integration,  but 
the  jet  streak,  which  can  be  considered  an  eddy  in  the  zonal  mean  flow,  is  free  to  move 
and  develop. 

The  first  jet  streak  experiment  is  conducted  without  an  initial  surface 
perturbation.  As  expected,  a  surface  system  is  quickly  induced  by  the  upper-level 
forcing.  This  development  is  used  as  a  baseline  value  for  evaluating  further  experiments 
that  include  low-level  perturbations  at  different  initial  positions  with  respect  to  the  jet 
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Figure  4.21.    Isotachs  (contours  at  5  m  s'^  intervals)  at  200  mb  for  the  initially 
non-divergent  jet  streak. 
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streak.  In  the  zonal  tilt  experiments,  the  initial  800  mb  perturbation  is  in  the  center  of 
the  channel  at  various  zonal  superpositions  relative  to  the  jet  streak.  Since  the  surface 
disturbance  induced  by  the  jet  streak  is  significantly  larger  (about  an  order  of 
magnitude)  than  the  disturbances  induced  in  section  C,  the  amplitude  of  the  initial  low- 
level  perturbation  is  increased  to  0.5  m  s"^  for  these  experiments.  After  the  most 
favorable  zonal  tilt  is  established,  cross-channel  superposition  with  the  jet  streak  is 
investigated. 

2.  Jet  Streak- Induced  Surface  Development 

With  no  initial  perturbation  at  800  mb,  the  jet  streak  quickly  forces  the 
formation  of  a  surface  system  below  the  left-exit  quadrant  of  the  streak.  The  v- 
component  amplitudes  of  the  upper  and  lower-level  disturbances  and  the  phase 
difference  between  the  jet  streak  and  the  surface  trough  through  48  h  are  shown  in  Fig. 
4.22.  As  in  previous  experiments,  the  v-component  of  the  wind  is  used  as  a  measure  of 
the  intensity  of  the  perturbations,  and  horizontal  distances  between  locations  of  the 
maximum  v-components  are  used  to  determine  the  phase  difference. 

This  case  is  analogous  to  Case  2,  in  which  a  large  upper-level  feature  induces 
perturbation  growth  in  purely  zonal-  low-level  fields.  Although  the  low-level  growth 
induced  by  the  jet  streak  (Fig.  4.22a)  is  an  order  of  magnitude  larger,  it  follows  the 
same  general  curve  as  the  Case-  2  growth  (Fig.  4.8).  Upper-level  growth  is  more 
vigorous  when  a  jet  streak  provides  the  forcing.  Whereas  the  lower  perturbation 
amplitude  exceeds  the  upper  after  24  h  in  Case  2  (not  shown),  the  upper  feature 
remains  larger  through  the  first  44  h  in  the  jet  streak  case.  The  vertical  structure  of  the 
developing  disturbance  is  also  slightly  different  between  the  two  cases.  In  the  jet  streak 
case,  the  phase  difference  changes  rapidly  ver>'  early  in  the  experiment  (Fig.  4.22b).  A 
favorable  westward  tilt  is  established  after  4  h,  and  the  phase  difference  converges 
toward  an  equilibrium  separation  nearly  the  same  as  that  in  Cases  1  through  5  (Fig. 
4.8b).  These  phase  difference  oscillations  with  a  jet  streak  present  are  less  pronounced 
than  in  previous  experiments.  However,  the  convergence  toward  the  equilibrium 
separation  is  slower.  The  separation  at  48  h  is  82.5*^,  but  is  still  decreasing  toward  an 
eventual  equilibrium  value  of  54°. 
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Figure  4.22.   (a)  Upper  and  lower  perturbation  amplitudes  and  (b)  phase  differences 
through  4S  h  for  the  jet  streak  case  without  initial  low-level  perturbation.    Positive 
(negative)  phase  differences  indicate  westward  (eastward)  tilt  with  height. 
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3.  Jet  Streak  Zonal  Tilt  Experiments 

The  efTect  of  initial  zonal  phase  diflerences  of  a  0.5  m  s'^  lower  perturbation 
relative  to  the  jet  streak  is  investigated  in  this  section.  Based  on  the  straight  jet  stream 
experiment  and  Fig.  4.22,  it  is  expected  that  a  westward  tilt  with  height  from  the  low- 
level  perturbation  to  the  streak  will  be  most  favorable  for  development.  The  phase 
difference  is  between  the  center  of  the  jet  streak  and  the  800  mb  trough.  A  full  range 
of  possible  zonal  superpositions  (0°  to  360°)is  tested,  and  the  most  favorable  position 
is  identified. 

As  in  the  straight  jet  stream  case,  the  initial  superposition  is  critical  for 
subsequent  development.  Depending  on  the  initial  tilt,  the  surface  perturbation  after 
48  h  may  be  larger  or  smaller  than  the  surface  perturbation  induced  by  the  jet  streak 
alone  (Fig.  4.23).  The  most  and  least  favorable  separations  (55^  and  -135°)  are  the 
same  as  in  the  previous  case.  Net  growth  in  the  -135°  case  is  approximately  40  percent 
less  than  in  the  55^  case.  The  phase  differences  also  follow  the  pattern  established  in 
previous  sections.  After  a  period  of  adjustment,  a  favorabk  westward  tilt  is  established 
in  every  case,  and  the  separation  between  the  800  mb  and  400  mb  features  gradually 
converges  toward  the  equilibrium  value  (52°).  As  in  the  jet  streak  case  without  a 
perturbation,  .the  phase  difference  oscillations  are  smaller  and  the  convergence  toward 
the  final  tilt  is  slower  than  in  the  straight  jet  case. 

The  feedback  between  the  upper  and  lower  levels  differs  for  the  jet  streak  case. 
In  previous  experiments,  the  superposition  that  resulted  in  the  largest  surface 
intensification  was  also  most  favorable  for  upper-level  growth.  However,  the  greatest 
perturbation  growth  above  500  mb  in  this  experiment  is  not  in  the  55°  initial  tilt  case, 
but  when  the  surface  trough  and  the  jet  streak  are  almost  vertical  (5°  westward  tilt 
with  height).  This  nearly  "barotropic"  stacking  places  the  maximum  low-level  southerly 
flow  approximately  a  quarter-wavelength  ahead  of  the  jet  streak,  with  the  left-exit 
quadrant  divergence  reinforcing  the  upward  vertical  motion  associated  with  the  low- 
level  warm  advection.  This  5°  westward  tilt  superposition  is  also  reasonably  favorable 
for  surface  intensification,  but  the  resulting  lower  perturbation  amplitude  at  48  h  is 
smaller  than  in  the  case  with  a  55'^  initial  tilt. 

4.  Jet  Streak  Cross-Channel  Tilt  Experiments 

The  sensitivity  of  the  jet  streak  case  to  cross-track  position  is  investigated  by 
placing  the  0.5  m  s'^  800  mb  perturbation  in  a  favorable  zonal  tilt  (55°)  and  at  a 
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Figure  4.23.   (a)  Amplitudes  of  low-level  perturbation  and  (b)  phase  diflerences 
through  4S  h  for  jet  streak  cases  with  initial  phase  dilTerences  of  6  =   1S0° 
(dash-dotted),  6  =  55°  (solid),  6  =  0°  (dashed),  6  =  -135°  (dotted),  and  no 
perturbation  (heav\-  solid). 
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variety  of  positions  north  or  south  of  the  jet  streak.  When  the  model  was  run  without 
a  low-level  perturbation,  the  jet  streak  induced  surface  cyclogenesis  at  a  position 
approximately  600  km  east  and  515  km  north  of  the  center  of  the  jet.  Presumably,  a 
perturbation  initially  placed  near  this  position  would  experience  the  most  favorable 
conditions  for  rapid  growth.  Indeed,  the  largest  energy  conversion  terms  at  36  and 
48  h  are  for  the  400  and  600  km  north  cases  (Table  8).  However,  the  largest  48  h 
amplitude  is  found  for  a  perturbation  initially  placed  just  south  (60  km)  of  the  center  of 
the  channel  (Fig.  4.24).  However,  the  increase  in  amplitude  from  the  least  favorable 
position  (800  km  north)  to  the  60  km  south  position  is  only  6  percent. 

All  of  the  developing  surface  centers  move  to  a  position  approximately  500  km 
north  of  the  jet  within  the  first  hour  of  the  integration.  Because  the  favorable  position 
is  shared  by  all  the  perturbations  after  1  h,  the  48  h  amplitude  depends  directly  on  the 
amplitude  of  the  perturbations  ver\'  early  in  the  integration.  After  the  first  hour,  the 
low-level  perturbation  in  the  60  km  south  case  is  largest  of  all  the  cross-channel  cases 
(approximately  10  percent  larger  than  the  500  km  north  case).  This  result  again  shows 
the  importance  of  the  surface  baroclinic  zone  to  the  very  early  growth  of  the 
perturbations.  Initial  growth  is  directly  related  to  the  low-level  temperature  gradient, 
which  is  maximum  at  mid-channel.  The  perturbation  that  starts  just  south  of  this 
position  is  shifted  northward  during  the  earliest  time  steps,  benefitting  most  from  the 
strong  baroclinity.  Maximum  growth  for  the  jet  stream  case  is  in  the  50°  westward 
tilt/60  km  south  superposition,  although  the  amplitude  is  only  marginally  larger  than 
the  55*^/60  km  south  combination. 

5.  Summar}'  of  Jet  Streak  Results 

A  non-divergent  jet  streak  is  added  to  the  zonal  basic  flow  conditions  in  the 
previous  sections  to  investigate  this  form  of  upper-level  cyclogenetic  forcing.  When  the 
model  is  run  without  a  lower  perturbation,  the  results  are  analogous  to  Case  2  (10 
cm  s'^  upper  perturbation  and  zero  lower  perturbation).  The  cyclogenesis  is  much 
more  energetic  in  the  jet  streak  case,  but  the  size  of  the  lower  perturbation  relative  to 
the  upper  feature,  and  the  phase  differences  between  them  are  quite  similar  to  Case  2. 
In  general,  the  jet  streak  experiments  support  the  results  from  previous  sections,  with 
only  minor  modifications. 

•      The  zonal  phase  difference  in  the  initial  conditions  is  extremely  important  to 

subsequent  growth.    The  most  favorable  initial  tilt  is  55*^  between  the  center  of 

the  velocitv  maximum  and  the  initial  SCO  mb  troueh. 
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TABLE  8 

Energy 

Conversions 

for  the  Jet 

Streak  Case 

in  units  of  10'^ 

Wm-2. 

Tilt 

Time  (h) 

CA 

CE 

CZ 

CK 

800  km  N 

12 

1.92 

-0.39 

0.01 

-0.06 

24 

3.27 

1.61 

-0.19 

0.13 

36 

8.20 

4:49 

-0.34 

0.45 

48 

21.39 

10.19 

-0.38 

1.74 

600  km  N 

12 

1.94 

-0.40 

0.01 

-0.06 

24 

3.25 

1.61 

-0.19 

0.13 

36 

8.21 

4.51 

-0.34 

0.44 

48 

21.53 

10.28 

-0.40 

1.74 

400  km  N 

12 

1.95 

-0.41 

0.02 

-0.06 

24 

3.21 

1.61 

.  -0.18 

0.13 

36 

8.19 

4.51 

-0.34 

0.44 

48 

21.58 

10.33 

-0.41 

1.73 

200  km  N 

12 

1.95 

-0.41 

0.02 

-0.06 

24 

3.17 

1.59 

-0.18 

0.13 

36 

8.13 

4.50 

-0.33 

0.44 

48 

21.53 

10.30 

-0.41 

1.72 

60  km  S 

12 

1,94 

-0.42 

0.02 

-0.06 

24 

3.10 

1.57 

-0.18 

0.13 

36 

8.02 

4.46 

-0.32 

0.43 

48 

21.32 

10.19 

-0.41 

1.69 

200  km  S 

12 

1.93 

-0.42 

0.02 

-0.06 

24 

3.07 

1.55 

-0.17 

0.12 

36 

7.95 

4.42 

-0.32 

0.43 

48 

21.15 

10.10 

-0.41 

1.68 
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Jet  Streak  Cross-Channel  Experiments 

Lower  Perturbation  Amplitude 
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Figure  4.24.     Lower  perturbation  amplitudes  after  I  h  and  48  h  for  the  jet  streak 
cross-track  phase  difference  cases. 
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•  Cross-track  superposition  is  not  extremely  important  in  the  jet  streak  case  since 
the  surface  cyclogenesis  occurs  at  nearly  the  same  location  relative  to  the  jet 
streak  regardless  of  the  initial  meridional  location  of  the  lower-level 
perturbation.  The  most  favorable  position  for  the  initial  low-level  perturbation 
is  at  the  location  of  the  strongest  low-level  baroclinity. 

•  Unlike  previous  experiments,  the  most  favorable  along-stream  tilt  for  upper 
perturbation  development  is  not  the  same  as  the  most  favorable  tilt  for  the 
surface  low.  When  a  jet  streak  is  present  in  the  model,  the  upper  perturbation 
grows  fastest  when  the  jet  streak  and  the  surface  trough  are  nearly  vertically 
stacked. 

E.       SUMMARY 

A  channel  version  of  NOR.'XPS  that  is  frictionless;  adiabatic,  quasi-linear  and  has 
a  constant  low-stability  lapse  rate  (8  K  km'h  is  used  to  investigate  the  superposition  of 
perturbations  centered  at  400  mb  and  800  mb.  The  simplified  model  is  chosen  so  that 
different  variables  can  be  reliably  interpreted.  This  approach  has  some  limitations  m 
direct  application  to  real-world 'cyclogenesis.  Some  additional  variables  that  are  not 
included  in  this  research  are  discussed  below. 

Sialic  siabiliiy.  The  stability  is  set  at  a  relatively  low  value  to  ensure  that  it  is 
favorable  for  cyclogenesis.  In  actual  cases,  the  low-level  static  stability  may  not  be 
favorable,  or  the  region  in  which  it  is  favorable  may  be  limited  in  horizontal  or  vertical 
extent.  Additionally,  the  stability  may  become  less  favorable  with  time.  A  more  stable 
atmosphere  would  reduce  the  growth  rate  for  all  cases  and  would  result  in  a  longer 
wavelength  of  maximum  growth,  but  would  not  change  the  conclusions  regarding 
superposition.   The  relative  growth  rates  in  the  various  cases  would  remain  the  same. 

Horizontal  scale.  A  match  in  horizontal  scales  of  the  upper  and  lower-level 
features  is  assumed.  Farrell  (1985)  used  an  initial  value  approach  to  show  that  jet-level 
and  surface  features  that  are  matched  in  size  exhibit  significant  growth  when  properly 
superposed.  However,  the  horizontal  scales  of  translating  jet-level  features  and  existing 
surface  low  pressure  centers  are  often  diflerent.  While  disturbances  with  wavelengths 
as  short  as  in  these  experiments  are  common  in  a  low  static  stability  maritime 
environment,  a  wide  range  of  upper-level  features  is  potentially  available  for  coupling. 
If  a  good  scale-match  between  the  upper  and  lower  perturbations  is  not  achieved, 
growth  will  be  reduced  in  all  cases.    Provided  a  reasonable  match  is  made  between  the 
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scales  at   the  two   levels,   the  basic  conclusions  with  respect  to   superposition  are 
unchanged. 

Timing  of  initial  superposition.  These  experiments  involve  the  development  of 
systems  in  which  the  upper  and  lower  perturbations  have  been  initially  placed  in 
certain  superpositions.  Many  cases  of  maritime  cyclogenesis  involve  the  coupling  of  a 
translating  upper  feature  with  a  quasi-stationary  surface  low  that  may  be  associated 
with  a  coastal  front  or  a  sea-surface  temperature  gradient  region.  This  scenario  would 
require  a  range  of  superpositions  from  unfavorable  to  favorable  as  the  upper  system 
moves  toward  the  surface  low.  At  the  same  time,  both  systems  may  be  changing  in 
size  and  intensity.  Presumably,  there  is  an  optimum  combination  of  timing  and 
development  that  will  minimize  negative  feedback  during  the  unfavorable  superposition 
period  and  maximize  positive  feedback,  when  the  systems  are  favorably  coupled. 

Nonlinear  diabatic  processes.  While  the  perturbations  in  these  experiments  are 
not  infinitesimal,  they  are  relatively  small  compared  to  the  "mean  flow.  Therefore,  the 
growth  rates  and  favorable  superpositions  indicated  by  these  results  are  directly 
applicable  only  to.  the  early  stages  of  cyclogenesis  cases.  As  cyclones  grow,  they 
interact  with  the-  mean  flow,  so  the  quasi-linear  approach  is  not  strictly  valid  for 
rapidly  deepening  cases.  Also,  diabatic  processes  are  most  important  when 
perturbation  amplitudes  and  associated  vertical  motions  are.  relatively  large.  Several 
experiments  were  repeated  with  physical  parameterizations  included  in  the  model 
formulation.  Although  the  growth  rates  of  the  perturbations  Were  increased  as 
physical  processes  became  important  late  in  the  integrations,  the  relative  sizes  of  the 
perturbations  in  the  various  cases  remained  the  same.  However,  it  is  expected  that  the 
convergence/divergence  pattern  associated  with  the  jet  streak,  would  be  even  more 
important  when  physical  processes  are  included.  In  that  case,  the  upper-level 
divergence  in  the  right-entrance  and  left-exit  quadrants  would  favor  development  below 
those  areas. 

With  these  limitations  in  mind,  these  experiments  suggest  a  number  of 
conclusions  that  are  applicable  to  actual  cyclogenesis  cases. 

•  For  a  given  mean  flow,  maximum  growth  rates  are  achieved  when  relatively 
large  perturbations  are  present  at  upper  and  lower  levels.  When  this  is  the  case, 
the  zonal  superposition  of  the  two  perturbations  is  crucial.  While  favorable 
coupling  results  in  the  maximum  growth  rate  possible,  an  unfavorable  phase 
difference  can  cause  decav  of  the  features  at  both  levels. 
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The  most  favorable  zonal  tilt  between  the  surface  low  and  the  upper-level  wave 
is  slightly  less  than  a  quarter  wavelength. 

The  most  favorable  position  of  the  surface  low  perpendicular  to  the  jet  stream 
depends  on  the  characteristics  of  the  mean  flow.  Several  factors  are  important 
in  the  determination  of  the  most  favorable  meridional  tilt,  including  horizontal 
shear  and  low-level  baroclinity.  Increasing  the  horizontal  shear  moves  the  most 
favorable  position  northward  while  enhancing  the  surface  temperature  gradient 
moves  it  southward.  The  location  of  the  most  rapid  growth  early  in  a 
superposition  episode  is  critically  dependent  on  the  location  of  the  maximum 
surface  temperature  gradient. 
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V.  SUPERPOSITION  EFFECTS  DURING  GALE  IOP-1 

In  preceeding  chapters,  simplified  models  were  used  to  investigate  superposition 
between  upper  and  lower-level  disturbances  in  different  mean  flows.  The  conclusions 
from  those  experiments  are  consistent  with  general  synoptic  reasoning  and  can  provide 
valuable  guidance  to  operational  forecasters  and  field  experimenters.  Cyclogenesis  in 
nature  can  be  thought  of  as  a  complex,  nonlinear  process  involving  interactions 
between  time-dependent  mean  flows  and  finite-amplitude  disturbances.  Diabatic  and 
adiabatic  processes  act  on  several  temporal  and  spatial  scales.  Since  the  final 
development  of  the  surface  cyclone  depends  on  the  combined  contributions  of  all  these 
factors,  it  is  difficult  to  isolate  the  effects  of  a  single  process. 

In  this  chapter,  a  quasi-Lagrangian  diagnostic  (QLD)  approach  is  used  to  discern 
the  various  contributions  in  one  case  of  cyclogenesis  in  the  eastern  United  States. 
While  a  single  case  study  is  not  sufficient  to  verify  whether  the  superposition 
conclusions  can  be  directly  transferred  to  operational  use,  a  thorough  diagnostic  study 
of. one  case  can  pro'vide  indications  of  the  relative  importance  of  superposition  in 
cyclogenesis.  The  case  chosen  for  this  research  occurred  during  the  first  intensive 
observational  period  (lOP)  of  the  Genesis  of  Atlantic  Lows  Experiment  (GALE)  from 
19  to  21  January  1986.  GALE  was  a  large,  cooperative  experiment  devoted  to  the 
study  of  the  mesoscale  and  air-sea  interaction  processes  involved  in  winter  cyclogenesis 
over  the  east  coast  of  the  United  States.  The  field  portion  of  the  experiment  (15 
January  to  15  March  1986)  involved  the  collection  of  extensive  observations  on  small 
time  and  space  scales  for  an  area  centered  on  North  Carolina.  However,  the  data 
available  for  this  research  are  only  those  observations  in  the  normal  synoptic  network. 
Supplementary  GALE  observations  are  not  included. 

lOP-I  was  a  case  in  which  strong  jet-level  forcing  existed,  but  the  expected  rapid 
surface  cyclogenesis  did  not  occur.  Among  the  several  possible  reasons  for  the 
relatively  weak  surface  cyclogenesis  is  the  lack  of  favorable  superposition  between  the 
upper-level  disturbance  and  the  incipient  surface  low(s).  Relatively  large  static  stability 
in  the  vicinity  of  the  surface  cyclogenesis  may  have  also  inhibited  transfer  of  the 
available  energy  from  the  upper-level  feature  to  the  developing  surface  cyclone.  The 
mass,  vorticity,  heat  and  moisture  budgets  for  IOP-1  can  be  examined  to  describe  the 
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reasons  for  the  weak  surface  low  development.  It  is  desirable  to  also  calculate  the 
budgets  based  on  model  forecasts  as  well  as  those  based  on  analyses  because  of  the 
additional  contributions  to  budget  residuals  that  arise  from  inherent  data  uncertainties 
in  the  over-ocean  observations  and  from  the  absence  of  direct  observations  of  subgrid- 
scale  processes.  When  forecasts  are  used,  the  budget  residuals  include  only  the  effects 
of  interpolation  and  truncation  errors,  in  addition  to  any  physical  processes  not 
included  in  the  budget  formulation.  Additionally,  forecast  data  are  available  as  often 
as  desired,  rather  than  at  just  12  h  synoptic  intervals.  In  this  research,  the  forecast 
budgets  are  calculated  every  6  h.  ' 

If  the  budget  based  on  the  forecast  is  to  be  used  with  confidence,  the  model 
forecast  must  have  accurately  predicted  the  development  of  the  system,  and  the 
common  terms  in  the  budget  should  be  similar  to  those  based  on  the  analysis.  In  the 
following  sections,  the  NORAPS  operational  forecast  based  on  1900  (for  brevity  all 
dates  and  times  will  be  represented  by  four  numbers  -  1900  denotes  00  GMT  19 
January  1986)  initial  conditions  will  be  verified.  The  forecast  budgets  will  be  compared 
to  analysis  budgets,  and  important  features  that  provide  clues  to  the  dynamics  of 
IOP-1  will  be  highlighted.  Finally,  the  1900  lower-tropospheric  conditions  will  be 
altered  to  simulate  conditions  that  would  be  more  conducive  to  cyclogenesis  based  on 
the  favorable  superposition  hypothesis.  The  NOR.'\PS  forecast  from  these  simulated 
initial  conditions  has  more  rapid,  low-level  cyclogenesis.  Although  this  is  only  a  single 
case,  this  simulation  does  suggest  that  absence  of  a  favorable  superposition  could  have 
contributed  to  the  absence  of  a  rapid  cyclogenesis  event  in  this  case  of  strong  upper- 
level  forcing. 

A.       NORAPS  ANALYSES  AND  FORECASTS  OF  lOP-l 

1.  Synoptic  Discussion 

This  section  contains  a  brief  description  of  the  significant  features  of  IOP-1. 
Pertle  (1987)  provides  a  more  detailed  discussion  of  the  NOR.'\PS  analyses  and 
forecasts.  Although  the  data  collection  phase  of  IOP-1  began  at  1812,  the  most 
significant  development  occurred  during  the  48  h  period  beginning  at  1900,  which  is  the 
interval  studied  in  this  chapter. 
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a.   00  GMT  19  January  1986 

The  NO  RAPS  forecasts  used  in  this  research  are  initialized  from  the  1900 
data  base.  At  this  time,  a  strong  short  wave  trough  is  located  over  the  Mississippi 
river  valley  from  the  Great  Lakes  to  the  Gulf  of  Mexico  (Fig.  5.1a).  A  polar  jet  streak 
(60  m  s"^  maximum  velocity)  is  near  the  upstream  inflection  point  and  a  subtropical  jet 
maximum  (also  60  m  s'^)  is  east  of  SC. 

The  short  wave  trough  is  clearly  evident  at  500  mb  (Fig.  5.1b),  where  the 
cyclonic  curvature  of  the  wave  combines  with  shear  to  create  an  elongated  absolute 
vorticity  maximum  (24  x  10"^  s"^)  over  lA  and  MO.  Positive  vorticity  advection  (PVA) 
is  occurring  in  and  downstream  of  the  short  wave  trough.  In  the  southeastern  U.S., 
southwesterly  flow  advects  vorticity  from  a  weak  positive  center  over  GA  northward 
toward  NC. 

The  minimum  sea-level  pressure  over  the  U.S.  (1008  mb)is  over  the  Great 
Lakes,  where  a  surface  cyclone  associated  with  the  Midwest  short  wave  is  moving 
eastward  (Fig.  5.1c).  However,  a  frontal  wave  over  KY  (1009  mb)  will  soon  become 
the  dominant  surface  system,  and  is  designated  the  "primary"  low  for  this  study.  The 
lower  tropospheric  batoclinic  region  is  in  the  cold  air  west  and  north  of  the  frontal 
system.  PVA  aloft  supports  intensification  of  the  primary  wave.  The  ridge  of  high 
pressure  over  the  Atlantic  extends  to  Cape  Hatteras,  which  results  in  southerly  flow 
and  warm  advection  .along  the  eastern  seaboard.  There  is  no  indication  at  the  surface 
of  the  500  mb  vorticity  maximum  over  GA. 

A  cross-section  from  northern  WI  across  Cape  Hatteras,  NC  to  32°N, 
71°W  (bold  line  in  Fig.  5.1b)  is  plotted  in  Fig.  5.2.  This  cross-section  is  chosen  to 
study  conditions  that  affect  the  subsequent  development  of  both  the  primary  low  and  a 
secondary  center  that  forms  by  1912.  The  cross-section  is  approximately  perpendicular 
to  the  mean  tropospheric  wind  at  1900,  so  the  perpendicular  wind  component  is  a  good 
estimate  of  the  actual  wind.  Two  jet  cores  are  present.  The  40  m  s'^  wind  maximum 
associated  with  the  Midwest  short  wave  on  the  polar  front  jet  is  at  325  mb  over  Lake 
Michigan.  The  other  maximum  (45  m  s"^)  at  250  mb  over  the  Gulf  Stream  is 
associated  with  the  subtropical  jet  streak  centered  near  32'^N,  74°W.  The  winds 
seaward  of  Cape  Hatteras  between  200  and  300  mb  do  not  appear  to  be  balanced  with 
the  analyzed  temperature  field.  Indeed,  estimates  of  the  maximum  geostrophic  wind 
speed  based  on  the  integrated  thermal  wind  equation  are  only  30  to  35  m  s"^  The 
apparent    imbalance    may    be    partly    due    to    an    ageostrophic    component    of  the 
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Figure  5.1.   XORAPS  analysis  at  00  GMT  19  Januar>'  1986.   (a)  250  mb  heights  (solid) 
m  gpm  and  isotachs  (dashed)  in  m  s"^   (b)  500  mb  heights  (solid)  in  epm  and  vorticitv 
(dashed)  x  10'^  s"'.   (c)  Sea-level  pressure  (solid)  in  mb,  1000-500  mb^iuckness 
(dashed)  in  gpm  and  frontal  positions  from  corresponding  NMC  analysis. 
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Figure  :^.2.   Perpendicular  wind  component  (solid)  in  m  s"^  and  potential  temperature 
(dashed)  m  degrees  K  m  the  00  GMT  19  Januar\-  19S6  XORAPS  analvsis  along  a 
vertical  cross-section  from  WI  to  31°X,  72°W. 
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anticyclonically  curved  subtropical  jet,  and  partly  to  analysis  errors  in  the  data-sparse 
region  over  the  ocean.  This  analysis  may  be  improved  later  when  the  dropsondes  from 
the  GALE  data  set  are  included.  The  low-level  baroclinity,  as  represented  by  the 
potential  temperature  gradient  and  the  vertical  wind  shear,  is  very  weak  ever\-where 
east  of  the  Appalachians.  In  particular,  cold  air  damming  and  a  low-level  jet  that  were 
present  in  classic  cyclogenesis  cases  such  as  the  Presidents'  Day  cyclone  are  not  present 
in  this  case. 

b.    12  GMT  19  January  1986 

After  12  h,  the  upper-level  short  wave  has  translated  eastward  and  become 
sharper  (Fig.  5.3a  and  b).  The  horizontal  scale  of  the  upper-level  system  (estimated 
from  the  distance  between  the  trough  and  the  downstream  ridge)  of  approximately 
3000  km  is  similar  to  the  systems  in  Chapters  III  and  IV.  Thus,  the  growth  rate  would 
be  expected  to  be  large  //  the  static  stability  is  sufficiently  low.  The  trough  axis  lies 
along  88^^  W  with  a  60  m  s'^  jet  maximum  upstream.  The  height  gradient  east  of  the 
trough  is  larger  than  it  was  at  1900,  with  50  m  s"^  southwesterly  winds  over  the 
Carolina  coast.  Another  jet  maximum  of  60  m  s'^  is  centered  east  of  the  Delmarva 
peninsula. 

The  500  mb  vorticity  pattern  reflects  the  sharper  trough  with  an  even- 
stronger   positive  vorticity  maximum  (28  x  10'^  s'^)   over   IL  (Fig.   5.3b).     A  weak 
vorticity  center  (14  x  10'^  s'^)  persists  just  east  of  the  Appalachians  over  the  NCSC 
border,  and  weak  PVA  continues  northeastward  across  VA. 

Both  of  these  500  mb  vorticity  maxima  have  associated  surface  features. 
The  frontal  wave  has  moved  from  KY  to  the  OH/WV  border  and  deepened  to  1004 
mb.  A  second  surface  low  pressure  center  that  is  associated  with  the  coastal  PVA  is 
located  in  southern  NC.  Since  this  coastal  system  (called  the  "secondary'"  low  in  this 
study)  is  in  a  position  that  is  theoretically  ideal  for  a  pre-existing  surface  cyclone,  it 
should  be  the  starting  point  for  rapid  intensification,  provided  other  factors  are 
favorable.  However,  this  system  forms  too  late  (after  the  short  wave  has  already 
induced  the  primar>'  low),  and  is  located  too  far  to  the  southeast  of  the  upper-level 
forcing  for  most  rapid  growth.  Based  on  a  3000  km  wavelength,  the  vertical  tilt 
between  the  coastal  low  and  the  Midwest  trough  is  approximately  105°.  Furthermore, 
the  secondare'  low  is  approximately  450  km  equatorward  of  the  axis  of  the  strongest 
PVA  associated  with  the  Midwest  trough.  As  a  result,  the  sea-level  pressure  pattern  is 
a  broad,  disorganized  area  of  low  pressure  rather  than  a  compact  center.    Southerly 
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Figure  5.3.   As  in  Figure  5.1,  except  for  12  GMT  19  January  19S6. 
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\\'inds  ahead  of  the  secondan-  low  cause  an  onshore  flow  of  warm,  moist  air  from  Cape 
Hatteras  to  Cape  Cod,  but  the  land  surface  is  too  warm  to  provide  the  contrast 
necessary-  for  a  significant  coastal  front/baroclinic  zone.  The  strongest  lower- 
troposphere  baroclinic  zones  are  west  of  the  cold  front  associated  with  the  primarv'  low 
and  much  farther  to  the  north  over  eastern  Canada. 

The  lower-tropo spheric  static  stability  is  calculated  for  the  layer  between 
1000  and  500  mb  using  the  formula  from  Holton  (1979): 

(7  =  -(a;e)^e/^p  (5.1) 

The  stability  in  the  vicinity  of  the  coastal  low  is  relatively  high  (Fig.  5.4),  since  the 
1000-500  mb  static  stability  for  a  standard  atmosphere  is  approximately  10.6  x  10"^. 
(For  comparison,  the  1000-500  mb  static  stability  of  the  low-stability  atmosphere  used 
in  Chapters  III  and  IV  is  4.7  x  10'-.)  The  minimum  stability  in  the  NOR.APS  domain 
is  less  than  10  x  10"-  cm^  g'^  mb"^  in  the  vicinity  of  the  primarv-  low  cold  front.  The 
entire  Atlantic  coast  region  north  of  Charleston  is  more  stable  than  the  standard 
atmosphere.  This  stable  environment  would  tend  to  reduce  the  feedback  between  the 
upper  and  lower  features,  and  thus  inhibit  intensification; 
c.   00  GMT  20  January  1986 

The  upper-level  trough  has  continued  to  move  eastward  with  little  change 
in  intensity  (Fig.  5.5a).  The  55  m  s"^  jet  maximum  in  the  downstream  southerly  flow. is 
centered  over  the  NJ  coast.  At  500  mb.  a  closed  low  has  developed  over  the  southern 
Appalachians  with  a  28  x  10"^  s'^  vorticity  maximum  that  is  150  km  to  the  south  of 
the  height  minimum.  PVA  continues  from  the  GA  coast  to  the  Delaware  Bay.  The 
significant  mid-troposphere  cyclogenesis  is  evident  in  the  GOES  water  vapor  imager}' 
(not  shown)  that  has  a  symmetrical  "doughnut"  oi"  relatively  dry  air  around  the  500  mb 
cut-off  low. 

The  strong  cyclogenesis  in  the  middle  and  upper  troposphere  is  not  evident 
at  the  surface  (Fig.  5.5c).  The  surface  centers  are  still  separated  by  450  km  with  the 
997  mb  primary  low  in  western  PA  and  the  996  mb  secondary'  low  over  the  southern 
Chesapeake  Bay.  It  appears  that  the  diffuse  lower  tropospheric  circulation  has  not 
been  effectively  organized  into  the  thermal  and  vertical  motion  fields  that  can  transfer 
significant  energy  into  the  low-level  system.  Additionally,  little  baroclinity  is  present  in 
the  vicinity  of  the  surface  lows.  Strong  onshore  flow  exists  from  NJ  to  Nova  Scotia, 
but  there  is  little  warm  advection  because  the  low-level  temperature  gradient  is  weak. 
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Figure  5.4.   Static  stability  in,the  1000-500  mb  layer  for  12  GMT  19  January  1986. 
Contours  are  2  x  10"^  cm   g' 
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d.  12  GMT  20  January  1986 

The  upper  cyclogenesis  has  continued  (the  500  mb  low  has  deepened  60  m 
in  the  previous  12  h)  as  it  moved  northeastward  to  the  Chesapeake  Bay  (Fig.  5.6b). 
Because  the  24  x  10'^  s'^  vorticity  center  is  nearly  concentric  with  the  height  minimum, 
only  small  mid-tropospheric  PVA  remains.  The  primar\'  surface  low  has  been  absorbed 
into  the  secondary-  system  circulation,  and  the  minimum  sea-level  pressure  is  989  mb  in 
northeastern  PA.  The  secondary  low  is  mature,  and  a  third  center  is  beginning  to 
organize  at  the  frontal  triple  point  over  RI. 

e.  00  GMT  21  January  1986 

By  the  end  of  IOP-1,  the  jet-level  short  wave,  the  500  mb  cut-off  low  and 
the  surface  cyclone  have  moved  northeast  to  New  England  and  have  begun  to  fill  as 
the  system  becomes  vertically  stacked  from  the  surface  to  the  tropopause  (Fig.  5.7). 
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Figure  5.5.   As  in  Figure  5.1,  except  for  00  GVIT  20  Januar\'  19S6. 
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Figure  5.6.   As  in  Figure  5.1.  except  for  12  GMT  20  Januar>-  1986. 
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Figure  5.7     As  in  Fig.  5.1,  except  for  00  GMT  21  January-  1986. 
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2.  NORAPS  Forecast 

The  NORAPS  forecast  from  the  1900  analysis  accurately  predicted  the 
movement,  development  and  filling  of  the  primary  low  during  the  first  24  h  of  the 
integration  (Figs.  5.8  and  5.9a).  It  also  predicted  the  development  of  the  secondary 
low  when  it  was  first  identified  as  a  low  pressure  center  on  NMC  analysis  (1918).  The 
forecast  of  sea-level  pressure  for  the  secondary  low  was  generally  excellent  (Fig.  5.9b). 
Deepening  was  slightly  slow  during  the  first  12  h,  but  the  prediction  of  the  minimum 
pressure  and  the  trend  of  intensification  followed  by  slow  filling  matches  very  well  with 
the  NORAPS  and  NMC  analyses.  The  forecast  of  the  secondary  low  movement  was 
reasonably  accurate  through  the  36  h  point  in  the  integration  (2012).  During  the  final 
12  h  of  the  prognosis,  NORAPS  forecast  the  surface  low  to  become  barotropically 
stacked  beneath  the  500  mb  cutoff  low  and.  decelerated  the  system  at  all  levels.  Thus, 
the  position  error  is  large  at  2100. 

The  36  h  forecast  fields  for  2012  are  shown  at  three  levels  in  Fig.  5.10a,  b  and 
c.  The  right  panels  (d,  e  and  f)  show  the  differences  between  the  forecast  fields  and  the 
NORAPS  analyses  valid  at  that  time.  The  errors  are  relatively  smdll  and  are  due  partly 
to  errors  in  central  height/pressure  and  partly  to  errors  in  position.  In  most  of  the  area 
of  interest,  the  forecast  cyclogenesis  was  too  weak  at  250  mb  and  500  mb  (Figs.  5.10d 
and  e).  The  surface  low  is  slightly  weaker  and  slightly  farther  south  than  the  analyzed 
center.  For  the  purposes  of  this  study,  NOR.'^PS  successfully  forecast  the  development 
of  a  complex  synoptic  pattern  with  respect  to  timing,  position  and  magnitude.  If  the 
common  terms  of  the  budgets  of  mass,  vorticity,  heat  and  moisture  based  on  these 
forecasts  also  have  a  close  resemblance  to  those  in  the  budgets  based  on  the  analysis 
fields,  then  it  can  be  assumed  that  the  forecasts  capture  the  basic  physical  processes 
that  were  important.  In  that  case,  some  inferences  can  be  made  regarding  the  relative 
importance  of  other  terms  than  can  be  only  estimated  from  the  forecast  fields. 

B.       QUASI-LAGRANGIAN  DIAGNOSTIC  BUDGETS 

The  development  of  the  primar>'  and  secondary'  IOP-1  lows  is  investigated  using 
the  technique  of  quasi- Lagrangian  diagnostics  (QLD),  which  was  developed  by 
Johnson  and  Downey  (1975)  and  was  subsequently  used  by  Wash  (1978),  Calland 
(1983)  and  Rau  (1986),  among  others.  Budgets  are  calculated  for  a  cylindrical  volume 
that  moves  with  the  cyclone,  which  eliminates  the  effects  due  to  translation.  The 
remaining  budget  terms  are  related  to  the  development  of  the  system. 
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Figure  5.S.   Analyzed  (solid)  and  forecast  (dotted)  tracks  of  the  two  lows  during  lOP- 
Labels  represent  day  and  hour  (i.e..  1900  is  00  GMT  19  Januan.'  19S6). 
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Figure  5.9.   Central  pressure  for  the  (a)  primar\-  and  (b)  secondar>-  lows.   NORAPS 
forecast  (solid).  NOR.'-XPS  analysis  (dashed)  and  NMC  analvsis  (d'otted).    X-axis  labels 
represent  day  and  hour  (i.e.,  1900  is  00  GMT  19  Januar\-  1986). 
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Figure  5.10.   (a,  b  and  c)  As  in  Figure  5.1,  except  for  NORAPS  36  h  forecast  for  1 
GMT  20  Januarv'  1986.   (d,  e  and  f)  Height/pressure  verification  (forecast  minus 
analysis)  at  (d)  250  mb,  (e)  500  mb  and  (0  surface.    Contours  in  (d)  and  (e)  are  60 

and  m  (0  are  4  mb. 
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For  this  study,  a  budget  radius  of  4°  lat.  is  chosen  to  capture  the  major  part  of 
the  circulation  involved  in  the  primary  and  secondary  lows  and  to  limit  the  influence  of 
one  center  on  the  other.  At  2000,  the  centers  are  separated  by  approximately  4°  lat. 
and  the  budget  volumes  have  a  significant  overlap.  Budgets  of  mass,  vorticity  and  heat 
are  calculated  for  both  lows  using  NORAPS  analyses  and  forecasts.  Moisture  budgets 
are  presented  for  the  forecast  fields  only  as  moisture  is  not  analyzed  in  NORAPS. 
Although  budgets  are  calculated  separately  for  the  primary  low  and  the  coastal  low, 
the  significant  development  in  this  case  occurred  in  the  secondary  system. 
Additionally,  the  life-span  of  the  primary  low  was  only  approximately  24  h.  Therefore, 
the  budgets  for  the  coastal  low  are  discussed  in  this  section  and  the  primary'  low 
budgets  are  presented  in  Appendix  C. 

1.  Mass  Budget 

The  mass  budget  is  the  simplest  QLD  budget  because  it  includes  no  sources  or 
sinks.  In  general,  the  mass  budget  for  the  coastal  low  has  a  secondary*  (vertical) 
circulation  in  which  mass  converges  into  the  budget  volume  in  the  lower  troposphere, 
moves  upward  in  the  volume,  and  then  is  exported  by  upper-level  divergence.  The 
circulation*  has  a  strong,  but  short-lived,  maximum  between  1921  and  2006  (Fig.  5.11a). 
Convergence  is  present  in  a  very  deep  layer  (surface  to  450  mb)  and  there  is  strong 
divergence  at  the  jet  level.  The  associated  upward  velocity  maximum  (not  shown)  is 
centered  near  the  level  of  non-divergence  (LND)  at  2006  and  occurs  during  the  most 
rapid  .deepening  phase  of  IOP-1.  However,  this  vertical  circulation  is  not  maintained, 
as  the  amplitude  is  diminished  by  50  percent  during  the  succeeding  12  h,  and  the  LND 
moves  toward  lower  levels.  In  the  forecast,  the  convergence  is  somewhat  stronger  near 
the  surface  and  the  LND  is  not  quite  as  high.  However,  the  rapid  increases  and 
decreases  in  vertical  circulation  are  similar  to  the  analysis. 

The  rapid  decrease  in  the  vertical  circulation  may  be  due  to  the  diffuse  nature 
of  the  surface  pressure.  The  two  lows  cause  low-level  convergence  over  a  wide  area, 
but  do  not  organize  it  into  a  well-defined  column  of  upward  motion.  The  jet-level 
divergence  maximum  at  2000  to  2006  decreases  rapidly  as  the  mass  circulation 
diminishes.  This  decrease  is  associated  with  the  movement  of  the  off-shore  jet  streak 
away  from  the  surface  low  by  2012  (Fig.  5.6a).  If  favorable  superposition  had  occurred 
between  the  available  upper-level  forcing  and  a  single  surface  system,  the  enhanced 
mass  circulation  probably  would  have  been  established  early  and  maintained  for  a 
longer  time. 
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Figure  5.11.   Horizontal  mass  flux  across  the  4*^  lat.  radius  in  the  (a)  analysis  and  (b) 
forecast  of  the  secondary*  low.   Solid  (dashed)  contours  denote  outflow  from  (inflow 
into)  the  4°  lat.  radius  budget  volume.   Contours  are  5  x  10'^"  g  s'^  mb'^ 
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2. .  Vorticity  Budget 

Development  of  the  low  can  be  described  in  terms  of  the  change  in  absolute 
vorticity  {C,^)  through  the  life-cycle  (Fig.  5.12).  The  coastal  low  has  relatively  low 
values  of  C,^  at  its  formation,  but  it  spins  up  rapidly  by  the  end  of  the  period.  The 
most  rapid  increase  occurs  between  2000  and  2012  and  is  largest  in  the  upper 
troposphere  (around  350  mb).  Although  the  vorticity  is  relatively  constant  in  the 
vertical  between  the  tropopause  and  850  mb,  the  spin-up  at  the  surface  is  only  about 
half  as  strong  as  that  aloft.  If  favorable  superposition  had  been  achieved  between  the 
PVA  aloft  and  the  surface  center,  the  surface  ^^  might  have  increased  nearly  as  fast  as 
it  did  aloft. 

Examination  of  the  terms  in  the  QLD  vorticity  budget  can  yield  clues  to  the 
processes  that  caused  the  spin-up  in  IOP-1.   The  vorticity  equation  can  be  written 


6(;^;6t    =    -VA_J_^  (V^-  Vq)  dl  -  d!dp  (cot;^)  __         _  (5.2) 

'i;^  (V  •  V)    +    k  •  (dVldp  X  Vco)    +    Fr   +    R, 

where  C,'^  is  absolute  vorticity,  A  is  the  area  of  a  pressure  surface  within  the  budget 
volume,  y^  is  the  normal  wind  component,  V^  is  the  cyclone's  normal  velocity 
component,  co  is  vertical  velocity  and  the  overbar  denotes  an  area-average.  The  term 
on  the  left  side  is  the  quasi- Lagrangian  vorticity  tendency.  On  the  right  side,  the  first 
term  is  the  horizontal  transport  (flux)  of  vorticity  and  the  second  term  is  the  vertical 
flux  convergence.  The  vorticity  divergence  term  (third  term)  is  considered  to  be  a 
source  of  positive  vorticity.  The  fourth  term  is  the  tilting  term  and  the  fifth  term  is 
frictional  dissipation.  The  last  term  on  the  right  is  the  residual,  which  includes  the 
sources  and  sinks  not  explicitly  resolved  by  the  other  terms,  plus  the  computational 
errors.  If  all  the  sources  and  sinks  contributing  to  vorticity  changes  have  been 
properly  accounted  for,  and  interpolation  errors  have  not  been  introduced,  the  residual 
should  be  small  compared  to  the  leading  terms.  Interpolation  and  calculation  errors 
are  the  major  source  of  uncertainty  in  budget  studies  based  on  numerical  model 
output.  For  the  analyses,  an  uncertainty  due  to  incomplete  observations  is  also 
present. 

The  leading  budget  terms  (friction  is  small)  are  averaged  through  the  lower 
troposphere  (1000  to  500  mb)  and  combined  in  Fig.  5.13.  The  leading  terms  causing 
vorticity  spin-up  are  the  divergence  term  and  the  horizontal  transport  of  i^^.    The 
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Figure  5.12.   Area-averaged  absolute  vorticity  over  4°  lat.  radius  in  the  (a)  analysis  and 
,'uTr^^o^ocr  nfthp  cprnndnrv  lov/.    Contours  are  2  x  10"-  s'  . 


;by  forecast  of  the  secondary  lov/.   Contours 
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Figure  5.13.   Average  over  the  lower  troposphere  (1000-500  mb)  of  the  vorticity  budget 
terms  for  the  (a)  analysis  and  (b)  forecast  of  the  secondary'  low.   Time  tendency  (solid), 
lateral  transport  (dashed),  vorticity  divergence  (dotted),  vertical  transport 
(dash-dotted),  tilting  term  (long-short  dashed)  and  residual  (dash-asterisk). 
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vertical  transport  and  the  residual  term  tend  to  decrease  the  lower  tropospheric 
vorticity. 

The  trends,  magnitudes  and  relative  importance  of  the  various  terms  are  quite 
similar  between  the  analysis  and  forecast  budgets.  These  vorticity  budget  calculations 
reinforce  the  earlier  conclusion  that  NORAPS  handled  this  case  well.  The  residual  in 
the  analysis  is  approximately  the  size  of  the  vertical  flux  convergence,  but  smaller  than 
the  horizontal  transport  and  divergence  terms.  In  the  forecast  budget,  the  residual  is 
smaller  than  all  three  leading  terms.  Although  smaller  than  the  leading  terms,  the 
residuals  are  still  relatively  large,  especially  in  the  forecast  budget.  This  error  may  be 
due  to  the  parameterization  of  friction  in  the  budget  formulation.  As  this  low  occurred 
over  land,  it  may  be  that  friction  is  much  larger  than  estimated  in  the  budget. 
Consequently,  a  large  negative  residual  is  required  to  balance  the  vorticity  budget. 

The  divergence  term  is  a  source  (sink)  of  positive  vorticity  below  (above)  the 
USD.  For  low-level  cyclogenesis.  it  is  one  of  the  leading  terms  in  the  spin-up  in  both 
the  analysis  and  forecast  (Fig.  5.14).  Convergence  in  the  lower  levels  concentrates  the 
ambient  (^^  over  a  deep  column.  A  local  maximum  in  the  magnitude  of  the  divergence 
.  term  exists  between  2003  and  2006.  This  maximum  corresponds  to  the  peak  in  the 
vertical  circulation  during  that  period.  After  2009,  the  magnitude  of  the  upper-level 
divergence  decreases  in  both  the  analysis  and  forecast.  The  divergence  term 
contribution  to  C  in  the  lower  troposphere  decreases  in  the  analysis,  but  the  spin-up  in 
the  lowest  layers  continues  until  2015  in  the  forecast.  If  the  superposition  between  the 
surface  center  and  the  upper-level  forcing  had  been  more  favorable,  it  may  be  expected 
that  the  magnitude  of  the  C,^  increases  due  to  the  divergence  term  would  have  been 
significantly  greater. 

3.  Heat  Budget 

The  heat  budget  during  the  life  span  of  a  cyclone  provides  direct  indications  of 
the  physical  properties  that  contribute  to  the  development.  The  area-averaged 
advective  form  of  the  thermodynamic  equation  is 


6T;6t    =    -1/A(\V-Vq)«   VT    -    (O^T/^p     4-     coa/Cp     +     Q.  (5.3) 

where  T  is  temperature,  a  is  specific  volume,  c    is  the  specific  heat  at  constant  pressure 
and  an  overbar  denotes  the  area-average  of  a  term.    The  term  on  the  left  is  the  quasi- 
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Figure  5.14.   Vorticity  divergence  term  within  4°  lat.  radius  in  the  (a)  analysis  and  (b) 
forecast  of  the  secondary'  low.   Solid  (dashed)  contours  indicate  vorticity  increase 
(decrease).    Contours  are  5  x  10'      s*-. 
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Lagrangian  temperature  tendency.  The  first  two  terms  on  the  right  side,  are  the 
horizontal  and  vertical  temperature  advections.  The  third  term  is  the  energy 
conversion  from  ZAPE  to  EAPE  (CE  in  the  energy  budgets  in  Chapter  IV).  Diabatic 
heating  (Q)  cannot  be  directly  calculated  from  the  analysis  fields,  so  it  must  be 
estimated  from  the  budget  residual  (last  term),  which  also  contains  the  computational 
errors  in  the  budget  calculation. 

The  energy  conversion  term  and  the  vertical  advection  term  have  large 
magnitudes  and  opposite  signs.  For  example,  the  energy  conversion  term  reaches 
nearly  40°  C  d"^  in  the  forecast.  These  two  terms  are  often  combined  to  form  an 
adiabatic  cooling  term,  co(r^  -  F),  where  V^  is  the  dry  adiabatic  lapse  rate  and  F  is 
the  environmental  lapse  rate  or  vertical  temperature  gradient.  (F^  -  F)  is  constrained 
to  be  always  positive  in  the  model  atmosphere,  due  to  the  dr>'  convective  adjustment 
scheme.  As  the  mean  vertical  motion  is  upward  in  the  region  of  the  cyclone,  i.e.,  (O  is 
less  than  zero,  this  term  is  negative  and  tends  to  cool  the  column. 

The  adiabatic  cooling  term  and  the  remaining  three  terms  in  (5.3)  are 
vertically  integrated  from  1000  to  500  mb  for  each  time  period  (Fig.  5.15).  The  time 
tendency  of  temperature  based  on  the  analysis  is  negative  throughout  the- period.  That 
is,  the  cooling  effects  of  the  adiabatic  term  override  the  weak  diabatic  heating.  The 
horizontal  advection  of  temperature  is  small  throughout  the  life-cycle  of  the  secondary 
low,  which  is  partly  due  to  the  lack  of  warm  advection  associated  with  the  onshore 
flow  ahead  of  the  system.  If  the  cyclone  had  formed  in  a  strongly  baroclinic  lower 
troposphere  (hypothesized  to  be  an  important  factor  in  rapid  cyclogenesis),  the 
horizontal  thermal  advection  would  probably  have  contributed  to  warming  the  1000  - 
500  mb  column. 

In  the  analysis,  the  dramatic  increase  in  the  vertical  circulation  that  occurred 
between  2003  and  2009  is  reflected  in  an  increase  in  the  magnitude  of  the  adiabatic 
cooling,  but  no  significant  change  in  the  diabatic  heating.  In  the  forecast,  the 
magnitudes  of  both  terms  double  between  1921  and  2003,  then  remain  relatively  large. 
The  vertical  velocity  maximum  in  the  forecast  is  slightly  greater  than  in  the  analysis, 
which  results  in  a  4°C  d'^  increase  in  the  magnitude  of  the  energy  conversion  term  (not 
shown).  Comparisons  of  the  relative  magnitudes  of  the  heat  budget  terms  between  the 
analysis  and  forecast  indicate  that  the  diabatic  heating  and  energy  conversion  terms  are 
slightly  greater  in  the  forecast.  However,  the  changes  compensate  each  other  and  the 
time  tendency  of  temperature  is  nearly  the  same  in  the  forecast  as  in  the  analysis.    The 
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diabatic  heating  term,  which  is  the  primary  source  of  heat  for  both  lows,  is  investigated 
in  more  detail  below. 

Since  diabatic  heating  is  determined  as  a  budget  residual,  all  calculation  and 
interpolation  errors  are  also  included  in  this  term.  Based  on  NORAPS  analyses,  the 
maximum  diabatic  heating  of  the  column  containing  the  coastal  low  was  centered  near 
500  mb  (Fig.  5.16a).  A  maximum  in  heating  at  this  level  does  coincide  with  the 
maximum  in  vertical  circulation  peak  value  at  2006.  The  contribution  to  the  lower- 
troposphere,  column-average  temperature  tendency  is  partially  countered  by  a 
maximum  in  diabatic  cooling  near  the  surface  at  the  same  time.  The  stratospheric 
diabatic  heating  maximum  is  attributed  to  computational  errors.  The  budget  based  on 
the  forecast  of  the  secondary  low  is  quite  similar  to  the  analysis  (Fig.  5.16b).  Diabatic 
heating  is  stronger  (5°C  d'^)  throughout  the  troposphere  after  2003,  but  the  general 
shape  of  the  field  is  the  same.  The  greater  diabatic  heating  during  the  last  half  of  the 
forecast  is  consistent  with  the  differences  found  in  previous  budgets.  Greater  diabatic 
heating  implies  a  larger  vertical  circulation,  greater  low-level  convergence  and  a  larger 
spin-up  of  C,^  in  the  lower  levels  in  the  forecast.  The  height  of  the  diabatic  heating 
maximum  is  lower  (600  mb)  in  the  forecast  and  the  layer  of  diabatic  cooling  near  the 
surface  is  absent  in  the  model. 

The  diabatic  heating  rate  estimated  from  the  heat  budget  residual  for  the 
forecast  can  be  compared  to  the  archived  heating  rates  in  the  model  output  (Fig. 
5.17a).  The  model-predicted  heating  maximum-is  slightly  higher  (500  mb  versus  600 
mb)  than  the  center  diagnosed  as  a  residual  in  the  budget.  The  model  predicts  diabatic 
cooling  below  900  mb,  similar  to  the  analysis  budget.  The  differences  between  the 
budget-diagnosed  diabatic  heating  and  the  model-predicted  "heating  are  generally  small 
(Fig.  5.17b).  Primar\'  areas  of  disagreement  are  near  the  surface,  which  arises  because 
sensible  heating  is  not  included  in  the  archived  heating,  and  in  the  lower  stratosphere, 
which  is  probably  due  to  computational  errors  in  the  budget. 

In  summary,  the  heat  budget  calculated  from  the  NORAPS  forecast  fields 
provides  a  fair  approximation  of  the  analysis  budget.  The  relative  contributions  of  the 
diabatic  and  adiabatic  processes  are  similar  in  both  budgets.  Based  on  the  verification 
of  basic  height  fields  and  the  comparison  of  mass,  vorticity  and  heat  budgets,  it  is 
concluded  that  the  NOR.APS  forecast  from  1900  provided  a  good  description  of  the 
subsequent  development  of  the  important  features  in  IOP-1.  The  moisture  budget 
(following  section)  can  only  be  made  from  the  forecast  fields,  but  it  can  be  assumed 
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Figure  5.16.   Diagnosed  diabatic  heating  within  4°  lat.  radius  in  the  (a)  analysis  and  (b) 
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Figure  5.17.   (a)  Diabalic  heating  from  archives  of  NORAPS  during  the  secondan.-  low 
forecast.    Solid  (dashed)  contours  indicate  heating  (cooling),  (b)  Difference  between 
diasnosed  and  model-predicted  diabatic  heating.   Solid  (dashed)  contours  indicate 
diagnosed  heating  is  larger  (smaller).   Contours  are  5°C  d'  . 
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with  reasonable  confidence  that  it  is  an  accurate  representation  of  the  moisture 
processes  involved  in  the  development  of  both  lows. 

4.  Moisture  Budget 

The  area-averaged  moisture  budget  equation  is  similar  to  the  heat  budget, 

6q;6t    =    -1/A    §    qCVn-V^)   dl  -   dldp(7aq)    +  E-P,  (5.4) 

where  q  is  the  specific  humidity,  E  is  evaporation,  P  is  precipitation  and  an  overbar 
denotes  the  area-average  of  a  term.  The  term  on  the  left  is  the  quasi-Lagrangian 
moisture  tendency.  On  the  right  side,  the  first  two  terms  are  the  horizontal  and  vertical 
moisture  fluxes.  Evaporation  and  precipitation  are  combined  in  the  residual,  or 
moisture  source/sink.  term.  The  sign  convention  for  the  residual  is  negative  (positive)  if 
precipitation  (evaporation)  is  occurring.  This  budget  residual  can  be  compared  to  the 
predicted  moistening  rate  from  the  archived  model  output. 

As  in  the  vorticity  and  heat  budgets,  the  terms  are  integrated  from  1000  to 
500  mb  (Fig.  5..18a).  The  horizontal  flux  moistens  the  lower  troposphere  and  the  . 
condensation/evaporation  term  acts  as  a  sink  (condensation  dominates).  The  vertical 
flux  of  moisture  also  dries  the  lower  troposphere,  but  the  magnitude  of  this  term  is 
relatively  small.  The  lower  troposphere  is  becoming  drier  throughout  the  life-cycle  of 
the  coastal  system  because  the  combination  of  vertical  flux  and  condensation  exceeds 
the  horizontal  influx  of  moisture.  The  horizontal  flux  reaches  a  maximum  at  2003,  but 
quickly  decreases  to  relatively  small  values  during  the  next  period.  This  indicates  that 
the  secondary  low  did  not  continue  to  tap  into  the  moisture  that  presumably  was 
present  in  the  onshore  flow  northeast  of  the  center.  Although  this  may  be  due  to  the 
difluse  nature  of  the  surface  pressure  distribution,  the  moisture  from  the  Atlantic  may 
have  been  too  far  northeast  of  the  surface  center,  and  did  not  provide  an  adequate 
source  of  latent  heat  to  maintain  a  strong  vertical  circulation. 

The  distribution  of  the  moisture  sink  in  the  vertical  (Fig.  5.18b)  coincides  with 
the  location  of  the  maximum  heating  in  the  forecast  as  diagnosed  in  the  heat  budget. 
The  condensation  maximum  is  centered  near  600  mb  at  2006.  The  model-predicted 
moistening  rate  (not  shown)  is  quite  similar  to  the  moisture  sink  diagnosed  from  the 
moisture  budget.  There  is  no  evidence  of  a  stratospheric  maximum  as  there  was  in  the 
diabatic  heating,  which  indicates  that  the  "heating"  at  this  level  was  probably  a  result 
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Figure  5.  IS.   (a)  Average  over  the  lower  troposphere  (1000-500  mb)  of  the  moisture 
budget  terms  for  the  forecast  of  the  secondar\'  low.   Time  tendency  (solid),  moisture 
sink  (dashed),  horizontal  flux  (dotted)  and  vertical  flux  (dot-dashed),   (b)  Moisture  sink 
within  4°  lat.  radius  for  the  secondary'  low  forecast.   Dashed  contours  indicate 
condensation  in  e  ks     d"  . 
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of  errors  in  the  heat  budget,  as  previously  stated.  The  net  moisture  sink  throughout 
the  lowest  levels  does  not  support  a  net  cooling  due  to  the  evaporation  of  rain  below 
850  mb  that  was  indicated  by  the  diabatic  heating  term  in  the  analysis  heat  budget. 
However,  the  moisture  sink  distribution  is  consistent  with  the  presence  of  diabatic 
heating  only  in  the  NORAPS  forecast.  In  summary,  the  XORAPS  forecast  has  more 
heating  and  moistening  of  the  coastal  low  than  is  inferred  from  the  analyses.  It  can 
not  be  determined  whether  this  more  vigorous  circulation  in  the  forecast  is  real  or 
whether  the  weaker  circulation  inferred  from  the  analysis  is  simply  due  to  incomplete 
observations. 

C.       SUMMARY  AND  HYPOTHESIS 

IOP-1  was  a  case  in  which  strong  upper-level  forcing  failed  to  cause  significant 
low-level  cyclogenesis.  The  short  wave  trough  intensified  rapidly  at  upper  levels,  which 
resulted  in  a  very  strong  vorticity  maximum  (28  x  10'^  s"^  at  500  mb).  However,  the 
associated  surface  development  was  poorly  organized.  A  low  west  of  the  Appalachians 
developed  and  filled  within  an  24  -h  period  as  the  center  of  action  shifted  southeast  to  a 
system  that  moved  along  the  Atlantic  seaboard.  It  appears  that  this  coastal  system 
was  close  .to  the  correct  horizontal  scale  to  provide  an  excellent  initial  circulation  for 
rapid  cyclogenesis,  as  in  the  classic  Presidents'  Day  cyclone.  However,  the  timing  and 
superposition  with  the  upper-level  forcing  did  not  appear  to  have  been  favorable,  and 
the  low-level  environment  was  not  conducive  for  intensification.  Baroclinity  along  the 
coast  was  weak  and  the  static  stability  was  relatively  high. 

There  are  many  factors  in  this  case  that  appear  to  have  not  been  favorable  for 
rapid  cyclogenesis.  In  addition  to  improper  superposition  and  a  relatively  stable  low- 
level  atmosphere,  the  location  of  the  short  wave  over  land  may  have  detracted  from  the 
development  of  the  surface  cyclone.  If  the  surface  low  had  been  located  over  the 
ocean,  it  would  have  been  subjected  to  smaller  surface  friction  and  had  a  direct  access 
to  the  energy  available  from  the  relatively  warm  water  via  fluxes  of  sensible  and  latent 
heat.  From  the  point-of-view  of  the  preceding  chapters,  the  coastal  low  developed  too 
late  and  too  far  southeast  of  the  upper-level  vorticity  source.  Therefore,  the  upper- 
level  forcing  induced  development  of  a  low  to  the  northwest  of  the  coastal  low.  When 
the  coastal  low  and  the  upper  feature  did  move  to  a  position  where  mutual 
reinforcement  was  possible,  the  coastal  system  was  too  weak  to  achieve  the  maximum 
deepening  that  would  have  been  possible  if  the  lower-level  system  would  have  already 
had  significant  vorticitv. 
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A  more  favorable  superposition  may  have  been  achieved  if  the  low-level 
baroclinity  along  the  Atlantic  coast  had  been  larger  on  19  January-  (i.e.,  a  coastal  front 
or  cold-air  damming  against  the  eastern  slope  of  the  Appalachians).  The  presence  of  a 
strong  thermal  gradient  would  have  encouraged  more  rapid  development  of  the  coastal 
low,  which  might  have  made  it  strong  enough  to  better  link  with  the  shortwave  aloft 
when  it  moved  into  a  more  favorable  superposition  around  00  GMT  20  January. 
Based  on  the  conclusions  in  the  "baroclinic"  case  in  Chapter  IV,  a  stronger  low-level 
thermal  gradient  should  shift  the  position  of  the  most  favorable  surface  development  of 
the  primar>'  low  equatorward  of  the  jet  core.  The  combination  of  a  stronger  coastal 
low  and  an  equatorward  shift  in  the  most  favorable  surface  position  may  result  in  a 
superposition  that  increases  the  surface  development.  Specifically,  the  energy 
conversion  and  diabatic  heating  terms  in  the  heat  budget  should  increase,  along  with 
the  mass  circulation.  This  hypothesis  will  be  tested  in  the  following  section  by 
artificially  increasing  the  low-level  baroclinity  along  the  mid-Atlantic  coast  to  provide 
modified  initial  conditions  for  NORAPS.  If  the  superposition  conclusions  are  correct, 
the  primary  low  should  shift  to  the  southeast,  the  secondary'  low  should  dominate  more 
quickly,  and  the  deepening  rate  of  the  coastal  cyclone  should  increase. 

D.       MODIFIED  FORECAST 

To  test  the  above  hypothesis,  the  initial  conditions  at  1900  are  adjusted  to  make 
the  lower  troposphere  less  stable  along  the  Atlantic  seaboard.  This  is  accomplished  by 
inserting  a  low-level  jet  (LLJ)  parallel  to  the  coast  form  Charleston,  SC  to  New  York. 
To  balance  the  changed  winds,  the  mass  field  is  altered  in  the  lowest  levels.  While  the 
forecast  from  the  new  initial  conditions  is  quite  similar  to  the  operational  NORAPS 
forecast,  several  significant  changes  support  the  hypothesis  in  the  previous  section. 

1.  Numerical  Procedure 

Analyzed  winds  between  700  and  1000  mb  from  Charleston  to  New  York  are 
replaced  by  a  LLJ  from  the  southwest.  The  jet  width  extends  from  the  Appalachians 
to  approximately  200  km  seaward  of  the  continent.  A  300  km-wide  "blend  zone"  on 
either  side  of  the  jet  allows  a  smooth  transition  to  the  analyzed  fields.  Divergence  is 
removed  from  the  new  wind  pattern  and  the  mass  field  is  adjusted  on  pressure  levels  to 
balance  the  non-divergent  winds. 

The  operational  NORAPS  initialization  scheme  removes  divergence,  but  does 
not  balance  the  mass  field  on  pressure  surfaces  as  was  done  for  the  modified  forecast. 
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To  insure  that  this  step  does  not  introduce  significant  changes  in  the  forecast,  the 
operational  forecast  is  re-created  with  mass  balancing  included.  Forecast  fields  at  all 
levels  are  nearly  indistinguishable  from  the  forecasts  discussed  in  section  B.  The 
forecast  minimum  pressure  (at  2018)  is  1  mb  higher  in  the  balanced  forecast  and  the 
maximum  deepening  rate  is  slightly  smaller.  Thus,  any  changes  induced  by  the 
initialization  scheme  that  includes  mass  balancing  would  tend  to  reduce  the  predicted 
intensification  rate. 

2.  Modified  Initial  Conditions 

Insertion  of  a  southwesterly  LLJ  along  the  east  coast  and  removal  of 
divergence  results  in  small  changes  in  the  winds  elsewhere  in  the  troposphere  (Fig. 
5.19).  The  polar  front  and  subtropical  jet  maxima  are  in  the  same  locations,  but  are 
somewhat  (3  to  7ms")  weaker.  Additionally,  a  local  minimum  near  500  mb  over 
Wheeling,  WV  is  more  pronounced.  The  LLJ  maximum  is  32  m  s"^  at  620  mb  over 
Cape  Hatteras.  The  strong  vertical  wind  shear  below  the  LLJ  indicates  a  lower 
Richardson  number  (R-)  between  Richmond  and  Cape  Hatteras.  To  balance  the 
increased  winds,  the  boundary'  layer  temperature  (not  shown)  is  approximately  4'-*C 
cooler  between  Richmond  and  Cape  Hatteras  than  in  the  analysis  (Fig.  5.2).  This 
weak,  cold  dome  simulates  conditions  that  may  result  from  moderate  cold-air  damming 
against  the  eastern  slopes  of  the  Appalachians.  These  revised  initial  conditions  should 
be  more  favorable  for  development  of  the  coastal  low.  Also,  'the  results  from  the 
"baroclinic"  case  in  Chapter  IV  indicate  that  increasing  the  low-level  shear  moves  the 
most  favorable  position  of  the  lower  perturbation  toward  the  equatorward  side  of  the 
jet.  Thus,  these  changes  in  the  initial  conditions  should  move  the  primary  low.  which 
is  induced  by  the  vigorous  upper-level  forcing,  toward  the  anticyclonic  side  of  the 
upper-level  jet. 

The  reduced  stability  is  evident  in  the  lower-tropospheric  static  stability  at 
1900  (Fig.  5.20a).  At  this  time,  the  coastal  low  is  forming  near  the  NC/SC  border 
(asterisk  in  Fig.  5.20b).  The  stability  in  the  modified  initial  conditions  is  lower  than  in 
the  operational  analysis  throughout  the  eastern  U.S.  (Fig.  5.4). 

These  modifications  are  expected  to  induce  changes  in  the  forecast  that  will 
lead  to  a  more  favorable  superposition.  Results  from  previous  chapters  indicate  that 
superposition  is  most  important  when  strong  features  exist  at  both  levels.  Although  a 
large  low-level  vorticity  center  is  not  added  in  this  case,  more  vigorous  development  of 
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Figure  5.19.   As  in  Fig.  5.2.  except  for  the  modified  initial  conditions  at  00  GMT  19 
Januan.-  1986. 
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the  coastal  low  is  expected  to  create  a  surface  center  that  interacts  more  strongly  with 
upper-level  forcing. 

3.  Changes  in  the  Forecast 

The  influence  of  the  modified  initial  conditions  is  evident  at  all  levels  and  at  all 
times  of  the  forecast.  In  the  lower  troposphere,  the  position  and  relative  intensification 
rates  of  the  primary  and  secondary  lows  in  the  modified  forecast  change  as  expected. 
The  track  of  the  primary-  low  is  displaced  to  the  southeast  compared  to  the  position  in 
the  operational  forecast.  At  1912,  the  lowest  surface  pressure  is  165  km  to  the  east- 
southeast  (Fig.  5.20b).  The  lower  stability  and  strong  low-level  shear  also  contribute 
to  more  rapid  development  of  the  coastal  low.  Cyclogenesis  occurs  approximately  6  h 
earlier,  and  the  secondary  low  becomes  the  dominant  system  at  2000  (compared  to 
2006  in  the  operational  forecast). 

The  cyclogenesis  in  the  middle  and  upper  troposphere  in  the  modified  forecast 
is  slightly  stronger  than  in  the  operational  forecast.  The  36  h  modified  forecast  (valid 
at  2012)  is  shown  in  Fig.  5.21  for  comparison  with  the  analysis  (Fig.  5.6)  and 
operational  forecast  (Fig.  5.10).  The  250  mb  trough  is  approximately  80  meters  deeper 
than  in  the  operational  forecast  (Fig.  5.2 Id),  and  is  actually  closer  to  the  analyzed 
height  pattern  at  2012.  The  modified  500  mb  low  (80  m  deeper  than  the  operational 
forecast)  is  also  quite  similar  to  the  2012  analysis.  Since  the  upper- level  trough  is 
deeper,  the  height  gradient  is  larger,  and  the  jet-level  wind  maxima  are  approximately  S 
to  10  m  s'^  stronger  than  in  the  operational  forecast. 

The  maximum  12  h  deepening  of  the  coastal  low  is  14  mb  in  the  modified 
forecast,  which  is  larger  than  the  rate  in  the  operational  run  (11  mb)  or  in  the  forecast 
with  balancing  included  (9  mb). .  However,  the  rapid  deepening  in  the  modified  forecast 
does  not  continue,  and  the  minimum  sea-level  pressure  of  989  mb  is  only  1  mb  deeper 
than  in  the  operational  forecast  (Fig.  5.21c).  Thus,  the  changes  induced  in  the  low- 
level  circulation  are  in  the  desired  direction,  but  are  not  sufficient  in  magnitude  to 
achieve  the  most  favorable  superposition  between  the  coastal  low  and  the  strong 
upper-level  positive  vorticity.  Despite  the  southeastward  shift  in  the  position  of  the 
primary  low,  the  surface  centers  are  still  separated  by  more  than  500  km.  Thus,  even 
after  modification  of  the  initial  conditions,  the  along-stream  and  cross-stream  distance 
between  the  upper-level  vorticity  center  and  the  coastal  low  is  too  great  to  achieve  the 
most  rapid  deepening.    For  an  optimum  superposition  to  occur  in  IOP-1.  the  entire 
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Figure  5.21.   (a,  h  and  c)  As  in  Fig.  5.1.  except  for  the  modified  36  h  forecast  at  12 
GVIT  20  Januar\-  19S6.   (d,  e  and  0  DifTerence  between  modified  and  operational 
forecasts  (operational  minus  modified)  at  (d)  250  mb.  (e)  500  mb  and  (0  surface. 
Contours  in  (d)  and  (e)  are  20  gpm  and  in  (0  are  2  mb. 
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middle  and  upper  troposphere  would  have  to  be  shifted  eastward  several  hundred  km 
with  respect  to  the  surface  low. 

E.       BUDGETS  FOR  THE  MODIFIED  FORECAST 

Quasi- Lagrangian  diagnostic  budgets  are  performed  on  the  forecast  from  the 
modified  1900  conditions  for  comparison  with  the  budgets  based  on  the  operational 
forecast.  Since  the  modified  forecast  is  quite  similar  to  the  original  forecast,  it  is 
expected  that  the  modified  budgets  will  also  be  very  much  like  the  forecast  budgets  of 
the  previous  sections. 

1.  Mass  and  Vorticity  Budgets 

The  mass  circulation  in  the  modified  forecast  of  the  coastal  low  (Fig.  5.22)  is 
very  similar  to  the  operational  forecast  (Fig.  5.11b).  The  magnitudes  of  the  low-level 
convergence  and  jet-level  divergence  maxima  are  nearly  unchanged.  The  major  change 
in  the  mass  circulation"  is  that  the  strongest  circulation  occurs  approximately  6  h  earlier 
in  the  modified  forecast.  This  earUer  development  in  the  circulation  is  consistent  with 
the  earlier  formation  of  the  coastal  low  (see  Fig.  5.11a),  and  will  be  evident  in  each 
budget.  The  mass  circulation  in  the  modified  forecast  also  decreases  more  rapidly  at 
the  end  of  the  forecast  period.  By  2018,  the  magnitude  of  the  mass  outflow  (Fig.  5.22) 
is  approximately  50  percent  of  that  in  the  operational  forecast  (Fig.  5.11b).  The  low- 
level  mass  convergence  is  also  smaller  than  in  the  operational  forecast.  By  the  end  of 
the  modified  forecast,  the  mass  circulation  is  minimal.  Weak  convergence  remains 
below  800  mb,  but  divergence  aloft  is  nearly  zero.  This  rapid  weakening  of  the  vertical 
circulation  occurs  as  the  upper-level  low  moves  directly  over  the  surface  center,  which 
reflects  the  lack  of  favorable  superposition  between  the  levels. 

The  earlier  formation  is  also  evident  in  the  time-height  cross-section  of  'C,^ 
(Fig.  5.23).  The  vorticity  in  the  modified  forecast  not  only  peaks  earlier,  but  it  also 
begins  to  decay  by  the  end  of  the  period.  Although  the  upper-level  vorticity  evolution 
in  the  modified  forecast  is  closer  to  the  analysis  (Fig.  5.12a),  the  desired  conditions  for 
a  more  favorable  superposition  between  the  upper-level  and  surface  were  not  achieved. 
The  earlier  spinup  and  subsequent  decrease  in  the  cyclogenesis  can  also  be  seen  in  the 
lower-troposphere  average  vorticity  budget  (Fig.  5.24).  The  terms  have  the  same 
general  shape  and  magnitude  as  they  had  in  the  operational  forecast  budget,  but  the 
maximum  values  occur  approximately  6  h  earlier  (as  in  the  analysis)  and  the 
magnitudes  decrease  rapidly  after  2015.    The  vorticity  divergence  and  lateral  transport 
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Figure  5.22.   Horizontal  mass  flux  across  the  4°  lat.  ra"dius  in  the  modified  forecast  of 
the  secondary'  low.    Solid  (dashed)  contours  denote  outflow  from  (inflow  into)  the  4° 
lat.  radius  budget  volume.   Contours  are  5  x  10'^^  g  s'^  mb*^ 


terms  decrease  so  much  that  the  quasi- Lagrangian  time  tendency  of  C^  actually 
becomes  negative  after  2015  as  the  lower- level  system  begins  to  spin  down.  The  spin- 
down  occurs  as  the  system  becomes  vertically  stacked  from  the  surface  to  the 
tropopause  at  2018  (not  shown).  The  decrease  in  i^^  generation  by  the  divergence  term 
is  consistent  with  the  weaker  mass  -circulation  at  the  end  of  the  forecast. 

2.  Heat  and  Moisture  Budgets 

Indications  of  an  earlier  intensification  and  earlier  decay  are  also  found  in  the 
column-averaged  budgets  of  heat  and  moisture  for  the  coastal  low  (Fig.  5.25).  The 
relative  sizes  and  maximum  magnitudes  of  all  the  budget  terms  are  nearly  the  same  as 
in  the  operational  forecast  (Figs.  5.15b  and  5.18a),  except  that  the  maxima  are 
achieved  approximately  6  h  earlier  and  the  magnitudes  are  smaller  by  the  end  of  the 
period.  Although  the  amplitudes  are  larger  and  are  achieved  earlier,  the  levels  of 
maximum  heating  and  dr>'ing  (Fig.  5.26)  are  nearly  identical  to  the  operational  forecast 
(Figs.  5.16b  and  5.18b). 

136 


888- 


962^ 


1912 


1918  2000 


2006 
TIME 


2012  2018 


2100 


Fisure  5.23.   Area-averaeed  absolute  vorticity  over  4°  lat.  radius  in  the  modified 
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Figure  5.24.   Average  over  the  lower  troposphere  (1000-500  mb)  of  the  vorticity  budget 
terms  for  the  modified  forecast  of  the  secondar\-  low.   Time  tendency  (solid),  lateral 
transport  (dashed),  vorticity  divergence  (dotted),  vertical  transport  (dash-dotted),  tilting 
term  (long-short  dashed)  and  residual  (dash-asterisk). 
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Figure  5.25.   Average  over  the  lower  tropospheric  (1000-500  mb)  of  the  budgets  of  (a) 
heat  and  (b)  moisture  for  the  modified  forecast  of  the  secondan."  low. 
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Figure  5.26.   (a)  Diabatic  heating  (°C  d"^)  and  (b)  moisture  sink  {g  kg"^  d"^)  within  4° 
lat.   radius  in  the  modified  forecast  of  the  secondare'  low. 
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F.       SUMMARY  OF  IOP-1  CASE  STUDY 

IOP-1  was  a  case  in  which  a  vigorous  upper-level  PVA  source  was  not  favorably 
superposed  with  respect  to  the  surface  lows  to  create  the  conditions  for  rapid 
cyclogenesis.  Additionally,  the  low-level  environment  was  not  favorable  for 
cyclogenesis,  as  the  stability  was  too  high  and  low-level  baroclinity  was  ver>'  weak.  As 
a  result,  intensification  of  the  surface  low  was  minimal.  The  sea-level  pressure  pattern 
was  broad  and  diffuse,  with  centers  co-existing  west  of  the  Appalachians  and  along  the 
Atlantic  coast  for  much  of  the  period.  Although  the  coastal  low  formed  in  a  position 
near  the  NC  coast  that  is  often  favorable  for  subsequent  rapid  deepening,  it  was  well  to 
the  east  of  the  most  favorable  surface  position  with  respect  to  the  upper-level  forcing. 
As  in  the  cases  discussed  in  Chapters  III  and  IV,  deepening  does  occur  with  these 
larger  vertical  tilts,  but  it  is  not  as  rapid. 

The  NORAPS  forecast  from  00  GMT  19  January  is  quite  credible  through  36  h. 
The  track,  intensification  and  rapid  filling  of  the  primary  low  center  west  of  the 
Appalachians  are  accurately  predicted,  and  the  genesis  and  deepening  of  the  coastal 
low  are  well-represented.  This  provides  justification  for  also  calculating  the  budgets  for 
the  forecast  lows  to  compare  with  the  budgets  based  on  the  analyses. 

Quasi-Lagrangian  diagnostic  budgets  from  the  NORAPS  analyses  and  forecast 
illustrated  the  lack  of  favorable  coupling  between  the  upper  and  lower  levels.  The  low 
west  of  the  Appalachians  was  induced  by  the  strong  upper-level  processes,  but 
developed  slowly,  and  failed  to  establish  a  significant  vertical  circulation  before  the 
upper  center  overtook  the  surface  circulation  and  moved  to  the  east.  The  resulting 
easterly  tilt  with  height  (unfavorable  for  intensification)  caused  rapid  weakening  of  the 
primar\'  low  after  2000. 

The  coastal  low  formed  to  the  southeast  of  the  primary  center  as  the  primary  low 
reached  maximum  intensity.  As  a  result,  the  low-level  convergence  and  positive 
vorticity  were  dispersed  over  a  wide  area,  and  the  coastal  cyclone  was  still  relatively 
weak  as  the  upper-level  center  moved  into  a  favorable  superposition  around  2000. 
Consequently,  the  maxima  in  the  mass  circulation,  vorticity  generation,  diabatic 
heating  and  precipitation,  which  were  observ'ed  between  2000  and  2006,  were  not 
maintained,  and  intensification  of  the  system  ceased. 

Modifications  to  the  initial  conditions  for  the  NORAPS  forecast  were  made  to 
simulate  earlier  development  of  the  coastal  low,  and  to  move  the  most  favorable 
position  with  respect  to   the  upper  forcing   toward  the  coast.    The   modified  initial 
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conditions  simulated  the  cold-air  damming  and  low-level  jet  that  evidently  contributed 
to  the  classical  cyclogenesis  in  the  Presidents'  Day  case.  Although  both  low  centers 
reacted  as  expected  to  the  modified  conditions,  and  the  intensification  rate  of  the 
coastal  system  was  measurably  increased,  the  modifications  were  not  sufiicient  to  bring 
about  a  truly  favorable  superposition  or  to  sustain  rapid  deepening  for  more  than  a 
short  time.  Since  creation  of  more  favorable  low-level  environmental  conditions  along 
the  Atlantic  coast  was  not  sufficient,  more  drastic  modifications  to  reduce  the  vertical 
tilt  between  the  upper  and  lower  circulations  would  be  necessary  to  achieve  optimum 
conditions. 
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VI.  SUMMARY  AND  CONCLUSIONS 

Rapid  intensification  of  mid-latitude  cyclones  requires  a  synergistic  interaction  of 
factors  on  several  scales  and  at  all  levels  of  the  troposphere.  A  favorable  superposition 
of  PVA  aloft  over  a  pre-existing  surface  low  facilitates  this  interaction.  Whereas 
cyclogenesis  can  occur  when  small  perturbations  are  present  in  a  baroclinically 
unstable  mean  flow,  it  is  hypothesized  that  rapid  cyclogenesis  begins  with  the 
superposition  of  relatively  large  perturbations.  Therefore,  this  study  investigates  this 
favorable  superposition  between  fmite-amplitude  perturbations  in  the  upper  and  lower 
troposphere.  The  horizontal  scales  of  the  perturbations  are  the  same  at  both  levels  and 
are  equal  to  the  most  unstable  wavelength  of  the  baroclinically  unstable  mean  flow. 

Experiments  with  idealized  models  (a  two-level  analytical  model  and  an  adiabatic, 
frictionless  channel  version  of  NORAPS)  demonstrate  that  the  initial  alongstream 
phase  difference  between  developing  features  is  crucial  when  both  the  upper  and  lower- 
level  perturbations  are  relatively  strong.  When  either  perturbation  is  significantly 
stronger  than  the  other,  the  phase  difference  has  less  effect  on  subsequent  deepening. 
The  dominant  perturbation  in  that  case  simply  forces  growth  at  the  other  level  in  a 
position  that  is  favorable  for  intensification,  and  the  (small)  pre-existing  perturbation 
at  that  level  is  of  little  consequence.  However,  the  most  rapid  growth  for  given  mean 
flow  and  static  stability  conditions  occurs  when  relatively  strong  perturbations  at  both 
levels  are  favorably  superposed.  The  most  favorable  phase  difference  (a  westward  tilt 
with  height  of  slightly  less  than  a  quarter-wavelength)  is  consistent  with  synoptic 
experience.  Westward  tilts  other  than  the  most  favorable  also  result  in  growth  of  the 
perturbations,  but  the  growth  rate  is  smaller.  By  contrast,  when  the  initial  tilt  is 
eastward  with  height,  the  perturbations  negatively  interact  and  the  amplitudes  decrease 
at  both  levels  until  the  perturbations  shift  to  a  westward-tilt  orientation. 

The  analytical  model  can  only  address  the  along-stream  tilt  between  the 
perturbations.  Thus,  the  simplified  version  of  NORAPS  is  used  to  test  variations  in 
the  north-south  position  of  the  surface  low.  The  cross-stream  tilts  only  modify  the 
optimum  growth  rates  that  is  determined  by  the  zonal  superposition.  If  the  zonal  tilt  is 
near  the  optimum  westward  tilt  with  height  (approximately  70°  between  400  mb  and 
the  surface),  the  perturbations  grow  rapidly  regardless  of  their  initial  placement  north 
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or  south  of  mid-channel.  Each  case  has  a  most  favorable  cross-track  position,  which  is 
determined  by  the  mean  flow  conditions.  If  the  low-level  environment  is  highly 
baroclinic  (vertical  wind  shear  is  large  in  the  lowest  layers),  the  most  favorable  pre- 
existing disturbance  position  shifts  toward  the  south.  If  the  horizontal  wind  shear  in 
the  middle  and  upper  troposphere  is  enhanced  (by  increasing  the  zonal  mean  velocity 
of  the  jet  stream),  the  most  favorable  disturbance  position  shifts  toward  the  north. 
Regardless  of  the  initial  disturbance  position,  the  developing  surface  center  shifts  in 
time  toward  the  most  favorable  cross-track  position  as  it  intensifies. 

When  the  meridional  tilt  is  varied,  the  relative  growlh  rates  of  the  surface 
disturbances  change  with  time.  The  perturbation  that  is  nearest  the  most  favorable 
position  for  a  given  mean  flow  has  the  largest  baroclinic  energy  conversions  after  24  to 
36  h,  but  the  growth  rate  early  in  the  integration  depends  directly  on  the  position  of 
the  disturbance  relative  to  the  zone  of  low-level  surface  baroclinity.  The  nearer  the 
perturbation  is  positioned  to  the  largest  surface  temperature  gradient,  the  larger  the 
growth  early  in  the  integration. 

When  the  upper-level  wave  perturbation  is  replaced  as  the  source  of  PVA  aloft 
by  a  jet-streak,  the  conclusions  with  respect  to  superposition  are  unchanged.  The 
growth  rates  of  the  surface  perturbations  depend  directly  on  the  initial  zonal  phase 
difference  between  the  streak  "and  the  surface  trough.  As  in  the  experiments  with  an 
upper-level  wave,  the  perturbation  that  has  the  same  meridional  position  as  the 
strongest  low-level  baroclinic  zone  experiences  the  most  rapid  growth  at  the  outset  of 
the  integration.  However,  the  cross-stream  position  is  less  important  when  a  jet  streak 
provides  strong  upper-level  forcing,  because  the  streak  effectively  shifts  the  surface 
perturbation  to  a  favorable  position  (under  the  left-exit  quadrant)  early  in  the 
integration,  regardless  of  where  it  is  initially  positioned. 

Since  these  experiments  are  conducted  with  simplified  (linear,  adiabatic, 
frictionless)  models,  caution  must  be  used  in  applying  the  results  to  the  rapid 
cyclogenesis  hypothesis.  However,  the  results  appear  to  be  consistent  with  synoptic 
experience. 

An  actual  case  of  east  coast  cyclogenesis  that  occurred  during  lOP-l  of  GALE 
(19-21  January'  1986)  is  investigated  with  the  conclusions  from  the  two-level  and 
NORAPS  channel  model  studies  in  mind.  Rapid  intensification  did  not  occur  at  the 
surface  in  this  case,  even  though  a  dramatic  upper-level  cyclogenesis  did  develop.  Two 
surface  lows  did  develop:    one  west  of  the  Appalachians  associated  with  a  strong  short 
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wave  and  PVA  at  upper  levels,  and  the  second  along  the  Carolina  coast.  Although  the 
coastal  low  subsequently  dominated,  it  did  not  establish  a  well-developed  vertical  link 
with  the  upper-level  forcing,  and  intensification  was  relatively  weak  and  short-lived. 
Although  this  absence  of  a  rapid  surface  cyclogenesis  may  have  been  due  to  a 
combination  of  factors,  the  hypothesis  tested  here  is  that  it  was  due  to  a  lack  of 
favorable  superposition. 

By  modifying  the  initial  conditions  in  the  IOP-1  case,  the  low-level  environment 
in  the  region  of  the  Atlantic  coast  was  made  more  favorable  for  cyclogenesis.  The 
forecast  from  the  modified  initial  conditions  included  earlier  development  of  the  coastal 
low  and  a  shift  of  the  Appalachian  low  toward  the  anticyclonic  side  of  the  upper-level 
jet.  Although  the  intensification  rate  of  the  coastal  low  was  increased  compared  to  the 
operational  forecast,  a  truly  favorable  superposition  with  the  upper-level  system  was 
not  achieved.  Although  the  separation  between  the  upper-level  vorticity  center  and  the 
surface  low  in  this  case  is  too  large  to  completely  simulate  favorable  superposition,  the 
modified  forecast  supports  the  conclusions  from  the  simplified  model  experiments.  The 
southeastward  shift  in  the  track  of  the  Appalachian  low  indicates  that  the  most 
favorable  position  for  a  surface  disturbance  moves  equaton^'ard  in  response  to  the 
increased  vertical  wind  shear  in  the  lowest  levels.  By  forcing  earlier  growth  of  the 
coastal  low,  the  lower  feature  is  stronger  in  the  modified  forecast  when  the  upper-level 
forcing  moves  eastward  into  a  slightly  more  favorable  superposition.  Consequently, 
the  intensification  rate  in  the  modified  forecast  is  significantly  larger. 

This  study  isolates  the  effects  of  variations  in  the  superposition  between  existing 
upper-level  and  lower-level  disturbances.  It  demonstrates  that  when  vigorous  upper- 
level  forcing  interacts  with  a  surface  cyclone  of  the  same  horizontal  scale,  the  relative 
position  between  the  two  is  critical  to  subsequent  development  at  both  levels. 
Therefore,  rapid  intensification  not  only  depends  on  the  strength  of  upper-level  and 
surface  features,  and  the  favorability  of  the  low-level  environment,  but  also  on  the 
proper  superposition  of  these  prerequisites. 
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APPENDIX  A 
TWO-LEVEL  MODEL  STANDARD  ATMOSPHERE  CASES 

1.        BAROCLINIC  INSTABILITY  IN  A  STANDARD  ATMOSPHERE 

The  two-level  model  is  tested  in  a  standard  atmosphere  to  insure  the  validity  of 
the  method  used  in  the  low  stability  cases  of  Chapter  III.  For  a  static  stability  which 
corresponds  to  a  lapse  rate  of  6.25  K  km"^  and  a  mean  flow  shear  of  20  ms'^  the 
growth  rate  curve  is  given  in  Fig.  3.2.  The  maximum  growth  rate  (2.38  day"^)  occurs 
for  a  wavelength  of  4000  km,  consistent  with  expectations  from  baroclinic  instability , 
theory.  The  shortwave  cutofl*  for  this  static  stability  is  at  2591  km.  Perturbations 
shorter  than  this  will  be  neutral  rather  than  unstable. 

The  first  standard  atmosphere  experiment  (Case  IS)  is  made  with  a  wavelength  of 
4000  km,  an  initial  streamfunction  perturbation  at  level  1  of  10  m^s'^  and  a  level  3 
perturbation  of  1  m^s'^  The  second,  Case  2S,  is  made  with  the  same,  perturbation 
amplitudes,  but  with  a  wavelength  of  2600  km.  The  amplitudes  and  phase 
relationships  of  the  growing  tipper  and  lower-level  waves  are  shown  in  Fig.  A.l. 
Positive  phase  differences  indicate  a  wave  which  tilts  westward  with  height. 

The  Case  IS  upper-level  wave  grows  from  10  m^s'^  to  more  than  30  m^s"^  in  48, 
hours.  The  corresponding  lower  perturbation,  which  is  initially  an  order  of  magnitude 
smaller,  grows  very  rapidly  and  becomes  nearly  as  large  as  the  upper  after  48  hours. 
The  growth  of  the  vertically-averaged  perturbation  from  5.5  m~s"^  to  29.4  m^s'^ 
coincides  with  the  expected  growth  rate.  The  perturbations,  which  were  initially  in 
phase,  shift  so  that  a  westward  tilt  with  height  (65*^)  is  established. 

Case  2S  uses  a  perturbation  wavelength  just  slightly  longer  than  the  shortwave 
cutoff.  The  reduced  instability  is  reflected  in  the  lower  growth  rate,  but  the  800  mb 
perturbation  again  grows  to  the  same  order  of  magnitude  as  that  at  400  mb  in  48 
hours.  As  explained  by  Hoskins  er  al.  (1985),  the  westward  tilt  for  this  shorter  wave  is 
not  as  great  as  in  Case  IS.  A  55°  tilt  is  established  after  8  hours,  but  it  quickly 
decreases  to  less  than  30°. 


146 


CASES  1S  and  2S 

Perturbation  Amplitude 


Phase  Difference 


90  n 


Case  1S 


16  24 

Time  (h) 


— 1 — 

32 


Case  28 


b 


40 


48 


Figure  A.l.   (a)  Amplitudes  of  400  mb  and  800  mb  perturbations,  and  (b)  phase 
difTerences  through  48  h  for  Cases  IS  (solid)  and  2S^ (dashed).    Initial  streamfunction 
perturbations  in  fa)  are  10  m-s*^  at  level  1  and  1  nrs"^  at  level  3. 
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2.  THE  EFFECT  OF  VARIATIONS  IN  THE  MEAN  FLOW  SHEAR 

Cases  3S  and  4S  investigate  the  importance  of  the  shear  in  the  mean  flow.  In 
both  cases,  the  amplitude  and  wavelength  of  the  initial  perturbations  and  the  static 
stability  are  identical  to  Case  IS.  In  Case  3S  the  shear  between  levels  1  and  3  is 
10  ms'^,  or  half  that  of  Case  IS.  Comparing  Fig.  A. 2a  to  Fig.  A. la,  it  can  be  easily 
seen  that  the  amplitudes  of  both  the  upper  and  lower-level  growing  waves  are  much 
smaller  when  the  mean  shear  is  reduced.  Conversely,  if  the  shear  is  doubled  (Case  4S), 
the  waves  grow  so  fast  that  by  48  hours  the  magnitudes  of  the  perturbations  are  nearly 
six  times  as  large  as  those  found  in  Case  IS.  Indeed,  when  the  shear  is  40  ms"^  it 
takes  only  24  hours  to  grow  a  wave  which  is  as  large  as  that  in  Case  IS. 

The  basic  patterns  of  the  growing  wave  are  similar,  regardless  of  the  shear  in  the 
mean  flow.  The  wave  tilts  westward  with  height  ("against"  the  shear)  and  the  phase 
separation  between  the  400  mb  and  800  mb  features  is  approximately  65°  (Fig.  A. 2b). 
The  magnitude  of  the  shear  only  affects  the  growth  rate,  as  seen  in  equation  (3.14). 
The  (constant)  shear  in  the  mean  flow  between  level  3  and  level  1  is  directly  related  to 
the  meridional  temperature  gradient,  dlidy,  by  the  thermal  wind  equation.  A  larger 
shear  indicates  a  larger*  meridional  temperature  gradient,  and  thus  a  greater  source  of 
APE.   Thus,  doubling  the  shear  will  result  in  a  two-fold  increase  in  the  growth  rate. 

3.  SHORT  NEUTRAL  WAVES 

The  development  of  waves  which  are  too  short  to  be  baroclinically  unstable  is 
also  investigated  using  the  two-level  model.  Initial  conditions  are  the  same  as  in  Case 
IS  (static  stability  corresponding  to  a  lapse  rate  of  6.25  K  km'',  U^  =  10  m  s"^ 
\|/j  =  10  m-^s'^  y-j  =  1  m-^s'^).  Perturbation  wavelengths  shorter  than  the  shortwave 
cutoff  (2591  km)  are  tested  and  96  hour  level  3  growth  and  phase  differences  are  shown 
in  Fig.  A. 3. 

The  neutral  waves  do  not  maintain  a  favorable  vertical  connection  between  the 
perturbations.  The  phase  continually  changes  (Fig.  A. 3b),  and  the  waves  grow  when 
the  tilt  is  westward  with  height.  As  soon  as  the  phase  difference  changes  sign,  the 
amplitude  of  the  wave  begins  to  decay.  The  closer  the  wavelength  is  to  the  critical 
value,  the  longer  the  initial  growth  period  (Fig.  A. 3b)  and  the  larger  the  maximum 
amplitude  (Fig.  A. 3a).  For  a  2000  km  wave,  an  oscillation  between  growth  and  decay 
is  established. 
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For  waves  which  are  just  shorter  than  the  shortwave  cutofT,  the  period  of  gro'^th 
can  be  quite  long.  The  2500  km  wave  does  not  shift  to  a  decaying  phase  until  92 
hours.  Naturally,  linear  model  assumptions  break,  down  long  before  this.  Still,  the 
two-level  model  indicates  that  waves  just  shorter  than  the  shortwave  cutoff  can  achieve 
significant  early  grovvih. 
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Figure  A. 2.   As  in  x^igure  A.l,  except  for  Cases  3S  (solid)  and  4S  (dashed). 
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Fieure  A. 3.   (a)  Amplitudes  ofSOO  mb  perturbations,  and  (b)  phase  difTerences  through 
96^h  for  wavelengths  (L)  less  than  the  shortwave  cutoff.    L  =  2000  km  (dotted),  2200 
km  (dot-dashed)r2400  km  (dashed)  and  2590  km  (solid). 
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APPENDIX  B 

NAVY  OPERATIONAL  REGIONAL  ATMOSPHERIC  PREDICTION 

SYSTEM-NORAPS 

I.       MODEL  CHARACTERISTICS 

NORAPS,  which  includes  a  regional  analysis  and  forecast  model,  produces  high 
spatial  resolution  (typically  80  km  in  the  western  Atlantic  version),  short-term  (36-48  h) 
numerical  forecasts  over  a  limited  domain.  The  principle  advantage  of  using  a  regional 
model  as  opposed  to  a  global  or  hemispheric  model  is  the  small  spatial  scale  features 
that  are  resolved.  This  resolution  is  very  important  for  simulations  of  maritime 
cyclones,  which  characteristically  are  smaller  than  continental  systems  and  depend  on 
interactions  between  fields  that  varv'  dramatically  over  short  (meso-a  scale)  distances. 
NOR.A.PS  has  the  additional  asset  of  flexibility,  as  the  grid  is  globally  relocatable,  and 
the  user  may  specify  the  forecast  area  dimensions  and  horizontal/vertical  resolution. 
Three  different  projections  (M creator,  Lambert  conformal,  or  polar  stereographic)  are 
available  to  minimize  distortion  in  the  tropics,  mid-latitudes  or  polar  regions.  An 
additional  feature'is  the  "terrain  enveloping"  concept,  in  which  topography  is  calculated 
at  a  high  horizontal  resolution  to  incorporate  the  effects  of  the  sub-grid  scale  features 
into  the  topographic  field.  A  thorough  discussion  of  NORi\PS  is  provided  by  Hodur 
(1982,1984). 

NORAPS  includes  analysis,  initialization,  forecast  and  output  components.  The 
analysis  component  begins  with  acquisition  of  different  types  of  data  (radiosonde, 
satellite,  land  and  ship  reports)  and  quality  control  checks  to  determine  data  validity. 
The  next  step  in  the  analysis  is  the  interpolation  of  the  observations  to  the  model  grid. 
The  data  fields  for  the  model  are  the  u  and  v  wind  components,  temperature, 
geopotential,  specific  humidity,  surface  pressure,  sea-surface  temperature  and  terrain 
height.  A  regional  update  cycle  is  used  in  which  the  12  hour  NORAPS  forecast  serves 
as  the  first  guess.  A  successive  corrections  technique  is  used  for  the  NORAPS 
objective  analysis  to  improve  the  first-guess  fields  of  wind,  temperature,  and 
geopotential.  An  exponential  weighting  function  takes  into  account  the  distance  of  the 
observation  from  the  grid  point  and  the  number  of  observations  surrounding  each  grid 
point.  The  analyses  are  performed  at  eleven  standard  levels  (1000,  925,  850,  700,  500, 
400,  300,  250,  200,  150  and  100  mb). 
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The  initialization  phase  approximately  balances  the  mass  and  wind  fields  and 
suppresses  the  growth  of  large  amplitude  inertial-gravity  waves  that  would  contaminate 
the  forecast  fields.  The  static  initialization  procedure  uses  diagnostic  constraints  to 
relate  the  wind  and  mass  fields. 

The  forecast  component  is  the  heart  of  the  NORAPS  model  and  requires  the 
majority  of  the  computation  time.  The  model  uses  the  flux  form  of  the  primitive 
equations  on  a  staggered  grid  scheme  C  (Arakawa,  1977).  This  grid  scheme  has 
excellent  geostrophic  adjustment  properties  and  group  velocity  characteristics  (Haltiner 
and  Williams,  1980).  The  vertical  coordinate  is  sigma,  which  orients  all  coordinate 
surfaces  parallel  to  the  surface  terrain.  Thus,  the  vertical  velocity  is  identically  zero  at 
the  lower  boundar\',  even  in  the  vicinity  of  mountainous  terrain.  The  vertical  structure 
of  the  atmosphere  is  normally  represented  in  12  discrete  layers,  although  as  many  as  20 
layers  may  be  specified.  All  variables  (u,  v,  q  and  T)  except  vertical  velocity  are  carried 
at  the  middle  of  each  layer.  A  split-explicit  time  integration  scheme  is  used  to  permit 
larger  time  steps  for  the  slower  meteorological  modes,  while  still  being  able  to  predict 
all  the  gravity  modes.  The  size  of  the  time  step  is  governed  by  the  computational 
stability  criterion  for  the  horizontal  resolution  selected.  Fourth-order  advection  is  used 
for  the  prediction  equation  set  to  reduce  errors  in  phase  speed. 

One-way  interactive  boundary  conditions  are  used  to  specify  the  time-dependent 
lateral  boundary  conditions  on  the  NORAPS  grid  from  the  coarser  mesh  NOGAPS 
predictions.  The  NOGAPS  solution  forces  the  fine-mesh  model,  without  the  fine  grid 
affecting  the  coarse  grid  solution.  Operational  time  constraints  require  that  these 
boundary  conditions  be  derived  from  the  most  recent  forecast,  rather  than  the  forecast 
from  the  same  time.  A  method  developed  by  Perkey  and  Kreitzberg  (1976)  is  used  to 
spatially  interpolate  the  solutions  near  the  boundarv'  of  the  finer  mesh.  The  NOGAPS 
time  tendencies  are  blended  with  the  NORAPS  time  tendencies  over  a  distance  of  five 
grid  points  to  dampen  spurious  reflections  due  to  the  change  in  grid  spacing  at  the 
regional  model  boundary'. 

The  output  phase  of  NOR.A.PS  prepares  the  forecast  data  for  interpolation  to 
standard  pressure  levels.  Output  fields  include  winds,  temperature,  specific  humidity, 
surface  pressure,  relative  humidity,  absolute  vorticity,  divergence,  surface  sensible  and 
latent  heat  flux,  and  precipitation. 
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2.        MODEL  PHYSICS 

The  model  physics  are  crucial  for  the  simulation  of  maritime  cyclones.  NORAPS 
includes  representations  of  the  following  physical  processes:  (i)  dry  convective 
adjustment;  (ii)  surface  friction;  (iii)  cumulus  parameterization;  (iv)  large-scale 
precipitation;  and  (v)  radiative  transfer  processes, 

a.  Planetan'  Boundary'  Layer 

The  planetary  boundary  layer  (PBL),  which  is  confined  to  the  bottom  three 
layers  of  the  model  atmosphere,  is  well  mixed  in  temperature,  momentum  and 
moisture.  Interactions  between  the  lower  boundary  and  overlying  atmospheric  layer 
provide  sources  and  sinks  for  momentum,  heat  and  moisture.  The  effects  of  the  PBL 
should  be  included  in  any  numerical  model  to  physically  simulate  maritime  cyclogenesis 
on  the  time  scales  of  more  than  a  few  hours  (Anthes  et  al.,  1983). 

The  NOR.A.PS  PBL  paranieterization  follows  Deardorfr(1972).  After  the  layer 
mean  values  of  V,  0  and  q  are  known,  a  bulk  Richardson  number  (Ri^)  is  computed  to 
determine  the  stability  of  the  PBL, 

Rib=  gh(e^  -  0^  )/U^2^  (B.l) 

'  m        s 

where  0^  is  the  virtual  potential  temperature,  subscript  s  denotes  surface  values  and 
subscript  m  denotes  mean  PBL  values.  For  unstable  conditions  (Ri^<0),  i.e.,  strong 
winds,  daytime  heating  over  land  or  strong  mixing,  a-  predictive  equation  for  the  PBL 
height  (h)  proposed  by  Stull  (1976)  is  used.  The  rate  of  change  of  the  PBL  height  is 
related  to  the  surface  sensible  heat  flux,  mean  PBL  wind  speed,  the  large-scale  vertical 
motion  and  cloud-induced  subsidence.  For  stable  or  neutral  conditions  (Ri|^>0).  a 
predictive  equation  for  the  PBL  height  after  Nieuwstadt  and  Tennekes  (1981)  is  used. 
Transfer  (drag)  coefficients  for  heat  (Cq)  and  friction  (C^)  are  computed  from 
empirical  formulas  that  include  the  stability  dependence  via  the  Ri^^.  Surface  fluxes  of 
moisture,  heat  and  momentum  are  computed  using  the  bulk  aerodynamic  formulas, 
which  assume  the  transfer  coefficients  are  functions  of  Ri^  and  h.  Seasonally 
dependent  climatological  values  of  albedo,  sea  ice,  ground  wetness  and  surface 
roughness  are  specified.  A  predictive  equation  for  the  ground  temperature  after 
Blackadar  (1979)  is  used  to  model  the  lower  boundary  condition  for  the  temperature 
over  land.    The  sea-surface  temperatures  are  assumed  to  be  constant  over  the  forecast 
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period,  which  is  valid  for  short-range  forecasts.  The  PBL  is  constrained  to  extend 
through  at  least  the  bottom  layer  of  the  model  (approximately  40  mb)  to  avoid 
extrapolation  problems  in  determining  mean  PBL  quantities. 

b.  Cumulus  Parameterization 

The  NORAPS  model  uses  a  modified  version  of  the  Kuo  (1965)  cumulus 
parameterization  scheme.  This  version  links  the  convection  to  the  PBL  by  requiring 
moisture  convergence  in  the  PBL.  By  contrast,  the  original  Kuo  version  required  net 
moisture  convergence  in  the  entire  column  before  convection  was  initiated.  The 
moisture  convergence  is 

M,=  1/g  V  .  (qn,7rV^)(l-<Tpbi)  +  P,{^^'q\   ,  '  (B.5) 

where  the  first  term  on  the  right  side  is  the  vertically-integrated  moisture  convergence 
and  the  second  term  is  the  surface  moisture  flux.  Convection  is  assumed  to  occur 
when  M^  >  4.0  x  10  gm  m*-^  s'^  and  the  equivalent  potential  temperature  decreases 
with  height  (conditionally  unstable)  from  the  PBL  to  the  first  model  layer  above  the 
PBL.  The  final  constraint  is  that  deep  convection  cannot  occur  if  the  lifting 
condensation  level  (LCL)  is  above  the  PBL.  -This  scheme  partitions  the  moisture 
transport  into  two  fractions:,  the  first  (bM^.)  serves  to  moisten  the  ^environment  to 
saturation  conditions  through  the  cloud  layer,  and' the  second  ((l-b)M^)  condenses  and 
falls  instantaneously  as  rain.  ,  • 

In  addition  to  large-scale  advection,  temperature  and  moisture  changes  at  any 
level  are  caused  by  convective  clouds.  It  is  assumed  that  the  temperature  of  the  cloud 
is  warmer  than  the  environment.  Cloud  production,  which  is  the  ratio  of  the  water 
vapor  available  to  the  water  vapor  needed  to  form  the  cloud,  is  computed  for  each 
gridpoint.  The  fractional  cloud  area  is  used  to  adjust  the  layer  mean  temperatures  and 
moistures  to  account  for  the  presence  of  clouds. 

c.  Large-Scale  Precipitation 

Large-scale  (non-convective)  precipitation  can  occur  when  supersaturation  is 
achieved  at  any  level.  The  excess  moisture  is  allowed  to  fall  into  the  next  layer  and 
increase  the  moisture  content  of  that  layer,  or  continue  to  fall  if  that  layer  is  already 
supersaturated.  Precipitation  occurs  only  when  the  air  is  saturated  from  the  cloud  to 
the  ground.    Physical  processes  evolve  more  slowly  and  are  not  as  severely  restricted  by 
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time  step  limitations  as  are  the  dynamics.  The  precipitation  routines  are  only  called 
every  eight  time  steps  (once  every  30  minutes)  for  computational  efficiency.  The 
heating  and  moistening  rates  are  distributed  evenly  over  subsequent  time  steps  until  the 
next  call  to  these  routines. 

d.  Radiation 

The  incorporation  of  solar  radiation  into  numerical  models  is  essential  for 
prediction  of  surface  temperatures  and  the  cooling  rates  at  cloud  tops  that  may  deepen 
cloud  layers.  The  radiation  parameterization  in  NORAPS  follows  Katayama  (1974)  for 
short-wave  radiation  and  Sasamori  (1968)  for  long-wave  radiation. 
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APPENDIX  C 
IOP-1  PRIMARY  LOW  BUDGETS 

The  low  pressure  center  which  formed  in  KY  at  1900  (denotes  00  GMT  19 
Januar\'  1986)  during  IOP-1  and  filled  in  western  PA  after  2000  showed  only  minor 
development  compared  to  a  second  low  that  formed  on  the  NC  coast  at  1912  (Chapter 
V).  Selected  figures  from  the  quasi-Lagrangian  diagnostic  budgets  of  the  first  low  are 
presented  for  completeness  in  this  appendix. 

1.  MASS  BUDGET 

In  the  mass  divergence  pattern  for  the  IOP-1  primary  low  (Fig.  C.l), 
convergence  is  limited  to  the  lowest  layers.  In  the  analysis  budget,  the  level  of  non- 
divergence  (LND),  which  is  initially  at  600  mb,  lowers  to  about  650  mb.  and  then  rises 
as  the  primar\'  low  merges  with  the  secondary  system  at  2000.  Divergence  aloft  is 
weaker  than  the  low-level -convergence,  but  extends  through  a- greater  depth.  The 
average  vertical  motion  (not  shown)  in  the  budget  volume  is  upward  at  all  times,  with 
the  peak,  velocity  at  the  LND  (near  680  mb).  There  are  only  small  changes  with  time 
in  the  mass  circulation  after  the  initial  increase  in  boundary'  layer  convergence. 

As  expected  from  the  synoptic  discussion  in  Chapter  V,  the  mass  circulation  in 
the  forecast  (Fig.  C.lb)  is  quite  similar  to  that  in  the  analysis.  Surface  convergence 
and  upper-level  divergence  are  slightly  stronger  in  the  forecast,  but'  the  general  pattern 
is  much  the  same.  The  LND  remains  below  600  mb  throughout  the  period,  but  an 
area  of  convergence  near  400  mb  is  present  at  2000,  as  the  primar\'  system  begins  to  be 
influenced  by  the  secondary'  low.  As  in  the  analysis,  there  is  little  change  in  the  pattern 
or- strength  of  the  circulation  from  the  beginning  to  the  end  of  the  budget. 

2.  VORTICITY  BUDGET 

Time-height  cross-sections  of  i^^  show  that  the  spin-up  at  all  levels  in  this  system 
was  small.  The  analysis  (Fig.  C.2a)  has  a  slight  (^  increase  near  the  tropopause  and  an 
even  smaller  increase  at  lower  levels.  The  forecast  (Fig.  C.2b)  has  only  a  small  increase 
below  700  mb.  In  both  budgets,  a  minimum  in  C,^  near  450  mb  persists  throughout  the 
life  of  the  low. 

The  column-averaged  vorticity  budget  for  the  primar\'  low  (Fig.  C.3)  has  a 
similar  pattern  as  that  for  the  secondary-  low  (Chapter  V),  but  the  magnitudes  oC  all  the 
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Figure  C.l.   Horizontal  mass  flux  across  the  4°  lat.  radius  in  the  (a)  analysis  and  (b) 
forecast  of  the  primar}'  low.   Solid  (dashed)  contours  denote  outflow  from  (inflow  into) 
the  4°  lat.  radius  budget  volume.    Contours  are  5  x  10'^-'  g  s"^  mb'^ 


158 


962- 


1900 


12 


B 

175- 

Fcst 

22S 

275- 

^    350 


oc 

3 
CO 
CO 
UJ 
flC 
CL 


450 


600 


775- 


8SS- 


96> 


1912 


2000 


1900 


1906 


1912 


1918 


2000 


Figure  C.2.   Area-averaged  absolute  vorticity  over  4°  lat.  radius  in  the  (a)  analysis  and 
(b)  forecast  of  the  primary'  low.   Contours  are  2  x  10'^  s'K 
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Figure  C.3.   Average  over  the  lower  troposphere  (1000-500  mb)  of  the  vorticity  budget 
terms  for  the  (a)  analysis  and  (b)  forecast  of  the  primary  low.   Time  tendency  (solid), 
lateral  transport  (dashed),  vorticity  divergence  (dotted),  vertical  transport 
(dash-dotted),  tiltmg  term  (long-short  dashed)  and  residual  (dash-asterisk). 
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terms  are  smaller.  Leading  spin-up  terms  are  the  horizontal  transport  of  ^^  and  the 
divergence  term,  while  the  vertical  flux  convergence  and  the  residual  term  tend  to 
decrease  the  lower  tropo spheric  vorticity  in  both  the  analysis  and  forecast. 

The  vorticity  generation  due  to  the  divergence  term  for  the  analysis  and  forecast 
of  the  primar\'  low  is  presented  in  Fig.  C.4.  With  the  exception  of  a  slight  increase  in 
the  lowest-level  spin-up,  the  analysis  has  only  minor  changes  with  time.  In  the  forecast 
(Fig.  C.4b),  the  upper-level  C^^  decrease  becomes  smaller  as  the  boundary  layer  spin-up 
continues. 

3.  HEAT  BUDGET 

The  relative  importance  of  the  different  terms  in  the  column-averaged  heat 
budgets  for  the  analysis  and  forecast  (Fig.  C.5)  are  similar  to  those  for  the  secondary 
low.  The  temperature  tendency  is  small  during  all  time  periods  in  both  the  analysis 
and  forecast,  which  is  consistent  with  the  slowly-changing  nature  of  the  low.  Adiabatic 
cooling  of  about  7°  C  d'^  in  the  analysis  is  balanced  by  small  amounts  of  warm 
advection  and  diabatic  heating.  In  the  forecast  budget,  the  magnitudes  of  the  adiabatic 
cooling  and  diabatic  heating  terms  are  increased  to  near  10°  C  d'^  This  adiabatic 
cooling  reflects  the  stronger  vertical  circulation  in  the  forecast  that  was  identified  in  the 
mass  budget. 

The  residual  (inferred  diabatic  heating)  in  the  primary  low  is  nearly  constant  in 
time  (Fig.  C.6).  The  level  of  maximum  heating  is  near  650  mb  in  both  the  analysis  and 
forecast.  The  secondary'  maximum  in  the  stratosphere,  which  is  evident  at  1906  in  the 
analysis  and  at  1915  in  the  forecast,  is  probai?iy  attributable  to  accumulated  errors. 
The  diabatic  heating  in  the  forecast  (Fig.  C.6b)  has  the  same  pattern  as  in  the 
analysis,  except  the  magnitudes  are  consistently  5  °  C  d"^  larger. 

4.  MOISTURE  BUDGET 

The  last  budget  for  the  lOP-l'  lows  is  the  moisture  budget,  which  can  only  be 
calculated  for  the  forecast  fields.  The  column-averaged  moisture  budget  has  a  pattern 
similar  to  that  in  the  secondar.'  low  (Chapter  V),  with  the  horizontal  flux  providing 
moisture  to  the  lower  troposphere,  and  the  moisture  sink  (precipitation)  term  providing 
the  greatest  amount  of  drying  (Fig.  C.7a).  The  moisture  tendency  is  positive  in  the 
early  period,  as  the  relatively  large  horizontal  flux  provides  more  moisture  than  is 
removed  by  condensation  and  the  vertical  flux.  As  the  horizontal  flux  decreases 
toward  the  end  of  the  budget  calculation,  the  tendency  becomes  negative  and  the 
budget  volume  begins  to  dr\'. 
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As  in  the  secondary  low,  the  distribution  of  condensation  in  the  vertical 
calculated  as  a  budget  residual  coincides  with  the  maximum  drying  as  diagnosed  in  the 
forecast  heat  budget.  The  largest  magnitude  in  the  moisture  sink  term  is  centered  near 
650  mb  (Fig.  C.7b). 
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Figure  C.4.   Vorticity  divergence  term  within  4°  lat.  radius  in  the  (a)  analysis  and  (b) 
forecast  of  the  primary  low.   Solid  (dashed)  contours  indicate  vorticitv  increase 
(decrease).   Contours  are  5  x  10"^°  s'-^. 
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Figure  C.5.  Average  over  the  lower  troposphere  (1000-500  mb)  of  the  heat  budget 
terms  for  the  (a)  analysis  and  (b)  forecast  of  the  primary'  low.   Time  tendency  (solid), 
horizontal  advection  (dashed),  adiabatic  cooling  (dotted)  and  diabatic  heating 
(dot-dashed). 
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Figure  C.7.   (a)  Average  over  the  lower  troposphere  (1000-500  mb)  of  the  moisture 
budget  terms  for  the  forecast  of  the  primary'  low.   Time  tendency  (solid),  moisture  sink 
(dashed),  horizontal  flux  (dotted)  and  vertical  flux  (dot-dashed),  (b)  Moisture  sink 
within  4^^  lat.  radius  for  the  primary'  low  forecast.    Dashed  contours  indicate 
condensation  in  g  kg'^  d"  . 
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